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Abstract 

This  dissertation  explores  optimal  path  planning  for  air  vehicles.  An  air  vehicle  exposed  to  illumination  by 
a  tracking  radar  is  considered  and  the  problem  of  determining  an  optimal  planar  trajectory  connecting  two 
prespecified  points  is  addressed.  An  analytic  solution  yielding  the  trajectory  minimizing  the  received  radar 
energy  reflected  from  the  target  is  derived  using  the  Calculus  of  Variations.  Additionally,  the  related 
problem  of  an  air  vehicle  tracked  by  a  passive  sensor  is  also  solved.  Using  the  insights  gained  from  the 
single  air  vehicle  radar  exposure  minimization  problem,  a  hierarchical  cooperative  control  law  is 
formulated  to  determine  the  optimal  trajectories  that  minimize  the  cumulative  exposure  of  multiple  air 
vehicles  during  a  rendezvous  maneuver.  The  problem  of  one  air  vehicle  minimizing  exposure  to  multiple 
radars  is  also  addressed  using  a  variational  approach,  as  well  as  a  sub-optimal  minmax  argument.  Local 
and  global  optimality  issues  are  explored.  A  novel  decision  criterion  is  developed  determining  the 
geometric  conditions  dictating  when  it  is  preferable  to  go  between  or  around  two  radars.  Lastly,  an 
optimal  minimum  time  control  law  is  obtained  for  the  search  and  target  identification  mission  of  an 
autonomous  air  vehicle.  This  work  demonstrates  that  an  awareness  of  the  consequences  of  embracing 
sub-optimal  and  non-globally  optimal  solutions  for  optimization  problems,  such  as  air  vehicle  path 
planning,  is  essential. 
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Abstract 

This  dissertation  explores  optimal  path  planning  for  air  vehicles.  An  air  vehicle 
exposed  to  illumination  by  a  tracking  radar  is  considered  and  the  problem  of  determining 
an  optimal  planar  trajectory  connecting  two  prespecified  points  is  addressed.  An  analytic 
solution  yielding  the  trajectory  minimizing  the  received  radar  energy  reflected  from  the 
target  is  derived  using  the  Calculus  of  Variations.  Additionally,  the  related  problem  of  an 
air  vehicle  tracked  by  a  passive  sensor  is  also  solved. 

Using  the  insights  gained  from  the  single  air  vehicle  radar  exposure  minimization 
problem,  a  hierarchical  cooperative  control  law  is  formulated  to  determine  the  optimal 
trajectories  that  minimize  the  cumulative  exposure  of  multiple  air  vehicles  during  a  ren¬ 
dezvous  maneuver.  The  problem  of  one  air  vehicle  minimizing  exposure  to  multiple  radars 
is  also  addressed  using  a  variational  approach,  as  well  as  a  sub-optimal  minmax  argument. 
Local  and  global  optimality  issues  are  explored.  A  novel  decision  criterion  is  developed  de¬ 
termining  the  geometric  conditions  dictating  when  it  is  preferable  to  go  between  or  around 
two  radars.  Lastly,  an  optimal  minimum  time  control  law  is  obtained  for  the  target  iden¬ 
tification  and  classification  mission  of  an  autonomous  air  vehicle. 

This  work  demonstrates  that  an  awareness  of  the  consequences  of  embracing  sub- 
optimal  and  non-globally  optimal  solutions  for  optimization  problems,  such  as  air  vehicle 
path  planning,  is  essential. 
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AIR  VEHICLE  PATH  PLANNING 


I.  Introduction 

1.1  Overview 

In  response  to  the  lessons  learned  from  recent  conflicts,  the  Air  Force  is  expanding  the 
uses  and  roles  of  the  Uninhabited  Air  Vehicle  (UAV).  Traditionally  used  for  reconnaissance, 
UAVs  are  now  being  considered  for  combat  operations,  including  autonomous  attack.  The 
Uninhabited  Combat  Air  Vehicle  (UCAV)  is  envisioned  to  carry  out  the  “dull,  dirty  and 
dangerous”  missions  of  suppression  of  enemy  air  defenses  and  the  prosecution  of  time 
critical  targets;  missions  that  often  involve  high  risk  to  pilots  and  demand  high  endurance 
[25].  Similarly,  air  weapons  designers  are  researching  standoff  miniature  munitions  that 
autonomously  search,  detect,  identify,  attack  and  destroy  time  critical  and  ground  mobile 
targets  of  military  significance  [21]. 

The  intention  to  develop  a  significant  military  force  consisting  of  UAVs,  UCAVs  and 
autonomous  munitions  provides  both  opportunities  and  imperatives  for  research.  If  these 
uninhabited  air  vehieles  are  to  make  effeetive  eontributions  to  air  power,  a  progression 
from  single  vehiele  autonomy  to  multiple  vehiele  eooperative  behavior  is  required.  The 
purpose  of  this  research  is  to  contribute  to  the  development  of  a  theory  of  path  planning 
for  autonomous  air  vehicles. 

1.2  Historieal  Overview 

1.2.1  Single  Vehiele  Path  Planning.  The  single  vehicle  path  planning  problem 
entails  finding  a  path  from  point  A  to  point  B,  such  that  one  minimizes  or  avoids  risk, 
while  subject  to  constraints,  e.g.,  limited  fuel,  minimum  time,  etc.,  see,  e.g..  Fig.  1.1. 

The  task  of  planning  optimal  and  sub-optimal  paths  for  single  vehicle  operations 
has  been  addressed  by  the  disciplines  of  Control  Engineering,  Operations  Research  and 
Computer  Science.  In  general,  path  planning  is  approached  using  four  methods:  trajectory 
optimization,  route  planning,  analogy-based  path  planning  and  game  theoretic  methods. 
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Figure  1.1  Simple  Path  Planning  Problem 

1.2. 1.1  Trajectory  Optimization.  One  approach  used  to  solve  air  vehicle 
path  planning  problems  is  trajectory  optimization.  Trajectory  optimization  seeks  time 
histories  of  state  and/or  control  variables  of  a  dynamic  model  which  minimizes  a  particular 
cost  function  or  performance  index.  Closed  form  solutions  of  these  time  histories  are  rarely 
obtained.  Thus,  numerical  techniques  are  usually  employed. 

Trajectory  optimization  problems  are  generally  posed  as  continuous  time  non-linear 
optimization  problems  that  are  addressed  using  non-linear  parameter  optimization,  Calcu¬ 
lus  of  Variations,  optimal  control  theory  or  dynamic  programming.  Discrete  time  and/or 
linear  approximations  are  commonly  made.  Many  classic  trajectory  optimization  problems 
have  been  solved  using  these  tools  and  techniques,  such  as  the  brachristichrone  problem, 
the  isoperimetric  problem,  numerous  minimum-time  performance  problems,  Zermelo’s  nav¬ 
igation  problem  and  others  -  see,  e.g.,  [16,  8,  7]. 

Numerical  codes  have  been  developed  for  trajectory  optimization  purposes.  Recently, 
Bryson  [7]  published  an  extensive  set  of  methods  and  codes  for  trajectory  optimization. 
Furthermore,  commercial  software  packages,  such  as  the  MATLAB  Optimization  Toolbox 
[13],  and  software  libraries  such  as  Numerical  Recipes  [36]  and  IMSL  [43],  contain  many 
tested  and  widely  used  routines.  Of  particular  interest  in  this  work  are  shooting  methods 
for  two-point  boundary  value  problems.  Roberts  and  Shipman  [40]  have  authored  a  classic 
text  on  the  subject.  FORTRAN  implementations  of  the  shooting  method  are  available 
from  many  sources  [11,  43,  36]. 
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1.2. 1.2  Route  Planning.  Generally,  in  operational  contexts,  route  planning 
is  performed  rather  than  trajectory  optimization.  A  route  consists  of  a  two  or  more  way- 
points  connected,  usually,  by  straight  lines.  Often  these  problems  employ  a  discretization 
or  gridding  scheme,  where  the  space  of  feasible  waypoints  is  predefined  and  in  general, 
bounded. 

A  subset  of  route  planning  problems  involve  the  situation  where  a  set  of  waypoints 
is  provided  and  one  must  determine  the  optimal  visitation  sequence  that  minimizes  some 
cost  function.  All  but  the  simplest  route  planning  problems  are  solved  numerically  via 
dynamic  programming,  network/graph  theory  or  combinatorial  optimization.  Efficient 
search  methods  and  heuristics  are  employed  to  make  the  computational  burden  reasonable. 
Unless  complete  enumeration  is  feasible,  optimal  solutions  cannot  be  guaranteed. 

1.2. 1.3  Analogy-based  Path  Planning.  Analogy-based  path  planning  in¬ 
volves  transforming  the  path  planning  problem  to  an  entirely  different  problem  that  has 
already  been  solved  or  to  one  that  has  a  convenient  method  of  solution.  Transformations 
can  be  found  from  a  variety  of  concepts  such  as  physics,  geometry  and  biology. 

Techniques  developed  for  robotic  motion  planning  include  a  method  of  path  planning 
using  potential  fields.  This  technique  has  been  used  by  several  researchers  [27,  17,  5]  at  Air 
Force  Research  Laboratory  (AFRL)  investigating  optimal  path  planning  for  autonomous 
air  vehicles.  The  AFRL  technique  relies  upon  a  defined  grid  and  a  system  of  connections 
between  nodes  on  the  grid,  based  upon  physical  laws  such  as  the  attraction  and  repulsion  of 
magnetic  fields.  Bortoff  [5]  employed  a  potential  field  method  by  representing  the  path  of  a 
UAV  through  a  radar  threat  zone  as  a  chain  of  spring-mass  systems.  The  radars  generated 
virtual  forces  proportional  to  the  l/R'^  distance  law  and  “pushed”  the  spring  mass  system 
away.  The  chain  is  initialized  in  some  arbitrary  configuration  and  through  simulation, 
reaches  it’s  potential  energy  minimum,  which  is  a  weighted  sum  of  the  path  length  and  the 
distance  from  the  radar.  Potential  fields,  however  popular,  have  shortcomings.  Koren  and 
Borenstien  [24]  identified  four  significant  problems  with  the  potential  field  approach:  trap 
situations  due  to  local  minima,  no  passage  through  closely  spaced  obstacles,  oscillations  in 
the  presence  of  obstacles,  and  oscillations  in  narrow  passageways. 
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Helgason,  et  al  [20],  developed  a  geometric  algorithm  that  finds  a  path  through  threat 
regions  modelled  as  circles.  The  routes  were  selected  by  moving  between  waypoints  in  a 
fashion  tangential  to  the  threat  circles.  The  waypoints  are  identified  by  circumscribing 
various  triangles  around  the  threat  circles  in  a  branch- and-bound  scheme.  In  a  similar 
manner,  Asseo  [4]  developed  an  algorithm  to  avoid  threat  zones  represented  as  circles. 
Asseo  generates  routes  consisting  of  linear  segments  tangent  to  the  threat  circles,  however 
Asseo  also  admits  circular  segments  as  well,  travelling  along  the  threat  circles  for  some 
segments.  Compared  to  the  grid  search  algorithms,  these  geometric  approaches  impose  a 
minor  computational  burden,  yet  they  make  no  claims  of  optimality. 

1.2.2  Game  Theoretic  Path  Planning.  Contained  mainly  to  applications  involv¬ 
ing  economics  and  warfare,  game  theory  emerged  in  the  late  1940s  from  the  fields  of  Applied 
Mathematics  and  Operations  Research.  Games  of  pursuit  and  evasion  are  commonly  stud¬ 
ied  using  the  methods  of  game  theory.  In  general,  these  games  consist  of  the  pursuer  P 
attempting  to  capture  the  evader  E  by  closing  the  distance  PE  to  within  some  prescribed 
positive  quantity.  For  example,  the  pursuer  could  be  a  torpedo,  a  fighter  aircraft  or  a 
guided  missile,  and  the  evader  could  be  a  warship,  a  bomber  aircraft  or  a  tank. 

The  simplest  games  often  involve  one  pursuer  and  one  evader,  and  the  dynamics  of 
their  motion  is  governed  by  differential  equations  with  bounds  on  speed  and  rates  of  turn. 
Consider  the  game  of  the  homicidal  chauffeur,  see,  e.g.,  [22].  The  pursuer  (the  chauffeur) 
and  his  quarry  are  restricted  to  movement  in  the  plane.  The  pursuer  moves  at  a  fixed 
speed,  wi,  with  a  bounded  turning  radius,  R.  The  evader  moves  at  a  constant  speed, 
W2  <  wi,  however,  the  evader  can  change  direction  at  any  moment.  Capture  occurs  when 
the  pursuer  and  evader  are  separated  by  a  distance  PE  <  1. 

Two  forms  of  the  homicidal  chauffeur  game  are  typically  studied.  In  the  game  of  kind, 
we  are  concerned  with  determining  the  conditions  under  which  the  pursuer  can  capture  the 
evader,  viz.,  what  values  of  the  parameters  wi/w2,R  and  I  guarantee  that  P  can  capture 
E?  Alternatively,  assume  the  game  is  designed  such  that  P  can  always  capture  E.  When 
we  are  concerned  with  optimizing  the  time  of  capture,  we  are  playing  a  game  of  degree. 
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Here,  the  evader  hopes  to  maximize  the  time  before  capture,  while  the  pursuer  is  working 
to  minimize  this  quantity. 


\ 


Figure  1.2  Isaac’s  Illustration  Concerning  the  Game  of  the  Homicidal  Chauffeur 

The  game  of  the  homicidal  chauffeur  has  many  interesting  scenarios  arising  from 
different  initial  starting  conditions.  For  example,  consider  Fig.  1.2,  from  Isaacs  [22].  At 
(a)  of  Fig.  1.2  we  have  the  pursuer  travelling  initially  up,  with  the  evader  ahead  and  to 
the  right.  The  optimal  strategy  is  for  P  to  turn  hard  to  the  right,  following  the  curve  of  a 
minimum  radius  turn,  until  the  pursuer’s  velocity  vector  is  aligned  with  E.  The  game  ends 
with  capture  of  E  at  the  point  C.  Now  consider  Fig.  1.2  at  (6),  where  P  is  still  proceeding 
initially  up,  however  E  is  to  the  rear  of  P.  If  the  pursuer  makes  the  same  hard  turn  to 
the  right,  E  has  an  opportunity  to  get  inside  the  pursuer’s  turning  radius,  thus  allowing 
E  to  avoid  capture.  To  capture  E,  the  pursuer  must  first  make  an  indirect  move.  In  (c), 
a  “swerve”  maneuver  is  depicted,  whereby  P  initially  turns  hard  to  the  left,  away  from  E, 
until  a  subsequent  hard  turn  to  the  right  will  result  in  the  advantageous  situation  of  (a). 

More  complicated  games  can  involve  multiple  pursuers  and  evaders,  such  as  the 
football  problem:  several  tacklers  against  a  ball  carrier,  or  the  case  of  a  few  fighter  aircraft 
against  a  fleet  of  bombers.  For  more  information  on  game  theoretic  path  planning,  the 
reader  is  referred  to  the  text  by  Isaacs  [22].  Isaacs  highlights  many  different  types  of  two 
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player  games  of  pursuit  and  evasion,  and  these  results  lend  themselves  naturally  to  the 
problem  of  path  planning  for  air  vehicles. 

1.2.3  Multiple  Vehicle  Control.  Path  planning  for  multiple  vehicles  will  use  one  or 
more  of  the  path  planning  methods  previously  described  for  single  vehicles.  Multiple  vehicle 
path  planning  can  be  organized  into  two  distinct  frameworks:  hierarchical  cooperative 
control  and  reactive  control. 

1.2. 3.1  Hierarchical  Cooperative  Control.  Hierarchical  control  systems  are 
typified  by  a  segmentation  of  functionality,  where  the  layers  of  control  coordinate  with 
each  other  in  a  predicable  and  deliberate  fashion.  Decisions  can  be  made  by  a  centralized 
leader  or  by  a  team  of  decision  makers  in  a  decentralized  manner. 

Centralized  hierarchical  control  can  be  achieved  through  encapsulation  or  recursion. 
Much  work  in  this  area  of  is  focused  on  formation  flight  control,  e.g.,  reducing  pilot  work¬ 
load  on  Special  Operations  Forces  missions  [9]  or  minimizing  fuel  expenditure  for  cargo 
aircraft.  Freund  and  Hoyer  [15]  describe  a  method  for  controlling  multiple  robots  where  the 
dynamic  equations  for  individual  robots  are  coupled  and  fed  to  a  hierarchical  coordinator 
which  handles  the  problems  of  collision  avoidance,  obstacle  avoidance  and  path  planning  - 
see,  e.g..  Fig.  1.3. 

McLain  [28]  details  an  approach  to  coordinating  the  rendezvous  of  multiple  UAVs 
using  a  decomposition  strategy.  This  decentralized  approach  was  driven  by  a  state  machine 
residing  on  each  UAV  whereupon  an  identical  “team-optimal”  estimated  time  of  arrival 
(ETA)  solution  is  determined  by  each  UAV.  Once  an  ETA  is  determined  and  coordinated, 
the  individual  UAVs  solve  local  optimization  problems  to  plan  paths  to  meet  the  ETA 
goal,  as  well  as  maximize  their  own  survivability. 

Hierarchical  control  schemes  often  attempt  to  take  advantage  of  high-level  knowledge- 
based  representations  or  maps  of  the  environment.  Typically  these  maps  are  generated  a 
priori  and  can  be  updated  with  on-line  sensor  information.  Paths  are  planned  based  upon 
the  knowledge  that  is  available  to  the  air  vehicles,  and  robustness  to  uncertainty  or  partial 
knowledge  is  an  active  area  of  research. 
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Figure  1.3  Freund  and  Hoyer’s  Hierarchical  Control  Structure 

1.2. 3. 2  Robotic  Navigation  and  Reactive  Control.  The  problem  of  safely 
navigating  robots  through  a  cluttered  environment  is  closely  related  to  the  problem  of 
routing  autonomous  aircraft.  Although  methods  involving  hierarchical  planning  in  robotic 
navigation  have  been  successful,  dissatisfaction  with  their  performance  led  to  the  emergence 
of  the  reactive  control  approach.  Influenced  by  the  Behaviorist  school  of  psychology,  Brooks 
developed  the  subsumption  architecture  [6] ,  in  which  a  model  or  representation  of  the  world 
is  eschewed  in  favor  of  a  system  of  simple  behaviors  (like  wander,  avoid  obstacles,  move  to 
goal).  In  the  subsumption  architecture,  the  intended  behavior  of  the  robot  emerges  from 
the  interaction  of  the  simple  behaviors  with  the  environment.  Reactive  control  research  is 
concerned  with  how  best  to  coordinate  these  simple  behaviors  to  achieve  desired  emergent 
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behaviors.  In  principle,  the  subsumptive  architecture  promotes  quick  response  to  changes 
in  the  world  and  hence,  adaptability  to  unplanned  circumstances. 

Arkin,  Ram  and  others  have  made  numerous  contributions  [39,  38,  37]  to  the  area 
of  reactive  control  for  robotic  navigation,  addressing  some  of  the  shortcomings  [24]  of  the 
purely  reactive  control  model.  Much  of  the  work  has  focused  on  how  the  robots  adapt  their 
behavioral  structure  and  control  parameters  using  machine  learning  techniques.  The  work 
on  reactive  control  of  robots  has  been  extended  to  include  cooperative  behaviors  amongst 
the  robots.  Arkin  and  Balch  give  a  short  history  of  the  development  of  multi-agent  robotic 
systems  [3].  Currently,  the  most  popular  application  of  multiple  robot  research  is  the 
RoboCup  robotic  soccer  initiative  [41]. 

Although  the  two  approaches  of  hierarchical  planning  and  reactive  control  seem  or¬ 
thogonal  to  each  other  in  terms  of  their  philosophy,  Arkin  eloquently  describes  the  need 
to  synthesize  the  two  approaches  in  “Towards  the  Unification  of  Navigation  Planning  and 
Reactive  Control”  [2].  He  states,  “The  reality  is  that  the  hierarchical  school  is  involved 
with  plan  formulation  while  the  reactive  school  copes  with  plan  execution  and  the  issues 
of  dynamic  replanning.”  Suggesting  that  the  two  approaches  are  addressing  two  different 
problems,  Arkin  suggests,  “we  look  not  towards  the  superiority  of  one  navigational  method¬ 
ology  over  the  other,  but  rather  toward  a  synthetic,  integrative  approach  that  applies  both 
of  these  paradigms  to  the  issues  of  navigational  planning...” 

1.3  Problem  Statement 

Previous  work  in  air  vehicle  path  planning  has  focussed  on  ad-hoc  or  heuristic  opti¬ 
mization  methods.  Thus,  extending  existing  path  planning  methods  to  multiple  vehicles, 
while  seeking  insight  into  the  nature  of  optimal  solutions,  is  difficult.  The  development  of 
a  comprehensive  theory  of  cooperative  control  for  air  vehicles  is  a  vast  problem,  beyond 
the  scope  of  a  single  dissertation.  A  realistic,  yet  challenging,  problem  statement  for  this 
dissertation  is  as  follows: 

The  objective  is  to  develop  a  mathematically  rigorous  approach  to  a  problem  in 
air  vehicle  path  planning,  seeking  fundamental  truths  concerning  autonomous 
air  vehicles  and  their  cooperative  control. 
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It  is  hoped  that  these  fundamental  truths,  will  guide  planners  in  the  development  of  tools 
for  more  complex  air  vehicle  path  planning  scenarios.  To  this  end,  the  scenarios  of  mini¬ 
mizing  exposure  to  a  threat  radar  and  target  classification  are  selected  as  the  basis  for  this 
research. 

1.4  Key  Results 

The  analysis  conducted  in  this  research  has  led  to  the  following  key  results: 

•  A  closed  form  solution  to  the  single  vehicle  radar  exposure  minimization  problem. 

—  The  solution  is  shown  to  satisfy  necessary  and  sufficient  conditions  for  a  weak 
local  minimum  as  well  as  the  Weierstrass  necessary  conditions  for  a  strong  local 
minimum. 

—  The  solution  is  shown  to  exist  if,  and  only  if,  the  angle  included  between  the 
departure  and  destination  points  is  less  than  60°. 

—  A  closed  form  expression  for  the  optimal  cost  is  obtained. 

—  An  expression  for  the  optimal  path  length,  not  available  in  closed  form,  is  de¬ 
rived  using  elliptic  integrals.  This  formulation  allows  for  efficient  and  accurate 
computation  of  the  optimal  path  length. 

—  A  method  for  determining  optimal  path  length  constrained  solutions  to  the 
single  vehicle  radar  exposure  minimization  problem  is  presented  and  conditions 
for  the  existence  of  these  solutions  are  identified. 

•  A  similar  analysis  is  conducted  for  the  single  vehicle  passive  sensor  exposure  mini¬ 
mization  problem  to  include  a  closed  form  solution  and  closed  form  expressions  for 
the  optimal  path  length  and  cost.  A  comparison  between  the  solutions  to  the  passive 
sensor  and  radar  exposure  minimization  problem  is  made. 

•  It  is  verified  that  the  Voronoi  edge  is  the  locally  optimal  solution  for  exposure  mini¬ 
mization  against  two  radars.  The  issues  of  global  versus  local  optimality  are  explored 
yielding  analytic  results  identifying  the  conditions  where  going  around  (between)  two 
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radars  is  preferable  to  going  between  (around)  two  radars.  A  suboptimal  algorithm 
for  n-radar  exposure  minimization  is  developed  and  examined. 

•  A  hierarchical  cooperative  control  algorithm  is  formulated  to  determine  optimal 
trajectories  minimizing  radar  exposure  for  two  (or  more)  air  vehicles  performing 
isochronous  rendezvous.  The  constructive  nature  of  the  algorithm  renders  the  exis¬ 
tence  of  these  solutions  readily  verifiable. 

•  The  automatic  target  recognition  process  is  modelled  and  a  novel  optimization  prob¬ 
lem  is  formulated  and  solved.  Minimum  time  trajectories  for  air  vehicles  with  a 
minimum  turning  radius  are  employed  to  integrate  these  concepts  and  a  solution  to 
the  autonomous  air  vehicle  search  problem  is  presented. 

1.5  Organization 

This  dissertation  is  organized  as  follows:  The  essential  elements  of  radar  technology 
and  optimization  theory  utilized  in  this  research  are  presented  in  Chapter  2.  In  Chapter  3, 
the  single  vehicle  radar  exposure  minimization  problem  is  addressed,  including  derivations 
of  the  optimal  trajectory,  its  cost  and  path  length.  Next,  an  analysis  is  conducted  for  the 
problem  of  minimizing  exposure  to  a  passive  sensor  in  Chapter  4.  Similarly,  closed  form 
solutions  for  the  optimal  trajectory,  cost  and  path  length  are  derived.  Chapter  5  contains  a 
description  of  the  numerical  methods  developed  and  employed  in  this  research,  including  a 
novel  discrete  formulation  of  the  problem,  as  well  as  an  exposition  on  the  shooting  method 
for  solving  two  point  boundary  value  problems.  Chapter  6  describes  an  extension  of  the 
results  of  Chapter  2  to  address  several  new  problems.  A  hierarchical  cooperative  control 
algorithm  is  formulated  for  the  problem  of  isochronous  rendezvous  of  multiple  air  vehicles. 
Additionally,  the  problem  of  a  single  vehicle  minimizing  exposure  to  two  radars  is  studied 
and  novel  analytic  results  are  derived.  A  suboptimal  technique  for  minimizing  exposure 
to  n-radars  is  presented.  In  Chapter  7,  the  process  of  automatic  target  recognition  is 
modelled  and  an  optimization  problem  is  suggested.  An  optimal  minimum  time  control 
law  is  developed  for  the  problem  of  autonomous  target  classification.  Chapter  8  provides  a 
discussion  of  these  results  and  recommendations  for  further  research.  Lastly,  appendices, 
a  succinct  list  of  relevant  references  and  a  vita  are  provided. 
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1.6  Summary 

The  forthcoming  application  of  UAVs  and  UCAVs  to  the  battlefield  motivated  the 
need  for  path  planning  tools  that  consider  autonomous  vehicles  and  the  cooperative  control 
of  multiple  air  vehicles.  The  multidisciplinary  nature  of  path  planning  has  been  illustrated 
and  the  different  disciplines  and  approaches  to  path  planning  were  discussed.  The  problem 
statement  motivated  the  need  for  mathematical  rigor  and  the  key  results  obtained  through 
this  research  were  summarized. 

In  the  next  chapter,  some  basic  theory  of  radar  technology  is  presented,  setting 
the  stage  for  the  development  of  the  cost  functionals  used  in  this  research.  Additionally, 
requisite  elements  of  the  Calculus  of  Variations  and  optimization  theory  will  be  covered. 
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II.  Underlying  Topics 


2. 1  Introduction 

The  dissertation  research  is  concerned  with  the  theoretical  nature  of  path  planning 
problems  for  air  vehicles.  This  section  provides  the  basic  theory  that  underlies  the  prob¬ 
lem.  Since  the  research  emphasizes  analytic  solutions,  only  the  basic  science  behind  radar 
technology  is  considered.  More  detailed  treatments  of  radar  systems  can  be  found  in  the 
literature  [42].  Most  of  the  analytic  work  will  involve  optimization  and  the  essential  ele¬ 
ments  of  the  Calculus  of  Variations  are  presented.  For  complete  treatments  of  the  Calculus 
of  Variations  and  Optimal  Control,  the  reader  is  referred  to  [16,  26,  18,  8]. 

2.2  Radar 

Radar,  or  “radio  detection  and  ranging” ,  is  one  of  the  primary  sensors  used  by  hostile 
forces  against  air  vehicles.  In  the  simplest  terms,  a  radar  transmits  an  electromagnetic 
signal  and  then  “listens”  for  an  echo  signal  reflected  from  the  target.  By  sensing  the  time 
delay  between  the  transmitted  pulse  and  the  reflected  echo,  the  range  to  the  target  can  be 
determined.  The  range  to  a  target  is  given  by 


where  c  is  the  velocity  of  propagation  of  the  radar  signal  and  At  is  the  time  elapsed  from 
transmission  to  reception  of  the  echo  signal.  In  general,  the  target  azimuth  is  detected 
by  rotating  the  radar  antenna  and  sensing  the  location  of  the  antenna  when  an  echo  is 
received.  The  radial  velocity  of  the  target  can  be  estimated  by  further  processing  the 
received  signal  to  sense  the  Doppler  frequency  shift. 

The  typical  radar  configuration  consists  of  a  collocated  transmitter  and  receiver. 
This  is  known  as  a  monostatic  radar  configuration.  When  the  receiver  and  transmitter  are 
separated  geographically,  this  is  known  as  a  bistatic  radar  configuration.  Bistatic  radar 
configurations,  while  important,  will  not  be  explicitly  considered  in  this  dissertation. 
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2.2.1  The  Radar  Equation.  The  basic  relationship  that  determines  the  maximum 
range  at  which  a  radar  can  detect  a  given  target  is  known  as  the  radar  range  equation. 
Consider  the  radar  transmitter  to  have  an  omnidirectional  antenna.  The  power  of  the 
signal  radiating  from  the  antenna  Pt  would  be  uniformly  distributed  about  a  sphere  with 
surface  area  dvrr^.  The  power  density  of  the  radar  is  given  by  dividing  the  signal  power  by 
the  surface  area.  Since  radars  use  a  directional  antenna,  we  can  adjust  the  signal  power 
by  a  gain  factor  Gt.  Thus,  the  power  density  at  the  target  is  given  by 

,  ,  PtGt 

power  density  at  the  target  =  - — — 5- 


where  R  is  the  range  from  the  radar  to  the  target. 

As  the  echo  propagates  away  from  the  target  and  toward  the  radar  receiving  antenna, 
the  power  density  of  the  echo  decays  with  the  same  At^R?  factor,  assuming  a  monostatic 
radar  configuration.  Thus  given  an  effective  area  Ag  of  the  receiving  antenna,  and  the 
radar  cross  section  of  the  target  a,  the  echo  power  at  the  radar  receiver  is  given  by 


PiGtcAg 

(47ri?2)^ 


(2.1) 


The  “radar  equation” ,  Eq.  (2.1),  can  be  written  as  a  ratio  of  the  received  radar  power 
to  the  transmitted  power, 

Pr  _  GtcrAg 

Pt  ~  {Anf  R4 

and  is  a  crucial  factor  in  the  detection  of  targets  by  a  radar  system  [42] . 


2.3  Optimization 

This  section  covers  the  optimization  theory  that  will  be  used  in  the  dissertation 
research.  Analytic  solutions  will  be  sought,  and  the  two  main  optimization  tools  to  be 
used  are  non-linear  programming  and  the  Calculus  of  Variations. 

2.3.1  Calculus  of  Variations.  The  Calculus  of  Variations  is  essentially  concerned 
with  finding  the  extremals  of  a  functional  and  classifying  them  as  minima  or  maxima. 
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Following  the  development  in  [16],  we  define  a  functional  as  a  correspondence  that  assigns  a 
definite  real  number  to  each  function  (or  curve)  belonging  to  some  class.  Thus,  a  functional 
can  be  thought  of  as  a  function  where  the  independent  variable  is  itself  a  function. 

We  concern  ourselves  with  functionals  of  the  form 

r-b 

J[y\=  j  F{x,y,y')dx,  y{a)  =  A,  y{b)  =  B  (2.2) 

J  a 

where  J[y\  is  defined  on  some  normed  linear  space.  We  define  the  increment  of  a  functional 
as 

AJ[h]  =  J[y  +  h]-  J[y] 

where  h  =  h{x)  is  the  increment  of  the  independent  variable  y  =  y{x).  If  y  is  fixed,  then 
AJ[h]  is  itself  a  functional.  In  general,  AJ[h]  is  non-linear.  Suppose  we  express  AJ[h]  as 

AJ[h]  =  ip[h]  -I-  e  \\h\\ 

where  ip[h]  is  a  linear  functional  and  e  ^  0  as  ||/i||  ^  0.  The  functional  J[y]  is  differentiable 
and  the  linear  functional  ip[h]  is  called  the  variation  of  J[h]  and  is  denoted  by  6J[h]. 

Theorem  2.3.1.  A  necessary  condition  for  the  differentiable  functional  J[y\  to  have  an 
extremum  for  y  =  y  is  that  its  variation  vanish  for  y  =  y,  i.e.,  that 

6J[y]  =  0 


for  y  =  y  and  all  admissible  h. 

Proof.  See  Theorem  2,  Section  3  in  [16].  □ 


2. 3. 1.1  The  Euler  Equation.  We  now  consider  one  of  the  most  well  known 
results  of  the  Calculus  of  Variations,  the  Euler  equation. 

Theorem  2.3.2.  Let  .J[y]  be  a  functional  of  the  form 


F{x,y,y')  dx 
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defined  on  the  set  of  funetions  y{x)  whieh  have  eontinuous  first  derivatives  in  [a,b]  and 
satisfy  the  boundary  eonditions  y{a)  =  A,y{b)  =  B.  Then  a  neeessary  eondition  for  J[y\ 
to  have  an  extremum  for  a  given  funetion  y{x)  is  that  y{x)  satisfy  Euler’s  equation 


(2.3) 


Proof.  See  Section  4.1  in  [16].  □ 

The  solution  of  Euler’s  equation  in  general  involves  solving  a  second  order  differential 
equation.  There  are  special  cases  where  Euler’s  equation  can  be  reduced  to  a  first  order 
differential  equation  or  another  simplified  expression. 

Case  1:  Assume  the  functional  is  of  the  form 


/  F{x,  y')  dx 

J  a 

That  is,  the  integrand  E(x,  y')  does  not  depend  on  y.  Thus  Fy  =  and  Euler’s  equation 
reduces  to 


We  can  integrate  to  obtain  a  first  order  differential  equation 


Fy'=C 


where  C  is  a  constant  of  integration. 

Case  2:  Assume  the  functional  is  of  the  form 


F(y,  y)  dx 


That  is,  the  integrand  F(y,y')  does  not  depend  on  x.  In  this  case  it  can  be  shown  that 
Euler’s  equation  can  be  reduced  to 


A 

dx 


(F-y'Fyfi 


=  0 


2-4 


Similarly,  we  can  integrate  to  obtain  a  first  order  differential  equation 


F  -  y'Fy,  =  C 


where  C  is  a  constant  of  integration. 

Case  3:  Assume  the  functional  is  of  the  form 

F{x,  y)  dx 

where  the  integrand  F{x,y)  does  not  depend  on  y' .  Thus  Fyi  =  0  and  Euler’s  equation 
reduces  to 

Fy  =  0 

This  equation  is  not  a  differential  equation. 

Case  4-  Assume  the  functional  is  of  the  form 

F{x,y)\Jl  +  y'^  dx 

In  this  case  it  can  be  shown  that  Euler’s  equation  can  be  reduced  to 

v" 

Fy  -  F^y'  -  F  ^  .  =  0. 

^  1  +  y'^ 

This  case  often  occurs  when  the  functional  involves  integration  with  respect  to  an  arc 
length  s,  where 

ds  =  \Jl  +  yF  dx 

Often,  but  not  always,  the  “simplifications”  offered  by  these  special  cases  expedite 
the  process  of  finding  a  solution  to  the  Euler  equation. 

2.3. 1.2  Constrained  Optimization.  Some  problems  posed  in  the  Calculus 
of  Variations  involve  constraints,  or  subsidiary  conditions,  that  are  imposed  upon  the 
admissible  solutions.  The  isoperimetrie  problem  is  given  as: 
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Find  the  curve  y  =  y{x)  for  which  the  functional 

J[y]  =  /  F'{x,y,y')  dx 
J  a 

has  an  extremum,  where  the  admissible  curves  satisfy  the  boundary  conditions 

y{a)  =  A,  y(b)  =  B, 
and  are  such  that  another  functional 

^[y]=  [  G{x,y,y')dx 

J  a 

takes  a  fixed  value  1. 

Assuming  the  functionals  F  and  G  have  continuous  first  and  second  derivatives,  then 
by  Theorem  1,  Section  12.1  in  [16],  there  exists  a  constant  A  such  that  y  =  y{x)  is  an 
extremal  of  the  functional 

[F{x,y,y')  +  XG{x,y,y')]  dx, 
i.e.,  y  =  y{x)  satisfies  the  differential  equation 

The  constant  A  is  analogous  to  the  Lagrange  multiplier  used  in  parameter  optimization 
problems.  The  solution  of  equation  (2.4)  will  result  in  two  unknown  constants  of  integration 
and  the  unknown  constant  A.  These  three  unknowns  are  solved  by  enforcing  the  boundary 
conditions  y{a)  =  A,  y{b)  =  B  and  the  constraint  K[y]  =  1.  One  can  also  specify  an  initial 
guess  for  the  value  for  the  parameter  A  and  iterate  on  the  solution  until  the  constraint 
K[y]  =  I  is  satisfied.  The  latter  technique  is  useful  when  numerical  solutions  to  the  two 
point  boundary  value  problem  are  sought. 
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2. 3. 1.3  Weak  and  Strong  Extremum.  Not  only  are  we  concerned  about 
the  existence  of  an  extremum  for  our  functional  J[y\,  but  we  also  wish  to  classify  our 
extremum.  Thus  we  obtain  the  following  definitions: 

Definition  2.3.1.  For  y  =  y,  a  weak  extremum  of  the  funetional  J[y\  exists  if  there  is 
an  e  >  0  sueh  that  J[y]  —  J[y]  has  the  same  sign  for  all  y  whieh  satisfy  \\y  —  y||i  <  s  where 

||y||;^  =  max  \y{x)\  +  max  \y{x)\ 

a<x<b  a<x<b 

denotes  the  norm  in  the  spaee  of  all  eontinuous  funetions  that  have  eontinuous  first  deriva¬ 
tives  on  some  elosed  interval  [a,  b]  . 

Definition  2.3.2.  For  y  =  y,  a  strong  extremum  of  the  funetional  J[y]  exists  if  there 
is  an  £  >  sueh  that  J[y]  —  .J[y]  has  the  same  sign  for  all  y  whieh  satisfy  \\y  —  y\\Q  <  e 
where 

||y||o=  rnax  |y(x)| 

a<x<b 

denotes  the  norm  in  the  spaee  of  all  eontinuous  funetions  on  some  elosed  interval  [a,  b] . 

Thus,  every  strong  extremum  is  also  simultaneously  a  weak  extremum  and  any  nec¬ 
essary  condition  generated  for  a  weak  extremum  would  also  be  a  necessary  condition  for  a 
strong  extremum. 

2.3. 1.4  Neeessary  and  Suffieient  Conditions  for  a  Weak  Extremum.  As 
shown  in  Theorem  2.3.1,  the  vanishing  of  the  (first)  variation  5J[hi\  of  the  functional  J[y\ 
was  a  necessary  condition  for  the  existence  of  a  (weak)  extremal.  Analogous  to  examining 
the  second  derivative  in  parameter  optimization  problems,  to  show  sufficient  conditions  for 
a  weak  extremal  we  consider  the  second  variation 

By  applying  Taylor’s  theorem  to  functionals  of  the  form  (2.2),  with  an  increment 
h{x)  satisfying  the  boundary  conditions 

h{a)  =  0,  h{b)  =  0 
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one  can  obtain  the  second  variation 


^  J  (^Fyyh^  +2Fyylhh'  +  Fylylh''^^  (IX  . 


(2.5) 


Applying  the  condition  that  the  increment  must  vanish  on  the  boundary  and  integrating 
by  parts  yields  a  convenient  form  of  the  second  variation 


52 j[h]  =  ^  (Ph'^  +  dx, 


(2.6) 


where 


P  =  P{x)  =  -Fylyl,  Q  =  Q{x)  =  -  ^Fyy  —  ^ ^ 

Theorem  2.3.3  (Legendre).  A  necessary  condition  for  the  functional  of  the  form 

rb 

J[y\  =  /  Fix,y,y')  dx,  y{a)  =  A,  y{b)  =  B 
J  a 

to  have  a  minimum  for  the  curve  y  =  y{x)  is  that  the  inequality 


Fylyl  >  0 

(Legendre’s  condition)  he  satisfied  at  every  point  of  the  curve. 

Proof.  See  [16]  Sec  25.  □ 

Similarly,  Legendre’s  condition  for  a  maximum  is  Fyiyi  <  0. 

In  order  to  demonstrate  sufficient  conditions  for  a  weak  extremum,  we  must  introduce 
the  the  concept  of  conjugate  points  and  a  means  to  determine  conjugate  points  in  a  closed 
interval.  First  we  define  the  Jacobi  equation  of  the  functional  (2.2). 

Definition  2.3.3.  The  Euler  equation 

—^{Ph')+Qh  =  0  (2.7) 

of  the  quadratic  functional  (2.6)  is  called  the  Jacobi  equation  of  the  original  functional 

(2.2) 
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Next  we  define  the  conjugate  point  and  introduce  an  important  theorem  relating  the 
functional  from  which  the  Jacobi  equation  is  derived. 

Definition  2.3.4.  The  point  a  a)  is  said  to  be  eonjugate  to  the  point  a  if  the  Jaeobi 
equation  (2.7)  has  a  solution  whieh  vanishes  for  x  =  a  and  x  =  a  but  is  not  identieally 
zero. 

Theorem  2.3.4.  The  quadratie  funetional 

i^Ph'^  +  Qh"^^  dx,  (2.8) 

where 

P{x)  >0  {a  <  X  <b) 

is  positive  definite  for  all  h{x)  sueh  that  h{a)  =  h{b)  =  0  if  and  only  if  the  interval  [a,  b] 
eontains  no  points  eonjugate  to  a. 

Proof.  See  [16]  Sec.  26  Theorem  3.  □ 

Showing  that  the  interval  [a,  b]  contains  no  points  conjugate  to  a  can  be  accomplished 
by  solving  the  Jacobi  equation  (2.7),  thus  satisfying  the  definition  (2.3.4),  or  by  proving 
equation  (2.8)  is  positive  definite  and  invoking  Theorem  2.3.4.  Proving  equation  (2.8)  is 
positive  definite  can  be  accomplished  by  showing  Q  >0  over  the  interval  [a,  b].  While  this 
is  sufficient  to  show  the  functional  (2.8)  is  positive  definite,  it  is  not  necessary. 

Finally,  we  are  able  to  state  sufficient  conditions  for  finding  a  weak  extremum. 

Theorem  2.3.5.  Suppose  for  some  admissible  eurve  y  =  y{x),  the  funetional  (2.2)  satisfies 
the  following  eonditions: 

1.  The  eurve  y  =  y{x)  is  an  extremal,  i.e.,  satisfies  Euler’s  equation 

2.  Along  the  eurve  y  =  y{x), 

P{X)  =  ^Fy'y'  >  0 
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(the  strengthened  Legendre  eondition). 


3.  The  interval  [a,  b]  eontains  no  points  eonjugate  to  the  point  a. 

Then  the  funetional  (2.2)  has  a  weak  minimum  for  y  =  y{x). 

Proof.  See  [16]  Sec  28.  □ 

2.3. 1.5  Neeessary  and  Suffieient  Conditions  for  a  Strong  Extremum.  Every 
strong  extremum  is  also  a  weak  extremum,  however  the  converse  is  not  generally  true. 
Therefore,  all  of  the  necessary  conditions  for  a  weak  extremum  are  also  necessary  conditions 
for  a  strong  extremum.  In  order  to  state  necessary  conditions  for  a  strong  extremum  as 
derived  in  [16],  we  need  to  introduce  the  Weierstrass  E-function  (excess  function). 

Definition  2.3.5.  The  Weierstrass  E-funetion  of  the  funetional 

rb 

J[y\  =  /  F{x,y,y')  dx,  y{a)  =  A,  y(b)  =  B 
J  a 

is  given  by 

E{x,  y,  z,  w)  =  E(x,  y,  w)  -  F{x,  y,  z)  -  {w  -  z)Fyfx,  y,  z) 

Necessary  conditions  for  a  strong  extremum  are  given  by  evaluating  the  Weierstrass 
E-function  for  a  candidate  extremal.  The  conditions  given  below  are  for  strong  minima, 
however  they  are  extended  to  strong  maxima  by  reversing  the  inequality. 

Theorem  2.3.6  (Weierstrass’  necessary  condition).  If  the  funetional 

r-b 

J[y\=  j  F{x,y,y)dx,  y{a)  =  A,  y{b)  =  B  (2.9) 

J  a 

has  a  strong  minimum  for  the  extremal  7,  then 


E{x,y,y',w)  >  0 


(2.10) 


along  7  for  every  finite  w. 

Proof.  See  [16],  Section  34.  □ 
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2.3.2  Numerical  Methods.  There  are  many  methods  available  for  solving  varia¬ 
tional  problems  numerically.  The  two  methods  discussed  here  are  the  direct  method,  where 
the  cost  function  is  approximated  with  discrete  increments  and  the  shooting  method  which 
iteratively  seeks  numerical  solutions  of  the  two  point  boundary  value  problem  resulting 
from  the  Euler  equation. 


2.3.2. 1  Direct  Methods.  One  technique  for  minimizing  the  cost  functional 
directly  is  the  method  of  finite  differences.  Here  an  approximation  to  the  cost  funcional 


J  = 


F{x,y,y)  dx 


is  performed  by  using  numerical  methods  of  integration  and  differentiation.  This  implies 
that  the  trajectory  consists  of  some  finite  number  of  segments,  e.g.,  N.  For  example,  the 
derivative  terms  in  the  integrand  of  the  cost  functional  can  be  replaced  with  the  central 
difference  approximation 

.  ^  Ui+l  Ui—l 
2{ti+i-ti-i) 

where  yi  represents  the  derivative  of  y  at  the  segment  i.  Integration  can  be  performed 
using  the  trapezoidal  rule,  Simpson’s  rule,  quadrature  methods  or  any  other  numerical 
technique.  A  non-linear  programming  routine  such  as  Sequential  Quadratic  Programming, 
can  be  used  to  choose  the  N  points  that  minimize  the  approximate  cost  function.  As  N 
grows  large,  the  sequence  of  N  points  should  converge  to  the  optimal  trajectory.  Issues 
such  as  numeric  truncation,  roundoff  error  and/or  computational  burden  can  limit  the 
range  of  N . 

Another  direct  method  is  that  of  piecewise  linear  approximation.  Consider  approxi¬ 
mating  the  optimal  trajectory  from  points  A  to  H  with  a  series  of  N  straight  line  segments 
-  see,  e.g..  Fig.  2.1.  Writing  one  of  these  N  segments  in  the  two-point  form  of  a  line  in 
polar  coordinates,  we  have 


m 


rir2sin(02  -  9i) 
ri  sin  [9  —  9i)  —  r2  sin  {9  —  ^2) 


(2.11) 
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which  has  the  first  derivative 


m 


rir2  sin  (6i  —  62)  [ri  cos  {6  —  61)  —  r2  cos  {6  —  62)] 
[ri  sin  {6  —  61)  —  r2  sin  {9  —  62)]'^ 


(2.12) 


Now  the  cost  to  travel  from  the  point  (ri,  9i)  to  the  point  (r2,  62)  can  be  written  as 


Figure  2.1  Line  Segment  Defined  by  Points  (ri,0i),  (r2,02) 


Jl,2 


r62 


lei 


F{9,  R,  R)  dO 


(2.13) 


Assume  the  substitution  of  Eqs.  (2.11)  and  (2.12)  into  the  cost  function,  Eq.  (2.13)  yields 
a  closed  form^  for  Ji^2-  Eliminating  dependence  upon  6  allows  the  cost  of  any  given  line 
segment  to  be  explicitly  determined  for  a  given  pair  of  points  (ri,0i),  (r2,02).  The  total 
cost  of  the  optimal  path  for  some  number  N  line  segments  is  the  summation  of  the  costs 
for  each  segment,  viz., 

N 


r  =  }_^  Ji,i+i 

i=l 

The  original  variational  problem  can  now  be  approximated  using  discrete  straight  line 
segments  and  the  resulting  problem  solved  using  non-linear  programming  techniques. 


2. 3. 2. 2  Shooting  Methods.  The  shooting  method  is  a  numerical  technique 
of  solving  two  point  boundary  value  problems.  Here  we  are  concerned  about  solving  n 

^This  could  be  relaxed  by  substituting  numerical  integration  for  a  closed  form  expression  of  Ji,2. 
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first  order  non-linear  ordinary  differential  equations  over  a  finite  interval  where  r  unknown 
boundary  conditions  exist  at  the  initial  point.  All  boundary  conditions  are  given  at  the 
final  point.  In  the  shooting  method,  the  missing  initial  conditions  are  guessed  and  the 
resulting  system  of  ordinary  differential  equations  is  propagated  forward  in  time.  At  the 
final  time,  the  endpoints  are  compared  to  the  known  final  boundary  conditions  and  if  the 
difference  is  not  within  some  tolerance,  a  correction  is  made  to  the  initial  guess  and  the 
process  of  solving  initial  value  problems  is  repeated. 

Shooting  methods  can  be  applied  to  equality  constrained  variational  problems  where 
a  Lagrange  multiplier  is  used,  e.g.,  the  isoperimetric  problem.  In  this  case,  in  addition 
to  solving  for  the  missing  initial  condition,  the  value  of  the  Lagrange  multiplier  must  also 
be  guessed  and  the  shooting  method  must  be  iterated  until  the  constraint  condition  is 
satisfied. 

Shooting  methods  are  sensitive  to  the  initial  guess  of  the  derivative  information. 
Thus  this  technique  requires  trial  and  error  and  can  be  difficult  to  automate  compared 
to  a  finite  difference  method,  for  example.  Some  modifications  to  the  shooting  codes 
can  improve  their  “robustness”.  Finite  difference  methods,  for  example,  can  be  used  as  a 
starting  point  for  the  initial  guess  needed  by  a  shooting  method.  Similarly,  continuation  is  a 
technique  where  the  shooting  problem  is  formulated  such  that  it  depends  on  a  parameter. 
The  shooting  problem  is  solved  for  one  value  of  the  parameter,  presumably  where  the 
solution  is  readily  obtained.  The  derivative  information  is  used  as  a  starting  point  for  the 
next  value  of  the  parameter,  and  so  on,  thus  the  problem  is  “continued”  over  the  range 
of  the  parameter.  Shooting  methods  can  also  be  employed  “backwards”  in  time,  where 
the  initial  boundary  conditions  are  known  and  the  final  boundary  values  are  guessed. 
Furthermore,  multiple  shooting  methods  that  employ  forward  and  backwards  shooting,  as 
well  as  unknown  free  parameters  at  both  ends  of  the  domain,  have  also  been  developed. 

For  a  more  detailed  description  of  the  shooting  method,  the  reader  can  consult  the 
text  by  Roberts  and  Shipman  [40],  the  texts  by  Bryson  [8,  7]  or  the  many  others  that  have 
been  written  on  this  topic.  FORTRAN  codes  that  implement  the  shooting  method  are 
available  in  [36,  11,  43]. 
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2-4  Summary 

In  this  chapter,  the  relevant  elements  of  radar  and  optimization  theory  were  pre¬ 
sented.  The  main  components  of  the  Calculus  of  Variations  were  discussed  and  the  nu¬ 
merical  treatment  of  variational  optimization  problems  was  addressed. 

In  the  next  chapter,  the  radar  exposure  minimization  problem  is  motivated  from 
the  radar  equation,  Eq.  (2.1)  and  is  formulated  in  the  Calculus  of  Variations.  A  rigorous 
derivation  and  analysis  of  the  solution  of  the  radar  exposure  minimization  problem  for  a 
single  air  vehicle  will  be  presented. 
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III.  Radar  Exposure  Minimization 

3. 1  Introduction 

Many  different  performance  objectives,  or  metrics,  can  be  conceived  when  considering 
the  problem  of  path  planning  for  threat  minimization.  This  work  does  not  consider  metrics 
such  as  probability  of  detection  or  probability  of  tracking.  While  such  metrics  may  have 
appeal  in  operational  contexts,  they  require  knowledge  of  specific  radar  systems  and  their 
targets  in  order  to  be  credible.  The  problem  formulation  that  follows  is  deterministic,  and 
is  based  upon  the  physics  of  radar  signal  propagation.  Thus,  the  cost  functional  developed 
here  applies  universally  to  all  monostatic  radars. 

In  this  chapter,  the  single  vehicle  radar  exposure  minimization  problem  is  addressed. 
The  problem  is  posed  in  the  Calculus  of  Variations  and  an  analytic  solution  is  derived. 
Expressions  for  the  optimal  cost  and  optimal  path  length  are  provided.  An  optimal  heading 
angle  control  law  is  developed  as  well  as  analytic  solutions  for  the  optimal  maximum  range 
attained  from  the  radar. 

3.2  Problem  Formulation 

Given  a  radar  located  at  the  origin  O  of  the  Euclidean  plane,  it  is  desired  to  find 
the  optimal  air  vehicle  trajectory  that  connects  two  pre-specified  points  A  and  B  in  the 
plane  such  that  the  received  Radio  Frequency  (RE)  energy  reflected  from  the  air  vehicle 
is  minimized;  see,  e.g..  Fig.  3.1.  According  to  the  radar  equation,  Eq.  (2.1),  the  ratio  of 
the  received  RE  power  to  the  transmitted  RE  power  reflected  from  the  target  is  inversely 
proportional  to  where  R  is  the  slant  range  from  the  target  to  the  monostatic  radar. 
The  cost  to  be  minimized  is  then 

L  iwA 

where  v  is  the  (constant)  speed  of  the  air  vehicle  and  I  is  the  path  length. 

Consider  the  trajectory  in  Fig.  3.1  to  be  given  in  polar  form,  as  R  =  R{9).  Assume 
R{0)  is  a  single  valued  function,  viz.,  we  do  not  allow  backtracking.  We  have,  v  =  ^,  i.e.. 
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dt  =  ^,  and  ds,  the  element  of  arc  length,  is  given  in  polar  coordinates  by 

Substituting  into  the  cost  equation  we  then  obtain  the  functional 


(3.1) 


(3.2) 

(3.3) 


j[Rm  = 


\/i?2  + 


dO 


The  boundary  conditions  are 


ii(0)  =  Ro 

R{ef)=Rf,  0<9<ef 


Figure  3.1  Illustration  of  the  Radar  Exposure  Minimization  Problem 


3.3  Unconstrained  Analytic  Solution 

Without  loss  of  generality,  assume  Rf  >  Ro  and  0  <  9 f  <  tt,  see,  e.g..  Fig.  3.1.  Polar 
coordinates  are  used.  We  have  the  following: 

Theorem  3.3.1.  The  optimal  trajectory,  connecting  points  A  and  B  at  a  distance  Ro 
and  Rf  from  the  radar  located  at  the  origin  O,  where  9f  is  the  angle  lAOB,  minimizing 
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exposure  to  the  radar  aeeording  to  Eqs.  (3.1)-(3.3),  is 


R*{e)  =  R, 


3/  sin(36*  +  cj)) 


sin  ( 


,  0  <  0  <  0/ 


(3.4) 


where  the  angle 


( 


=  Arctan 


sin  30 


/ 


(3.5) 


\[ro)  -cos30/^ 

Moreover,  the  length  of  the  optimal  path  is  given  by  the  elliptie  integral  of  the  first  kind 


I*  = 


Rn 


^sin  (p  . 


_  2 

[sin(30  +  (p)]  3  dO 


(3.6) 


and  the  eost  funetion  explieitly  evaluates  to 


1  sin  30/ 
3Ro^  sin(30/  +  (p) 


(3.7) 


This  result  holds  provided  0  <  0/  <  |.  However,  */  f  <  0/  <  vr,  then  an  optimal  path  does 
not  exist  and  a  eonstraint  on  the  path  length  must  he  ineluded  to  render  the  optimization 
problem  well  posed. 


Proof.  We  have  obtained  a  variational  problem  with  an  integrand  which  is  not  explicitly 
dependent  upon  the  independent  variable  0.  In  this  case,  the  first  integral  form  of  the 
Euler  equation  [16]  can  be  employed: 


F  -  RF^  =  C 


where  F  is  the  integrand  of  the  cost  functional,  Eq.  (3.1),  i.e., 

^  Vr'^  + 


This  results  in  the  first  order  differential  equation 


1 

W 


cVr^  +  R^ 


(3.8) 
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where  C*  is  a  constant.  Thus, 


R  =  ±- 


'XjC’^  - 


where  ^  >  i?®  >  0.  Hence,  we  have  obtained  the  non-linear  ordinary  differential  equation 


The  integration  constant  C  will  be  determined  by  the  terminal  condition,  R{9 f)  =  Rf  = 
\OB\. 

Momentarily  assume  R{0)  is  properly  unimodal  on  0  <  0  <  9f,  viz.,  3  0  G  (O,0j] 
such  that  R{9)  is  monotonically  increasing  (decreasing)  on  (0,0],  and  is  monotonically 
decreasing  (increasing)  on  [0,0/].  At  0  =  0,  R{9)  is  maximal  and  ^1^  =  0. 

Let  R{9)  be  monotonically  increasing  on  0  <  0  <  0  and  let  R{9)  be  monotonically 
decreasing  on  9  <  9  <  9f.  Consider  0  <  0  <  0  where  R{9)  is  monotonically  increasing,  and 

dR  _  v^l/C2  -i?6 


Thus, 


d9  =  ,  dR 


The  solution  of  this  ODE  entails  an  integration.  To  this  end,  define  the  new  variable 


Hence, 


u  =  CR^,  i.e.,  du  =  3CR^  dR 


3  y/l  —  V? 


Integration  yields  u  =  sin(30  -|-  (/>),  where  (p  is  the  integration  constant.  Hence, 


R^{9)  =  —  sin(30  -f  (p) 


(3.10) 
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Therefore  on  0  <  0  <  0, 


R{e) 


1 

~W 


■^sin(36'  +  0), 


0  <  (/> 


Similarly,  on  9  <  9  <  9f, 


m 


1 

W 


^sm{39  -  Ip), 


0  <  '0 


(3.11) 


(3.12) 


where  C  >  0. 

We  have  three  unknowns:  </>,  V')  and  9,  and  three  conditions:  i?(0)  =  Ro,  R{9f)  =  Rf 

and  R{9)  =  max  R(9).  The  latter  yields  -  see,  e.g.,  Eqs.  (3.11)  and  (3.12): 
o<e<ef 

30 +  </>=!  (3.13) 

and 

3S-V.  =  -| 

Thus,  combining  Eqs.  (3.13)  and  (3.3)  yields 


(/>  +  V’  =  TT 


i.e.. 


^  =  TT  —  (/> 


Hence,  for  0  <  0  <  0j,  inserting  Eq.  (3.14)  into  (3.12)  yields 


R{9)  =  —  \/ sin(30  +  (p  —  tt) 

yC 

=  ^^sin(30  +  (/)) 


(3.14) 


Therefore,  the  formula 

^  ^sin(30  +  0) 

is  unique  and  applies  on  the  complete  domain  of  i?(0),  viz.,  it  applies  for  0  <  0  <  0j. 
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Subsequently,  we  use  the  boundary  conditions  i?(0)  =  Rq  and  R{9f  )  =  Rf  to  deter¬ 
mine  C  and  (j),  respectively,  viz., 

1 


Ro  =  R{0)  =  -ji=  ^sin( 

VC 


Solving  for  C  yields 


C  = 


sin  (j) 


Thus,  the  extremizing  trajectory  is  explicitly  given  by 


(3.15) 


R{e)  =  R, 


3/ sin(30  +  (j)) 


sin  ( 


(3.16) 


In  addition. 


which  yields 


Rf  =  RiOff  =  Ro 


( 


3  _  D  3  Sin(30/  +  f 


sin  (j) 


sin  30 


/ 


(j)  =  Arctan 

The  extremal  R{9),  Eq.  (3.16),  satisfies  the  necessary  and  sufficient  conditions  for  a  weak 
local  minimum,  as  well  as  the  Weierstrass  necessary  conditions  for  a  strong  local  minimum; 
see  Appendix  A. 

Once  the  optimal  path  R[9)  =  R*{9)  has  been  explicitly  determined,  it  is  possible  to 
calculate  the  path  length  of  the  optimal  trajectory.  The  path  length  is  given  by 


rd 


1=1^  fR‘^{9)+R‘^{9)  d9 
Jo 


(3.17) 


Substituting  (3.8)  into  (3.17)  yields 


rO 


(3.18) 
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Using  Eqs.  (3.4)  and  (3.15)  yields 


I*  = 


Ro 


sin  (j) 


Ro 


/sin(30  +  (/>) 


-2 


sin  (j) 


(10 


Ro 


■^sin  Ip 


[sin  (30  +  <(>)]  ^  dO 


(3.19) 


In  Sec.  3.8,  we  explicitly  show  that  the  path  length  integral  (3.19)  evaluates  into  an  elliptic 
integral  of  the  first  kind. 

The  cost  function  Eq.  (3.1)  can  be  simplified  by  substituting  Eq.  (3.8)  to  obtain 


1 


de 


Substituting  for  the  previously  determined  integration  constant  (3.15),  and  optimal  tra¬ 
jectory  (3.4),  yields 


J* 


sm  ( 


Ro^ 

Jo 

sin  (j) 

f 

CO 

o 

CO 

J<f> 

sin 

3 

0 

sin  (p 

\  nnt 

de 


z0f+4> 


1 


cot  X 


sin^(a;) 

39f+<l> 


de 


3i?  3 


[cot  (j)  —  cot  (30 J  -|-  - 


sin  30 


/ 


3i?o^  sin(30/  -h  (j)) 


(3.20) 


The  optimal  trajectory  is  parameterized  by  two  non-dimensional  parameters,  viz., 
Rf/Ro  and  0/.  The  optimal  trajectory  represents  the  trade-off  between  minimizing  the 
time  of  exposure  (path  length)  and  the  RE  power  received  by  the  radar  over  time,  due 
to  the  air  vehicle’s  proximity  to  the  radar.  The  optimal  trajectory  no  longer  exists  as 
the  angle  0/  ^  a  critical  angle.  In  other  words,  the  path  length  becomes  infinite  at 
this  critical  angular  separation  of  the  segments  OA  and  OB.  Eurthermore,  as  the  path 
length  I  oo,  the  cost  approaches  a  finite  value  -  see,  e.g.,  Eig.  3.2.  Beyond  the  critical 
angle  0/  =  |,  there  does  not  exist  a  path  that  minimizes  our  cost  function.  That  is,  the 
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aforementioned  trade-off  breaks  down  and  it  is  advantageous  for  the  air  vehicle  to  travel 
away  from  the  radar.  Thus,  for  a  path  length  constraint  must  be  included  to 

render  the  optimization  problem  well  posed. 

J*  Cost  Path  Length 


Figure  3.2  The  Cost  for  the  Radar  Problem  is  Bounded  as  the  Path  Length  Grows 
Without  Bound 


□ 


3.3.1  Relationship  to  Rose  Funetions.  The  optimal  path  given  in  Eq.  (3.4)  can 
be  written  in  the  form 

R*{9)  =  3/^^  \/sin(36>  (p) 

S/sin  (j) 

which  is  very  similar  to  the  equation  of  the  three-leaved  rose  function  or  Rhodonea  [14], 
given  in  polar  coordinates  as 

r  =  a  cos  36 


The  rose  function  amplitude  parameter,  a,  is  the  maximum  distance  the  rose  function 
attains  from  the  origin.  The  rose  function  is  rotated  about  the  origin  by  adding  a  parameter 
(j),  such  that 

r  =  a  sin  (36  +  (p) 


Fig.  3.3  shows  the  three-leaved  rose  function  plotted  for  (p  =  {0,  ^}.  For  the  radar  problem, 
we  are  only  concerned  about  the  leaf  in  the  first  quadrant  of  the  rose  function  plot,  viz., 
when  0  <  0  <  60°. 


In  our  expression  for  R*(6),  the  rose  function  parameter  a  is  represented  by  the  term 
,  and  the  sine  term  is  taken  to  the  third  root.  Since  the  distance  from  the  origin 
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Figure  3.3  Three  Leaved  Rose  Functions 


in  polar  coordinates  is  real  valued,  and  since  complex  values  of  R{0)  are  not  physically 
realizable  in  this  problem,  the  sign  of  the  original  rose  function  is  preserved.  In  other 
words,  we  are  only  concerned  about  the  first  “leaf’  of  the  rose  function  as,  evidently, 
0  <  9f  <  60°.  The  shape  of  the  rose  function  is  distorted  somewhat  by  the  cubed  root,  as 
shown  in  Fig.  3.4. 


3.3.2  Special  Cases.  Several  interesting  special  cases  concerning  the  optimal 
trajectory  given  by  Eq.  (3.4)  are  now  considered. 

In  the  case  where  Of  =  0  and  Rf  >  Rq,  the  origin  O  and  the  points  A  and  B  are 
colinear,  and  the  optimal  trajectory  is  a  straight  line,  as  shown  in  Fig.  3.5. 
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X 


Figure  3.5  Optimal  Trajectory  for  the  Special  Case  where  Of  =  0 


Also,  the  following  holds. 

Corollary  3.3.1.  The  optimal  trajectory  which  connects  points  A  and  B  at  a  distance 
Ro  =  Rf  from  the  radar  located  at  the  origin  O,  and  minimizes  exposure  to  the  radar 
according  to  Eqs.  (3.1)-(3.3),  is 


R*{e)  =  R, 


\ 


cos(30  — 


cos 


30 


0  <  0  <  0/ 


(3.21) 


where  Of  is  the  angle  lAOB.  The  length  of  the  optimal  trajectory  is  then  given  by 


I*  = 


Ro 


rdf  r 


3/  39f  Jn 

'  '  COS  — ^  ^ 


cos  30  — 


30  i 


dO 


(3.22) 


which  is  an  elliptic  integral  of  the  first  kind.  The  cost  function  explicitly  evaluates  to 


T*  2  .  f30f 

J  =  - o  sin  — k 

3Ro^  V  2 


(3.23) 


This  result  holds  provided  0  <  Of  <  ^. 


Proof.  When  Rf  =  Ro  we  can  write  Eq.  (3.5)  as 


(j)  =  Arctan 
=  Arctan 

=  Arctan 


/  sin30j  \ 
\  1  —  cos  30  f  ) 


vr  30  f 
2~  ^ 


(3.24) 
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The  optimal  trajectory  R*{6)  is  then  obtained  by  substituting  (3.24)  into  (3.4),  whereupon 
we  obtain 


R*{e)  =  Ro 


=  R 


o 


\ 

3 


sin(36<  +  f  —  ^) 
sin(|  -  ^) 

cos  (30  —  ^) 

Wf 

cos 


Similarly,  by  substituting  Eqs.  (3.15)  and  (3.24)  into  the  equation  for  the  path  length 
(3.18)  yields 


I*  = 


Ro^ 


sm 

Ro^ 


E  —  ^ 

2  2  I 


1 


COS 


rOf  ^ 

W  L  W 


de 


Substituting  the  expression  for  the  optimal  trajectory,  R*{9),  developed  for  this  special 
case,  we  obtain 


I* 


Ro^ 

Wf 


cos 


Wf 

)  Ro 


Ro  n 

^ cos  ^  -^0 


3 

Wf 

cos 

2 

\ 

cos  (30  —  ^) 

COS  1 

V  2  )\ 

2 

3 

2 

de 

de 


Finally,  the  cost  for  the  optimal  trajectory,  J*,  is  calculated  by  subsituting  Eq.  (3.24)  into 
Eq.  (3.20),  which  yields 


J* 


1 

1 

3^ 

2 

3R^ 


sin  39  f 

sin(30/  +  f  — 
sin  39  f 


cos 


2  ) 


(3.25) 

□ 
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In  the  symmetric  special  case  where  Rf  =  Rq,  it  is  interesting  to  note  that  the 
relationship  between  4>  and  Of  is  linear.  This  relationship,  given  by  Eq.  (3.24),  is  depicted 
in  Fig.  3.6.  Furthermore,  the  angle  cf)  is  evaluated  for  some  interesting  Of  in  Table  3.1.  We 
note  that  the  angle  0/  — >  ^  as  >  0. 


<^[deg] 


0f[deg] 


Figure  3.6  The  Constant  (f)  as  a  Function  of  0/  for  the  Special  Case  Rf  =  Rq 


Table  3.1  Interesting  Values  of  (p  for  the  Special  Case  Rf  =  Rq  =  1 


</> 

sin  (p 

0°+ 

90° 

1 

10° 

Arctan 

1 

2a/2-V3 

15° 

Arctan 

1 

V‘-* 

20° 

60° 

V3 

2 

O 

O 

CO 

45° 

I 

V2 

45° 

Arctan 

1 

60° 

0° 

0 

The  optimal  trajectory  (3.4)  is  shown  in  Fig.  3.7  for  the  case  where  Rf/Ro  =  1  and 
Of  =  45°.  By  inspection  of  Fig.  3.7  we  see  that  the  extremal  trajectory  is  indeed  symmetric 
when  Ro  =  Rf,  as  expected. 
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0.8 


Figure  3.7  An  Example  of  the  Optimal  Trajectory  R*(9)  for  the  Symmetric  Case 
Rf/Ro  =  l,df  =  45° 


3-4  Optimal  Heading  Angle 

In  this  section  we  are  concerned  with  calculating  the  optimal  heading  angle,  ip*  {6), 
of  an  aircraft  flying  the  radar  energy  minimization  trajectory  specified  by  Eq.  (3.4).  The 
heading  angle,  for  example,  could  be  useful  as  a  command  generator  in  an  autopilot  mech¬ 
anization. 


Figure  3.8  Angles  of  a  Small  Increment  of  6 
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Consider  Fig.  3.8,  where  we  travel  along  a  small  increment,  dO,  of  the  trajectory.  We 


can  write 


Rde 


sm  ip  = 


dS 


R 


\/R^  +  R? 

We  now  substitute  Eq.  (3.8)  into  Eq.  (3.26)  to  obtain 


(3.26) 


sin  v?  =  C  R^ 


(3.27) 


where  the  constant  C  is  given  by  Eq.  (3.15).  Substituting  Eq.  (3.10)  into  Eq.  (3.27)  yields 


sin  (p  =  sin(30  +  (j)) 


Thus,  we  have  obtained  a  useful  expression  for  the  optimal  heading  angle. 


p*{e)  =  ?,e  +  (t)  (3.28) 

Remark  3.4.1.  The  constant  (f>,  which  depends  only  upon  the  problem  parameters  ^  and 
6f  and  is  specified  in  Eq.  (3.5),  is  the  angle  of  departure,  pd- 

Proof.  The  angle  of  departure  is  given  under  the  condition  9  =  0.  Thus,  Eq.  (3.28)  becomes 

Pd  =  P*{0)  =  fi  (3.29) 

□ 

Likewise  the  angle  of  arrival,  pa,  is  given  under  the  condition  0  =  6 j,  thus 

Pa  =  p*{ef)  =  2,ef  +  (t>  (3.30) 


It  is  interesting  to  revisit  the  plot  of  (/>  vs.  0/,  given  in  Fig.  3.6,  for  the  symmetric 
special  case  where  Rf  =  Ro-  As  the  final  angle  Of  approaches  zero,  the  distance  between 
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Ro  and  Rf  becomes  small  and  the  optimal  trajectory  approaches  a  straight  line.  Thus,  the 
angle  of  departure,  or  (/>  as  depicted  in  Fig.  3.6,  approaches  90°.  Similarly,  as  6f  approaches 
60°,  the  path  length  of  the  extremal  approaches  infinity  and  thus  the  departure  angle  is 
shallow  and  approaches  zero. 


3.5  Solution  Triangle 


In  this  section  we  derive  an  alternate  form  of  the  extremal  path  formula  and  present 
a  graphical  interpretation  of  the  relationship  between  the  angles  </>,  0/  and  V’-  Assume  a 
properly  unimodal  trajectory,  viz.,  2)9 f+cj)  >  and  recall  the  optimal  trajectory,  Eq.  (3.4), 
written  here  as 


R*{e)  =  Ro 


I  sm{39  +  (f)) 


sm  ( 


(3.31) 


over  the  interval  0  <  0  <  0,  where  R*{9)  is  monotonically  increasing,  and  which  we  know 
extends  in  reality  to  the  domain  9  <  9  <  9f. 


Consider  the  interval  9  <  9  <  9f,  and  take  the  point  of  view  that  the  aircraft  flies 
from  point  B  toward  point  A  until  the  maximal  range  from  the  radar  is  reached.  The 
latter  is  larger  than  Rf  and  thus  Eq.  (3.31)  applies 


y  sin^; 


(3.32) 


for  some  >  0. 

Obviously,  for  9  =  9,  Eq.  (3.31)  yields 


R*{9)  =  R, 


=  Rn 


3/ sin(30  +  (j)) 


sm  ( 


■^sin  (j) 
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Similarly,  Eq.  (3.32)  yields 


R*{0)  =  Rf 


3  sm(3{9f  -6)+ Ip) 


sin  pj 


=  Rf 


\/smPp 


Therefore 

Ro^  ^  Rf^ 
sin  (p  sin  ■0 

Moreover, 

30  +  0=1 

and 

3(0y-0)  +  0  =  | 

Combining  Eqs.  (3.34)  and  (3.35)  yields 


0  +  0  +  30j  =  7r 


(3.33) 


(3.34) 


(3.35) 


(3.36) 


Hence,  Eqs.  (3.33)  and  (3.36)  imply,  in  accordance  with  the  “law  of  sines”,  that  the 
angles  0  and  0  are  the  angles  of  the  triangle  shown  in  Eig.  3.9. 


Figure  3.9  Solution  Triangle 
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The  “solution  triangle”  is  determined  by  its  sides  Rq^  and  Rf^  and  the  included 
angle  2>9f.  Consider  the  solution  triangle  in  Fig.  3.9.  We  see  that  indeed, 


sin  30  f 

tan  p  =  - g - 

(£)  -cos  30/ 

(3.37) 

sin  (p  = 

Ro^  sin  30/ 

(3.38) 

^JRo^  +  Rf^  -  2Ro^R/  cos  39  f 

cos  (p  = 

Rf^  —  Ro^  cos  30/ 

(3.39) 

^JRo^  +  Rf^  -  2Ro^R/  cos  39  f 

Similarly, 

sin  30  f 

tan  p  = - 3 - - - 

(^)  -cos  30/ 

Clearly,  in  the  symmetric  case. 

viz.,  Ro  =  Rf,  the  solution  triangle  is  isosceles  and  the 

angles  p  and  p  are  equal. 

Remark  3.5.1.  The  angle  p  is 

related  to  the  angle  of  arrival,  ipa-  Speeifieally, 

(Pa  =  TT 


Proof.  The  angle  of  arrival  is  given  in  Eq.  (3.30)  as  (pa  =  36 f  +  cj)-  Substituting  cpa  into 
Eq.  (3.36)  yields  the  result.  □ 

Using  Eq.  (3.36),  we  can  express  Eq.  (3.32)  as 


R*{e)  =  Rf, 


lsm{3{9f  -9)  p-ip) 


sin 


=  R, 


=  R. 


3  /  sin(7r  —  (f>  —  39) 
y  sin(7r  —  39 j  —  (f) 


3  /  sin(30  +  (j)) 
Y  sin(36'j  +  (f>) 


(3.40) 


Comparing  the  expression  for  the  optimal  trajectory  obtained  in  Eq.  (3.40)  with  the  form 
in  Eq.  (3.4),  we  notice  that  the  constant  terms  in  Eq.  (3.4)  involve  a  parameterization  in 


3-17 


the  initial  range,  Rq,  and  the  angle  of  departure,  (pd  =  4>-  Similarly,  the  constant  terms  in 
Eq.  (3.40)  involve  a  parameterization  in  the  final  range,  Rf,  and  the  angle  of  arrival,  pa- 

Remark  3.5.2.  Eq.  (3.40)  holds  for  the  complete  range  of  9,  viz.,  0  <  9  <  9f,  for  all 
admissible  Rq  and  Rf,  and  for  the  case  where  39 j  +  (/>  <  |.  Hence,  we  obtain  an  alternate 
expression  Eq.  (3.40)  for  the  optimal  trajectory. 


Proof.  Recall  the  optimal  trajectory,  Eq.  (3.4),  is  of  the  form 


R*{9)  =  ^^sm{39  +  j>) 


holds  for  the  complete  range  of  9.  Since  R{9f)  =  Rf,  we  have 


R’(e,)  =  Sf  =  ^  {/sin  (3  9/ +  0) 


Solving  for  C  yields 


C  = 


sin  {39f  +  (p) 


Rf^ 


Thus, 


R*{9)  =  Rf, 


I  sin(30  +  (p) 
sin(30/  +  (p) 


(3.41) 


□ 


Remark  3.5.3.  An  alternate  expression  for  the  optimal  cost,  J*,  given  in  Eq.  (3.20),  is 


r  = - — oJRo^  +  Rf^  -  2Ro^Rf^cos39f 

SRo^Rf^  V  °  1  °  1  1 


(3.42) 


Proof.  Recall  the  trigonometric  identity 


sin  {39 f  +  (p)  =  sin  39 /  cos  (p  +  cos  39 /  sin  (p 
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Substituting  Eqs.  (3.38)  and  (3.39)  yields 


sin  {36 f  +  (/))  = 


Rf^  sin  36  f 


RJ'  +  Rf'°  -  2Ro^Rf^  cos  36  f 


(3.43) 


Thus,  substituting  Eq.  (3.43)  into  the  optimal  cost,  Eq.  (3.20)  results  in 


J*  = 


3Ro^Rf^ 


Ro^  +  Rf^  —  2Ro^Rf^  cos  36  f 


We  can  also  express  the  optimal  cost  from  Eq.  (3.42)  as 


□ 


3Ro^Rf^J*  =  JrJ^  +  Rf^  -  2Ro^Rf^cos36f 


where  we  see  that  the  right  hand  side  is  now  the  distance  between  the  two  lines  having 
lengths  Ro^  and  Rf^  separated  by  the  angle  36f.  Thus,  one  can  visualize  the  relationship 
between  the  geometry  of  a  radar  problem  and  its  cost  -  see,  e.g.,  Eig  3.10. 


X 


Eigure  3.10  Graphical  Representation  of  the  Optimal  Cost  3Ro^Rf^  for  the  Case  Rq  = 
Rf  =  1.5  and  6f  =  30° 
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3.6  Alternate  Extremal  Trajeetories 

Recall  that  when  =  0,  the  optimal  trajectory  is  the  line  segment  AB.  The  cost 
of  this  line  segment  can  be  obtained  from  Eq.  (3.42)  by  letting  9f  =  0,  see,  e.g.,  Fig.  3.11. 
We  see  that  the  optimal  (minimal)  cost  for  a  straight  line  trajectory  is  given  by 


r 


3  \Ro^ 


(3.44) 


When  9f  ^  60°,  the  path  length  approaches  infinity  and  the  cost  is  given  by 


r 


1 

3 


(3.45) 


Letting  9f  =  tt/3,  we  find  the  cost  for  a  trajectory  of  infinite  length.  This  trajectory  is  the 
“go-around”  trajectory,  in  which  the  air  vehicle  attempts  to  minimize  it’s  exposure  to  the 
radar  by  beginning  at  the  point  A  and  running  straight  to  infinity,  then  returning  from 
infinity  in  a  straight  line  to  the  point  B,  see,  e.g..  Fig.  3.12.  We  have  the  following: 


Remark  3.6.1.  The  east  of  the“go  around”  trajeetory  has  the  same  value  as  the  eost 
funetion  for  tt  >  9f  >  ^,  i.e.,  Eq.  (3.45). 


Proof.  From  Eq.  (3.44),  evidently,  the  cost  of  progressing  along  a  straight  line  from  the 
point  to  oo  is 


Ja=J. 


1 

3^ 


and  similarly  from  oo  to  the  point  B 


TOO  _  jB  _ 


3  Rf 


See,  e.g..  Fig.  3.11  and  Fig.  3.12. 


□ 


Remark  3.6.2.  The  optimal  trajeetory,  R*{9),  in  Eq.  (3.4)  always  has  a  lower  eost  than 
the  “go  around”  trajeetory,  for  0  <  9f  <  60°. 
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Figure  3.11  Cost  for  a  Straight  Line  Trajectory,  0/  =  0 


Figure  3.12  Cost  for  a  Trajectory  Of  >  60° 


Proof.  In  Appendix  A,  the  extremal  trajectory  is  shown  to  be  a  local  minimum.  The 
“go-around”  trajectory,  has  a  cost  given  in  Eq.  (3.45).  For  all  Of  £  (0,7r/3), 

3  ^  SRo^Rf^ 

Thus,  the  cost  of  the  extremal  trajectory  is  less  than  the  “go-around”  trajectory.  □ 

Proposition  3.6.1.  The  optimal  trajectory,  R*{0),  in  Eq.  (3.4),  is  globally  optimal. 

Proof.  The  optimal  trajectory  R*(0)  is  a  unique  solution  of  the  Euler  equation.  There 
are  no  constraints  imposed  on  the  solution.  There  are  no  other  trajectories  satisfying 
the  boundary  conditions  having  a  lower  cost.  The  locally  optimal  solution  is  globally 
optimal.  □ 
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3. 7  Maximum  Range 

The  maximum  range,  attained  by  an  air  vehicle  on  an  optimal  radar  avoid¬ 

ance  trajectory  is  sought.  We  have  the  following: 


Remark  3.7.1.  Without  loss  of  generality,  let  Rf  >  Rq-  For  the  optimal  trajeetory  de- 
seribed  in  Theorem  3.3.1,  the  maximal  range  from  the  radar  is  explieitly  given  by 


R 


* 

max 


R*{e) 


3/  yjRf^+Ro^-2Rf^Ro^  cos  30/ 
sin  30/ 


for 


Rf  ,  for 


3^/  +  </>  >  f 

30 f  (p  —  ^ 


(3.46) 


Proof.  Consider  the  case  where  39f  +  (p  >  The  optimal  trajectory  is  then  properly 
unimodal.  Hence, 

~dB 

and  our  extremal  evaluated  at  the  point  9  =  9  \s  given  by 


R*{9) 


R 


O 


■^sin  (p 


From  the  solution  triangle.  Fig.  3.9,  we  see 

Ro^  sin  39  f 

sin  (p  =  —  - 

Y  -|-  i?o^  “  2Rf^Rj^  cos  39 f 

Thus 


^max(0)  =  R*iO)  =  Rc 


^R/  +  Ro^  -  2Rf^Ro^  cos  30 f 


Ro^  sin  30  i 


Rf^  -|-  Ro^  —  2Rf^Ro^  cos  39  f 


sin  30 


/ 


Now  consider  the  case  where  30/  In  this  case  the  trajectory  R*(0)  is 

entirely  monotonic  (increasing)  and  the  condition  ^1^  =  0  does  not  occur.  In  this  case, 
the  maximal  range  is  achieved  at  R{9 f)  =  Rf.  □ 
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Normalizing  the  expression  for  R*{9)  we  obtain 


R*{9)  " 

V 

'm 

1+1-2  1 

(fi) 

1  cos  30/ 

Ro  \ 

sin  30 

/ 

(3.48) 


Remark  3.7.2.  For  the  symmetric  special  case,  the  inequality  2>6 f  +  (f)  >  |  holds  for  all 
Ro  =  Rf  >  0  and  for  all  9 f  G  (0,  |).  Hence,  for  the  special  case  where  Rf  =  Ro,  the 
optimal  trajectory  described  in  Theorem  3.3.1  is  always  properly  unimodal. 


Proof.  Recall  Eq.  (3.24)  which  was  developed  for  the  symmetric  special  case  Ro  =  Rf, 


(p  = 


TT 

2 


% 

2 


Thus,  we  can  write 


39f  +  (j)  =  39f  +  ^  - 

39  f  IT 

~  ^  2 


% 

2 


TT 


for  any  admissible  9f 


□ 

Therefore,  for  the  symmetric  case  Ro  =  Rf,  we  can  write  Eq.  (3.48)  as 


,/x/2\/(l-C05  3l7T  „  ,  ^  ^  >r 


(3.49) 


Table  3.2  shows  Eq.  (3.49)  evaluated  at  several  interesting  values  of  9f  when  Ro  =  Rf.  A 
plot  of  the  maximal  distance  from  the  radar,  Rj^ax)  is  shown  as  a  function  of  0/  for  the 
cases  ^  =  1,2, 3, 4. 


3.8  Path  Length  Calculation 

Accurate  calculations  for  the  path  length  of  the  optimal  trajectory  are  desired.  In 
this  section,  we  show  that  the  optimal  path  length,  I*,  can  be  expressed  as  function  of 
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Table  3.2  for  Interesting  Values  of  0/  when  Ro  =  Rf  =  1 


Of 

R*m.JRo 

0°+ 

1 

10° 

^  (2  -  V3)  ® 

15° 

(4  -  2  v/2)  ® 

20° 

s/X 

V  V3 

O 

O 

CO 

1 

26 

45° 

(2  (2  +  ^2))® 

O 

O 

oo 

R*[6i] 


Figure  3.13  R'^^^/Rp  as  a  Function  of  6f  for  the  Cases  ^  =  1,2, 3, 4. 

elliptic  integrals.  While  not  closed  form,  this  formulation  leads  to  efficient  and  accurate 
calculation  of  path  length  without  resorting  to  numerical  integration. 

Consider  the  parametric  equations  for  an  elliptic  arc,  given  by 

X  =  a  cos  9 
y  =  b  sin  6 
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where  a  is  the  semiminor  axis  and  b  is  the  semimajor  axis.  The  length  of  an  elliptic  arc  is 
given  by 


I  =  /  _|_  y2 

Jo 


^ cos^  0  +  6^  sin^  0  dO 


=  b  \/l  —  kJ  sin^  9  d9 


(3.50) 


where  0  <  0  <  7  and  k  is  the  eccentricity  of  the  ellipse  given  by 


'==V'-F 


The  integral  given  in  Eq.  (3.50)  cannot  be  evaluated  in  closed  form  and  the  class  of  problems 
given  by  integrals  of  the  form  in  Eq.  (3.50)  are  called  elliptic  integrals  [1]. 

Specifically,  an  elliptic  integral  of  the  first  kind  is  defined  as 


F{i,k)=  r 

Jo 


—  k"^  sin^  9 


,  0  <k  <1, 


(3.51) 


and  when  7  =  7r/2,  Eq.  (3.51)  is  considered  a  complete  elliptic  integral  of  the  first  kind, 


denoted 


f^/2  ng 

K{k)=  /  ,  0<fe<l. 

Jo  vl  —  sin^  9 


(3.52) 


Values  of  the  elliptic  integral  functions  are  tabulated  in  mathematics  references  [10,  29]. 
Elliptic  integrals  can  also  be  evaluated  as  function  calls  in  popular  software  packages  such 
as  Mathematica  [45]  or  the  EORTRAN  IMSL  libraries  [43] . 

We  claim  that  the  path  length  integral  given  in  Eq.  (3.6)  as 


r  =  [  [sin(30  +  (/>)]  3  d9. 

Vsm  cj)  Jq 
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is  an  elliptic  integral.  This  is  shown  by  performing  a  change  of  variables.  The  first  trans¬ 
formation  9  =  —  (f)  yields 


1  Ro  1 

3  4  sini  X 


(3.53) 


For  the  second  transformation,  let  =  sinx.  Thus 


^7,2 

dx  =  - du 

cos  X 


(3.54) 


and 


srr 

u  =  vsm  X 


(3.55) 


Substituting  Eqs.  (3.54)  and  (3.55)  into  Eq.  (3.53)  yields 


I*  = 


Ro 


■  ^sin(36»7'+(/>) 


■^sin  (j) 


sec  X  du 


/sin  <f> 


We  use  the  trigonometric  identities 


to  obtain 


where 


a  =  sin  X  sec  x  = 
a  =  sin  X  ^  sec  x  = 


1 

-1 

TT^’ 


0  <  X  < 

TT 


TT 

2 


—  <  X  <  TT 
2 


r 


Ro  J  ^sin  4> 


■  ^sin(30/+</.) 


\/l  — 


Ai  = 


0,  0<36i/-h(/)<  f 
1  ,  I  <  36*/  -I-  (/)  <  TT 


(3.56) 


(3.57) 

(3.58) 


(3.59) 


3.8.1  Monotonically  Increasing  Trajectories.  When  0<36*/-|-(/><7r/2,  we  have 
the  case  of  a  monotonically  increasing  trajectory,  viz.,  example  is  shown 

in  Fig.  3.14,  where  we  plot  the  extremal  R*{9)  for  Rj/Rq  =  10,  and  Of  =  20°. 
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Figure  3.14  Monotonically  Increasing  Extremal  for  the  Case  Rj/Rq  =  10  and  Of  =  20° 

The  path  length  integral,  Eq.  (3.59),  for  a  monotonically  increasing  trajectory  can 
be  written  as 


il 

Ro 


1  (  l‘^sm{3ef+<P) 

^sincj)  \Jo 


(3.60) 


From  [10],  the  following  useful  relation  can  be  found  in  a  table  of  elliptic  integrals  of 
the  first  kind. 


rY 


dr 


1 


Vl-T^  2^ 


Fi^P,k) 


where 


2-^/3 

4 

Arccos 


i-(i  +  \/3)yA 
i  +  {V^-i)Yy 


Thus,  we  recognize  Eq.  (3.60)  as  the  difference  of  two  elliptic  integrals  of  the  first 
kind  which  has  a  solution  of  the  form 


1 


1 


I* 


Ro  2^ 


[F{'ipi,k)  -  F{ijj2,k)] 
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where  the  parameters 


i’2 


2-^/3 

4 

Arccos  ( 


Arccos 


1  -  (1  +  a/s)  [sin(30/  + 

1  +  (a/s  -  1)  [sin(30/  +  J 

l  +  {V3-l)[sm<J)fy 


(3.61) 

(3.62) 

(3.63) 


For  our  example  shown  in  Fig.  3.14,  where  Rq  =  1,  Rf  =  10,  and  Of  =  20°,  we 
compute  jf-  using 


sm  0  =  - 
2 


999001 


=  0.000866458 


■01  =  Arccos 


■02  =  Arccos 


l-5O055lim{l  +  ^) 


1  +  50 


23  - 


999001 

999001 


2s  +  ^ 999001  ^)  , 


F(0i,/c)  =  2.7229 
F(02,  fe)  =  0.250934 


2.66935  rad 


0.25076  rad 


Thus,  =  9.85108. 

3.8.2  Properly  Unimodal  Trajectories.  Now  consider  evaluating  the  path  length 
extremal  when  7r/2  <  30/  +  0  <  vr,  i.e.,  the  case  of  a  properly  unimodal  extremal.  In 
Sec.  3.7,  we  showed  that  30/  +  0  is  always  greater  than  7r/2  for  the  symmetric  case.  An 
example  of  a  properly  unimodal  trajectory  is  given  in  Fig.  3.7,  which  depicts  the  symmetric 
case  where  Rf/Ro  =  1 ,  and  Of  =  45° . 
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The  path  length  integral  expressed  in  Eq.  (3.59),  can  be  rewritten  as 


I*  _  1  (  f  du  \ 

Ro  ^s'mcp  yia/siitf  \/l  -  j 

v^sinc/)  y  Jo  \/l  -  n®  Jo  Vl-u^ 


+ 


du 

\/l  —  u® 

^sin(3ej-+(/)) 


du 

\/l  —  tt® 


(3.64) 


From  [10]  we  recognize  Eq.  (3.64)  as  the  difference  of  three  elliptic  integrals  of  the  first 
kind  and  has  a  solution  of  the  form 


I*  1  1 


Ro  2^  ^sincj) 


(3.65) 


where  k,  ■01  and  02  are  given  in  Eqs.  (3.61)-(3.63)  and 


03  =  Arccos 


'l-(l  +  ^/3)[sin(7r/2)]"/"^ 
l  +  (^/3-l)[sin(7r/2)]2/3 


=  TT 


(3.66) 


Using  some  fundamental  relations  of  elliptic  integrals  [10],  namely 


F{rmr  ±  7,  /c)  =  2mK{k)  ±  ^(7,  k) 
F{0,k)  =  0 


we  recognize 


F{'^3,k)  =  F{'K,k) 
=  2K{k) 
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Thus,  we  can  write  Eq.  (3.65)  as 


TT  =  -  ^(^2,  k)]  (3.67) 

Tto  2  V  3  V  sin  cp 

Remark  3.8.1.  For  all  symmetric  cases  when  Rf  =  Rq, 


sin  4>  =  sin  (30j  +  (/>) 


(3.68) 


Furthermore,  the  angle  of  departure  is  related  to  the  angle  of  arrival  by 


ipd  =  7r  -ipa 


Proof.  Using  Eq.  (3.24), 


=  sin  {36 f  +  (f) 


Similarly, 


TT  -  ipa  =  TT  -  {36 f  +  (/)) 

,  ^  TT  36  f. 

=  >r-(3«/+2  -— ) 

TT  36  f 

~  2~  ^ 

=  </> 

=  Pd 


□ 
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For  our  example  shown  in  Fig.  3.7,  where  Rf/Ro  =  1,  and  Of  =  45°,  we  show  how 
to  compute  the  path  length,  Using  Eq.  (3.68)  we  see  that  ipi  =  V’2  and 


sin  (j) 


i’l  =  '02 


=  0.382683 


\/4  +  2\/2 

'l-  (1  +  ^/3)  (4  +  2\/2) 


Arccos 


,1  +  (^/3-l)  (4  +  2^2) 


03  =  TT 


K{k)  =  1.59814 
F(0i,  k)  =  F{'ip2,  k)  =  1.93226 


1.89393  rad 


Thus,  £  =  1.3229. 

3.8.3  General  Result.  We  have  the  following: 

Remark  3.8.2.  The  optimal  path  length,  I*,  for  the  trajeetory  minimizing  the  exposure 
to  a  monostatie  radar  given  in  Theorem  3.3.1,  formulated  using  the  Legendre  form  of  the 
elliptie  integrals  is 


I* 


[AK{k)  -  F(0i,  k)  -  F(02,  A:)]  , 


Rp 

2?I 


[F(0i, /c)  -  F(02,  A:)]  , 


2>6 f  +  0  >  ^ 

30/  +  0<f 


where 


k 

01 

02 


2  -  v/3 


Arccos 


Arccos 


/ 1  -  (1  +  \/3)  [sin(36>/  + 

\1  +  (\/3  -  1)  [sin(36'/  +  0)]^^^/ 
A-(l  +  v/3)[sin0]"/"\ 
l^l  +  (V3-l)[sin0]2/£ 


□ 

Fig.  3.15  demonstrates  how  the  length  of  the  extremal  aircraft  trajectory  I*  varies 
with  increasing  Of;  the  path  lengths  of  optimal  trajectories  are  plotted  foi  Rf/Ro  =  1,  2,  3, 4 
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6»f[deg] 

Figure  3.15  I*  (j)f,  for  the  Cases  where  Rf/Ro  =  1,2, 3, 4 

and  Of  G  [0,59°].  In  Fig.  3.15,  I*  is  plotted  and  the  asymptotic  behavior  of  the 

path  length,  viz.,  I*  ^  oo  as  Of  ^  60°  is  evident. 

3.9  Summary 

In  this  chapter,  the  radar  exposure  minimization  problem  was  posed  in  the  Calculus 
of  Variations  and  a  closed  form,  globally  optimal,  solution  was  presented.  The  solution  was 
shown  to  exist  if  the  final  angle  included  between  the  endpoints  of  the  trajectory,  0  f  <  60°. 
The  solution  was  shown  to  be  related  to  the  class  of  3-leaved  rose  functions.  An  optimal 
heading  angle  control  law  was  determined  from  the  closed  form  solution.  Alternative  forms 
of  the  solution  were  presented  via  construction  of  a  geometric  tool,  i.e.,  a  solution  triangle. 
An  analytic  expression  for  the  optimal  cost  was  obtained  and  is  shown  to  approach  a  finite 
value  as  Of  ^  oo.  Finally,  the  optimal  path  length  was  derived  in  terms  of  elliptic  integral 
functions  and  the  asymptotic  behavior  of  the  path  length,  viz.,  I*  ^  oo  as  Of  ^  60°,  was 
demonstrated. 

In  the  next  chapter,  a  similar  analysis  will  be  performed  for  the  case  of  minimizing 
exposure  to  a  passive  sensor.  Closed  form  expressions  are  sought.  In  particular,  the 
relationship  between  the  solutions  of  the  passive  sensor  and  radar  exposure  minimization 
problems  is  revealed.  Lastly,  the  path  length  constrained  problem  is  posed  for  the  passive 
sensor  problem  and  the  solution  is  shown  to  be  a  function  of  an  elliptic  integral. 
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IV.  Minimizing  Exposure  to  a  Passive  Sensor 

4 . 1  Introduction 

The  problem  minimizing  the  energy  received  by  a  passive  sensor  or  emitter,  e.g.,  an 
IR  sensor  or  RF  jammer,  is  formulated.  In  the  following  sections,  a  closed  form  solution  to 
the  passive  sensor  problem  will  be  sought,  as  well  as  solutions  for  the  optimal  path  length 
and  optimal  cost.  A  comparison  will  be  made  between  the  structure  of  the  solution  of 
the  passive  sensor  problem  and  the  solution  of  the  radar  exposure  minimization  problem 
presented  in  Chapter  III.  Finally,  the  path  length  constrained  problem  is  addressed. 

4-2  Unconstrained  Analytic  Solution 

For  a  passive  sensor,  the  received  energy  is  inversely  proportional  to 
obtain  the  variational  problem 

\/R‘^  +  R'^ 

Jim]  =  I  de 

with  the  familiar  boundary  conditions 


o 

II 

o' 

(4.2) 

R{ef)  =  Rf,  o<e<ef 

(4.3) 

Without  loss  of  generality,  assume  Rf  >  Ro  and  0  <  Of  <  n.  Polar  coordinates  are  used. 
We  have  the  following: 

Theorem  4.2.1.  The  optimal  trajectory  which  connects  points  A  and  B  at  a  distance  Ro 
and  Rf  from  the  passive  RF  sensor  located  at  the  origin  O,  where  0 f  is  the  angle  lAOB, 
and  minimizes  the  exposure  to  the  RF  sensor  according  to  Eqs.  (4.1)-(4-3),  is 

R{6)  =  Ro  cos  9  +  — — 0<0f<'7r  (4.4) 

sm  0  f 

This  trajectory  is  the  arc  AB  of  the  circle  which  circumscribes  AOAB. 


R? .  Hence,  we 


(4.1) 
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Proof.  A  candidate  extremal  which  minimizes  the  cost  functional  (4.1),  can  be  found  by 
solving  the  differential  equation  that  results  from  the  Euler  equation 


(4.5) 

where  F  is  the  integrand  of  the  cost  functional,  Eq.  (4.1),  i.e.. 


and  we  have  the  derivatives 

^  -R{9f-2R{ef 

=  - , 

R{ef^R{ef  +  R{e) 

= m 

^  R{ef  yjR{ef  +  R{9f 
_  -3  R{9f  R{9f  -  2  R{9f  +  R{9f  R{9) 

R{9f(^R{9f  +  R{9fy 

Substituting  Eqs.  (4.6)  and  (4.7)  into  the  Euler  equation  (4.5)  yields 

-R{9)  -  R{9)  _  ^ 

(R{9f  +  R{9ff 

which  simplifies  to 

R{9)  +  R{9)  =  0 

Thus,  we  have  obtained  a  second  order,  homogeneous,  linear  differential  equation,  which 
has  the  solution 

i?(0)  =  Cl  cos  0  +  (72  sin  6*  (4.8) 


(4.6) 


(4.7) 
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Applying  the  boundary  conditions  (4.2)  and  (4.2)  we  obtain 


Ci  =  Ro 

Rf  —  Ro  cosOf 

C2  = - - 

sin  tif 

Thus  our  candidate  extremal  is 

R{e)  =  Ro  cos  e  +  ~  sin  0 

sin  Of 

Recall  the  equation  of  a  circle  passing  through  the  origin  in  polar  coordinates,  given  by 

r  =  acosO  +  bs'mO  (4-11) 

where  a  is  the  x-axis  intercept  and  b  is  the  y-axis  intercept.  We  recognize  Eq.  (4.4)  as  a 
circle  in  the  form  of  Eq.  (4.11).  The  circle  can  be  inscribed  by  a  triangle,  consisting  of 
three  points,  given  in  Cartesian  coordinates  as: 

(xi,yi)  =  (0,0)  ,  the  origin 

{x2,y2)  =  {Ro,  0)  ,  the  x-axis  intercept  and  the  point  A 

(xstVs)  =  {Rf  cos  6  f ,  Rf  sin  6  f)  ,  the  point  B 


(4.9) 

(4.10) 


Eigure  4.1  Triangle  Inscribing  the  Circle  Suggested  by  the  Optimal  Trajectory 
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The  circle  inscribing  the  triangle  AOAB  is  depicted  in  Fig.  4.1.  The  center  and 
radius  of  this  circle  can  be  identified  by  assigning  coefficients  to  the  quadratic  curve  [44] 


where 


+  dx  +  ey  +  f  =  0 


Xl 

yi 

X2 

02 

3^3 

2/3 

=  RoRf  sin  Of 


yi 

1 

d=- 

X2^  +  y2^ 

02 

1 

xs'^  +  y3^ 

03 

1 

=  —RfRo^sinOf 

xi^  + 

Xl 

1 

e  =  — 

X2^  +  y2^ 

X2 

1 

X3^  +  y3^ 

X3 

1 

=  RoRf  {Ro  cos  0 

f  - 

Rf) 

xi^  +  yi^ 

Xl 

01 

/  =  - 

X2‘^  +  y2^ 

X2 

02 

X3^  +  y3^ 

X3 

03 

=  0 


Remark  4.2.1.  For  the  optimal  trajectory  described  in  Theorem  4.2.1,  the  circle  circum¬ 
scribing  the  triangle  AOAB  is  centered  at 

f Ro  Rf  —  Ro cos 9 f\ 

VT’  2  sin  Of  J 
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and  has  radius 


ro  = 


\J Rf"^  +  Ro^  —  2RoRf  cos  Of 


2  sin  9  f 

Proof.  The  center  of  the  circle  can  be  identified  in  Cartesian  coordinates  as 


xo  = 

2/0  = 


d  Ro 

“  T 

e  Rf  —  Ro  cos9f 
2a  2  sin  9 f 


or  equivalently  in  polar  coordinates 


ro 


9o 


2RoRf  cos  9j 


Arctan 


2  sin  9f 

Rf  —  Ro  cos9f\ 
Ro  sin  9  f  J 


Thus,  the  radius  of  the  circle  is  given  by  tq.  □ 

The  candidate  extremal  is  defined  for  0  <  9  <  180°.  A  sample  trajectory  is  plotted 
in  Fig.  4.2  for  the  case  where  Rj/Ro  =  1  and  9f  =  60°.  Recall  that  for  the  radar  exposure 
minimization  problem,  the  path  length  of  the  optimal  trajectory  for  9f  =  60°  was  infinite 
and  for  9f  >  60°  an  optimal  solution  does  not  exist. 


Figure  4.2  Passive  Emitter  Case:  Extremal  R*{9)  for  Rj/Ro  =  1,  9f  =  60° 
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Corollary  4.2.1.  The  optimal  trajectory  which  connects  points  A  and  B  at  a  distance 
Ro  =  Rf  from  the  passive  RF  sensor  located  at  the  origin  O,  and  minimizes  the  exposure 
to  the  RF  sensor  according  to  Eqs.  (4-l)-(4-3)>  is 


(4.12) 


□ 
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The  center  of  the  circle  can  be  obtained  by  a  geometric  construction.  Consider 
Fig.  4.3,  where  the  angle  0/  is  bisected.  We  know  that  the  point  C,  halfway  along  the  arc 
AB  will  be  the  maximum  distance  from  the  origin.  Since  the  circle  passes  through  the 
origin,  the  distance  OC  must  be  equal  to  the  diameter,  or  twice  the  radius  of  the  circle. 
Furthermore  if  OC  is  the  diameter,  then  the  center  of  the  circle  exists  on  OC.  From 
Theorem  4.2.1,  we  know  the  center  of  the  circle  lies  on  the  line 


(4.13) 


Thus  the  intersection  of  the  bisector  OC  and  the  line  in  Eq.  (4.13)  is  the  center  and  is 
given  by 


Ro  Ro  .  Of 
— ,  —  tan  — 
2  2  2 


Remark  4.2.2.  For  the  optimal  trajectory  described  in  Theorem  the  optimal  path 

length,  viz.,  the  length  of  the  arc  AB,  is 


I 


^  Ro^  +  Rf^  —  2RoRf  cos  Of 


sin  0/ 


Proof.  The  path  length  is  given  by 


where 


I 


FR{e)  +  R‘^{e)  dO 


R{0)  =  —Cl  sin  0  +  C2  cos  6 
R^{0)  =  Ci^  sin^  0  —  2C1C2  sin  0  cos  0  +  (72^  cos^  0 
R?{0)  =  Ci^  cos^  0  +  2C1C2  sin  6  cos  0  +  €2^  sin^  9 


(4.14) 


(4.15) 

(4.16) 
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Adding  (4.15)  and  (4.16)  and  substituting  into  (3.17)  yields 


I  =  yjci^  +  c^^  de 

=  \lci^  +  C2^  de 
Jo 

=  0/  \ICi^  +  C2^  (4.17) 


Substituting  (4.9)  and  (4.10)  into  (4.17) 


Ro^  +  Rf^  —  2RoRf  cos  Of 
sin  Of 


we  obtain  the  path  length  of  the  extremal.  We  notice  that  the  path  length  is  now  expressed 
solely  as  a  function  of  the  boundary  conditions  Rf/Ro,  and  9f.  □ 


Remark  4.2.3.  For  the  optimal  trajectory  described  in  Corollary  4.2.1,  the  path  length, 
viz.,  the  length  of  the  arc  AB,  when  Ro  =  Rf,  is 


I 


sinf 


sin  Of 


Proof.  Let  Ro  =  Rj.  We  can  write  Eq.  (4.14)  as 


l  =  0f 


sin  Of 


smdf 


sinf 


=  2Roef—^ 
sm6f 


□ 
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Remark  4.2.4.  For  the  optimal  trajectory  described  in  Theorem  4-2.1,  the  cost  function 
explicitly  evaluates  to 

\J Ro^  +  Rf^  —  ^RoRf  cos  6*/ 


J  = 


RoR 


oXLf 


(4.18) 


Proof.  It  is  also  of  interest  to  obtain  a  closed  form  expression  for  the  cost  function,  given 
by 


J  = 


/o 


R2 


d9 


Substituting  Eqs.  (4.15)  and  (4.16)  we  obtain 


J  = 


0  (Cl  cos  6*  +  (72  sin  0)2 
Of 


=  \Ci^  +  CT 


de 

1 


/o  (Cl  cos  0  +  C2  sin0)2 

1^/ 


dO 


=  \Ci^  +  CT 


sin  0 


Ci^  cos  0  +  C1C2  sin  0 


=  \Ci^  +  C2 


sin  0 


/ 


Ciscos  0/  +  CiC2sin  0/ 


(4.19) 


Substituting  (4.9)  and  (4.10)  into  (4.19) 


Rf  —  Ro  cos  0/  \  ^  sin  0 f 

sin0/  J  R^^cos9f  +  Ro^i^^f^smef 

\J Ro^  +  Rf^  —  2RoRf  COS0/ 

RoRf 


□ 

We  have  obtained  an  expression  for  the  cost  which  is  only  a  function  of  the  boundary 
conditions  Rf/Ro  and  0/.  This  expression  can  be  useful  in  comparing  the  cost  of  the 
optimal  trajectory  to  sub-optimal  trajectories 

Remark  4.2.5.  For  the  optimal  trajectory  described  in  Corollary  4-2.1,  i.e.,  when  Ro  = 
Rf,  the  cost  function  explicitly  evaluates  to 

tan  ^ 

J  =  2 - (4.20) 

Ro 
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Proof.  Let  Ro  =  Rf-  We  can  write  Eq.  (4.18)  as 


J=  4 

=  2 

=  2 


2Ro^  —  2Ro^  cos  0/ 


Ro" 

1 

J— 

cos  9  f 

V 

2 

Ro 

tan  ^ 

Ro 

□ 


Using  the  expressions  obtained  for  the  path  length  I  in  Eq.  (4.14)  and  explicit  cost 
J  in  Eq.  (4.20),  we  see  that  as  the  final  angle  Of  tt/S,  the  path  length  I  oo.  The 
optimal  cost  remains  finite  and  approaches  J-;  see,  e.g.,  Eig.  4.4. 


Eigure  4.4  Asymptotic  Behavior  of  the  Optimal  Cost  and  Path  Length  for  Ro  =  Rf 


4-3  Passive  Sensor  Maximum  Range 

Remark  4.3.1.  For  the  optimal  trajeetory  deseribed  in  Theorem  4-3.1,  given  by 


R*{9)  =  Ro  cos  0  +  ^ sin  9 

sm  9f 

and  without  loss  of  generality  Rf  >  Ro,  the  maximal  range  from  the  passive  RF  sensor  is 
explieitly  given  by 


R 


* 

max 


R*{9) 


Ro^  +  —  2RoR f  cos  9 f 

sin  9  f 


(4.21) 
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Proof.  Assuming  the  extremal  is  properly  unimodal,  the  maximum  range  from  the  RF 
sensor  achieved  by  the  aircraft  is  obtained  when  =  0.  Thus, 


/ TTs  „  I  Rf  —  Ro  cos  0 f 

R{6)  =  0  =  —Ro  sin  0  +  (  — — ^ - -  |  cos  ( 


sm  ( 


7 


From  which  we  can  obtain 


tan  9  = 


Rf  —  Ro  cos  6 f 
Ro  sin  6  f 


Using  trigonometric  identities  we  obtain 


sin  9 


cos  9 


Rf  —  Ro  cos  9 f 
Ro^  +  Rf‘^  —  2RoR  f  cos  9  f 
Ro  sin  9  f 

^  Ro"^  +  Rf"^  —  2RoR  f  cos  9  / 


Substituting  (4.22)  and  (4.23)  into  R{9)  yields 


(4.22) 

(4.23) 


max  R{9)  =  R{9) 
o<e<ef 


Ro^  +  Rf^  —  2RoRf  cos  9f 
sin  9f 


We  also  recognize  that  we  can  now  write  the  path  length,  I,  as  a  function  of  the  maximum 
range  to  the  sensor 

I  =  R{9)  9f 


Clearly,  if  the  extremal  is  monotonically  increasing,  the  maximum  range  is  Rf,  as  Rf  >  Ro 
is  given.  □ 


4-4  Alternate  Derivation  for  Passive  Sensor 

The  variational  problem  we  have  obtained  in  Eqs.  (4.1)-(4.3)  can  also  be  solved  using 
the  first  integral  of  the  Euler  equation  [16],  since  the  integrand  of  (4.1)  does  not  depend 
upon  the  variable  9.  The  first  integral  of  the  Euler  equation  is  given  by 

F-R'  Fr>  =C  (4.24) 
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where  C  is  a  constant.  Evaluating  Eq.  (4.24)  yields 


yjR{ef  +  R'{ef 

Thus,  R{9)  =  ^  —  R?{6)^  where  ^  >  0.  We  have  obtained  a  non-linear  ODE, 

^  =  ±^^-R\  ii(0)  =  Ro  (4.25) 

where  the  integration  constant  C  will  be  determined  by  the  terminal  condition  R{9f)  =  Rf. 

We  assume  Ri9)  is  unimodal  on  0  <  9  <  9f.  Hence,  3  0  <  9  <  9f  such  that 
R{9)  is  monotonically  increasing  (decreasing)  on  [0,0],  and  is  monotonically  decreasing 
(increasing)  on  [9,9f].  At  9  =  9,  R{9)  is  maximal  and  ^1^  =  0. 

Consider  0  <  0  <  0  where  R{9)  is  monotonically  increasing,  and 

dR  _  v^l/C2  -i?2 

Hb  ~  ^  i  ■ 


Thus, 


d9  = 


dR 


V^l/C2  -  i?2 

The  solution  of  this  ODE  entails  an  integration.  Eq.  (4.26)  is  of  the  form 


(4.26) 


dx  =  Arcsin  — 
a 


—  x^ 


Let  a  =  1/C.  Thus  the  integration  of  (4.26)  yields 


9  +  (f)  =  Arcsin 


R 


where  (/>  is  an  integration  constant.  Hence, 


R{9)  =  jaj  sin(6l  -|-  (f)) 
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Similarly  on  9  <  6  <  0f, 

R{9)  =  — |a|  sin(0  —  ip) 

We  have  three  unknowns:  ip^  and  0,  and  three  conditions:  R{0)  =  Ro,  R{df)  =  Rf  and 

R(9)  =  max  R(9).  The  latter  yields 

0<9<9f 

e  +  (j)=- 
^  2 


and 


thus 


e-ip 


TT 

2 


(j)  +  1p  =  TT 
Ip  =  TT  —  (p 


Hence,  for  0  <  6*  <  6f^ 


R{9)  =  — |a|  sin(6l  +  (p  —  tt) 
=  |a|  sin(6*  +  (p) 


Therefore,  the  formula 


R{9)  =  |a|  sin(6l  +  (p) 


applies  for  0  <  6*  <  0/.  To  determine  the  integration  constants,  we  substitute  the  boundary 
condition  (4.2)  to  obtain 

R{0)  =  Ro  =  \a\  sin  cp 


Thus, 


a 


Ro 


sin  (p 


(4.27) 
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Applying  Eqs.  (4.3)  and  (4.27)  we  obtain 


m)  =  Rf 


Rq 

sin  (f) 

Ro 


sin 
Ro  sin  9 


sin(6'/  +  (j)) 

(sin  9 f  cos  (j)  +  cos  9 j  sin  (p) 


f 


tan  (j) 


+  Ro  cos  9  f 


Solving  for  tan  (p 


Which  implies, 


tan  (p  = 


i?o  sin  9f 
Rf  —  Ro  cos  9f 


sin  (p 


cos  (p 


Ro  sin  9  f 

\J Ro^  +  Rf^  —  2RoRf  cos9f 
Rf  —  Ro  cos  9 f 
\J Ro^  +  Rf^  —  2RoRf  cos  6*/ 


(4.28) 


(4.29) 

(4.30) 


Thus,  we  can  give  an  alternate  expression  for  the  extremal  given  in  Eq.  (4.4)  as 


\l Ro^  +  —  2RoRf  cos9f 

R(9)  =  T - i  sin(0  +  cP) 


sm  t)f 


(4.31) 


where 


=  Arctan 


i?o  sin  9f 


(4.32) 


—  Ro  COS  9f  ^ 

We  can  show  that  the  trajectory  from  Eq.  (4.31)  is  explicitly  the  same  as  Eq.  (4.4).  Using 
a  trigonometric  identity,  Eq.  (4.31)  becomes 


^  Ro^  +  Rf^  —  2RoRf  cos  9f 

R{9)  =  - ^ — - (sin  9  cos  (p  +  cos  9  sin  cp) 


sin  9  f 

Substituting  (4.29)  and  (4.30)  we  obtain 


^  Rf  Ro  cos  9 f  ^  ^  ^ 

sm9f 


(4.33) 
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4-5  Comparison  of  Radar  and  Passive  Sensor  Cases 

The  optimal  trajectories  for  minimizing  the  energy  reflected  to  a  monostatic  radar 
and  passive  RF  sensor  share  a  very  similar  structure.  Consider  the  equation  for  the  optimal 
trajectory  for  the  monostatic  case,  given  by 


R*{e)  =  R, 


3/ sin(30  +  (j)) 


sm  ( 


where 


=  R, 


3  /  sin  36  cos  <p  +  cos  36  sin  < 


sm  ( 


„  n  sm36  ^ 

=  Ro  \  - -  +  cos  36 

y  tan  cp 


=  Arctan 


sin  30 


/ 


t/  I  —  cos  36 


\\Ro 


f 


Substituting 


tan  6  = 


R„^  sin  36 1 


Rf  —  Ro^  cos  36 f 


into  Eq.  (4.34)  we  obtain 


(4.34) 


R*{6)  =  \  — - -^o  cos  36 f  ^  (4.35) 

y  sin  30/ 

Comparing  the  structure  of  Eq.  (4.35)  and  (4.33),  we  see  that  both  are  essentially  the  sum 
of  a  sine  and  cosine  term.  The  constants  that  multiply  the  sine  and  cosine  terms  are  very 
similar  in  structure  as  well.  Table  4.1  compares  the  equations  for  trajectories,  path  length 
and  maximal  range  for  the  monostatic  radar  and  passive  sensor  cases. 


4-6  Constrained  Analytie  Solution 

Eor  the  passive  sensor  exposure  minimization  problem,  an  optimal  unconstrained 
solution  exists  for  0  <  0/  <  vr.  Although  this  is  a  great  advantage  over  the  radar  exposure 
minimization  problem,  it  is  still  useful  to  to  consider  a  path  length  constraint  for  the 
case  of  a  passive  sensor.  The  path  length  constrained  formulation  naturally  applies  when 
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Table  4.1  Comparison  of  Radar  and  Passive  Emitter  Equations 


Radar  (1/i?^ 


Passive  Emitter  (1/i?^ 


R*{e) 


sSr  ~ 

sin(30 +  (/.)-§ 


i? 


*  rll 


3/  ^Rf^+Ro^-2Rf'-^Ro^  COS  30  / 
sin  30/ 


^i-^g^sin0  +  RoCos0 

^i?o^+/?/^  — 2Ro^?/  cos  0/  ^ 
sin0/  f 

COS  0/ 

sin0/ 


dealing  with  limited  fuel  resources  or  when  attempting  to  coordinate  the  actions  of  multiple 
vehicles. 

Consider  the  problem  of  minimizing  exposure  to  a  passive  sensor  located  at  the  origin 
O,  starting  from  some  initial  location  (i?o,  0),  terminating  at  the  point  {Rf,9f),  where  the 
length  of  the  trajectory,  I,  is  given.  This  path  length  constraint,  I,  can  be  greater  than 
or  less  than  the  optimal  unconstrained  path  length,  I*,  determined  previously.  Clearly, 
the  path  length  constraint  is  bounded  below  by  the  minimum  possible  path  length,  viz. 
I  >  ||^i?||;  see,  e.g.,  Eig.  4.5. 


Eigure  4.5  Notional  Comparision  of  a  Unconstrained  Optimal  Path,  Constrained  Op¬ 
timal  Path  and  a  Minimum  Length  Path  for  the  Passive  Sensor  Exposure 
Minimization  Problem 


^For  properly  unimodal  case,  i.e.,  when  30/  +  cj>  >  l‘2 
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Assuming  a  constant  velocity  for  the  vehicles,  we  wish  to  minimize  the  cost  functional 


j[R{e)]  = 


1  A2 


de 


such  that  the  admissible  paths  satisfy  the  conditions 


i?(0)  =  Ro 

R{ef)  =  Rf,  o<e<ef 

L[R{e)]  =  [  ^  Vr^  +  A2  de  =  i 
Jo 


where  L[R{9)]  is  a  functional  representing  the  path  length  constraint.  The  variational 
problem  with  subsidiary  conditions  is  an  isoperimetric  problem  -  see,  e.g..  Section  2. 3. 1.2. 
We  can  express  the  augmented  cost  functional  as 


J  = 


(  VR^  +  A2 


+  xVr^  +  A2 1  de 


(4.36) 


where  A  is  a  parameter,  sometimes  called  a  Lagrange  multiplier  [16] . 

Since  our  functional  J  does  not  explicitly  depend  on  0,  we  can  use  the  Euler  equation 
of  the  form 

F  -  RF^  =  C 


where  C  is  a  constant  of  integration  and 

F  =  +  xVr^  +  A2^ 

p  _  AA  ^  R 

^  Vr^  +  A2  r^Vr^  +  r^ 

Thus,  the  Euler  equation  becomes 


1  +  XR^ 
\/iF  +  A2 


C 
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Solving  for  we  obtain 


Thus, 


where 


R^  =  -^  {l  +  \R^f  -R 


C2 


R  =  ±\l^{l  +  \R^f  -R‘^ 


=  ±^^{1+  2Ai?2  +  \^R^)  -  i?2 


d9 


C2  +  C'VC'2  -  4A  -  2A 
2A2 

C2  -  C'v/C2  -  4A  -  2A 
2A2 


For  the  moment,  consider  only  ^  >  0.  We  have 


de  = 


dR 


^V^(i?2_a2)(^2_52) 


Integrating  both  sides 


de  = 


di? 


^V(ii2-a2)(d?2_52) 

dR 

V^(i?2  -  a2)  (i?2  -  62) 

dR 

V^(i?2  _  a2)  (i?2  _  62) 
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Now  let  t  =  ^  dR  =  a  dt,  thus 


^(0  +  <^)  = 


a  dt 


\J —  a?)  —  6^) 

a  dt 

y{_„2)(l_j2){_62)(i_|J(2) 

f  dt 


bj  y/(l  -  t2)  (1  _  ^2^2) 


where  /c^  =  ^,  which  is  an  elliptic  integral  of  the  first  kind,  denoted  F((p,  k),  where 


(p  =  Arcsin  (  — 
a 


Thus, 

A{9  +  «  =  F(Arcsin(^),?) 
Taking  the  inverse  of  the  elliptic  integral  and  solving  for  R  yields 


R{0)  =  asn 


bXje  +  cP) 

c 


a 

b 


(4.37) 


where  sn(<^,  k)  is  the  Jacobi  elliptic  function  which  is  the  inverse  of  the  elliptic  integral 
F{ip,k).  In  this  formulation,  the  desired  path  length  I  is  not  explicitly  present.  Instead, 
the  parameter  A  is  chosen  in  an  iterative  fashion  such  that  the  constraint  L[R(9)]  =  I  is 
satisfied. 

The  author  is  aware  that  the  solution  of  the  constrained  passive  sensor  problem 
was  obtained  independently  by  Zabarankin,  et.  ah,  [46],  and  publication  of  their  result  is 
forthcoming. 


4-7  Summary 

In  this  chapter,  the  passive  sensor  exposure  minimization  problem  has  been  addressed 
and  a  closed  form  solution  has  been  obtained.  The  optimal  solution  was  shown  to  be  a  circle 
and  its  radius  and  center  were  determined.  Additionally,  analytic  solutions  for  the  optimal 
path  length  and  optimal  cost  were  provided.  A  comparison  between  the  solutions  of  the 
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passive  sensor  and  radar  exposure  minimization  problems  reveal  an  interesting  relationship 
in  the  structure  of  the  two  solutions.  Finally,  the  path  length  constrained  problem  is 
addressed  and  the  solution  is  shown  to  be  a  function  of  an  elliptic  integral. 

Unlike  the  passive  sensor  problem,  a  tractable  analytic  expression  for  the  constrained 
path  length  problem  appears  out  of  reach.  In  the  next  chapter,  numerical  methods  are 
applied  to  the  solution  of  both  the  unconstrained  and  constrained  radar  exposure  mini¬ 
mization  problem.  Two  methods  of  solution  are  demonstrated  and  specific  limitations  of 
the  methods  are  discussed. 
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V.  Numerical  Methods  of  Solution 

5. 1  Introduction 

In  Chapter  III  we  showed  how  the  Calculus  of  Variations  is  used  to  find  the  optimal 
trajectory  that  minimizes  exposure  to  a  radar.  A  two-point  boundary  value  problem  is 
arrived  at,  which  has  a  closed  form  solution  given  by  Eq.  (3.4). 

Although  numerical  techniques  are  readily  available  for  solving  boundary  value  prob¬ 
lems,  e.g.,  [40,  7,  11],  the  insights  gained  from  analyzing  the  closed  form  solution  are  invalu¬ 
able.  For  example,  consider  the  analyst  who  immediately  employs  a  numerical  technique, 
such  as  a  shooting  method,  and  wishes  to  evaluate  trajectories  for  >  60°.  Countless 
hours  could  be  saved  if  the  analyst  had  the  insight  that  the  two-point  boundary  value 
problem  does  not  have  a  solution  for  Of  beyond  60°.  Even  more  pathetic  is  the  analyst 
whose  numerical  experiments  lead  him  to  the  incorrect  conclusion  that  no  solution  then 
exists,  for  any  value  of  Conversely,  direct  numerical  optimization  methods  that  employ 
gridding  or  coarse  discretization  of  the  problem  space  can  result  in  the  incorrect  conclusion 
that  there  always  exists  an  optimal  solution,  for  any  value  of  Of ,  including  Of  >  60°.  While 
numerical  results  will  be  obtained  for  Of  >  60°,  the  ensuing  discretized  optimal  trajec¬ 
tory,  barring  numerical  precision  issues,  will  invariably  ride  the  boundary  of  the  search 
space;  the  latter  is  arbitrarily  imposed  by  the  analyst.  Hence,  the  analytical  result  that  an 
optimal  trajectory  does  not  exist  manages  to  manifest  itself  in  the  numerical  work. 

In  this  chapter,  we  present  the  results  of  using  numerical  methods  to  obtain  the 
solution  to  the  problem  of  minimizing  the  radar  exposure  cost  functional,  and  compare 
them  to  the  analytical  results  obtained  in  the  previous  sections.  This  comparison  will  aid 
in  the  validation  of  the  numerical  methods.  The  numerical  methods  will  then  be  applied 
to  the  synthesis  of  optimal  trajectories  for  Of  >  60°  where  a  path  length  constraint  is 
imposed. 

5.2  Finite  Difference  Approximation 

This  direct  method  employs  a  piecewise  linear  approximation  as  described  in  Sec.  2.3.2. 
Consider  approximating  the  optimal  trajectory  from  points  A  to  H  with  a  series  of  N 
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straight  line  segments.  Neglecting  the  departure  and  destination  points,  the  path  can  also 
be  defined  as  a  series  of  A^  —  1  waypoints,  equally  spaced  over  the  domain  (0,  Of)  -  see,  e.g.. 
Fig.  2.1.  One  of  these  N  segments  can  be  expressed  as  in  the  two-point  form  of  a  line  in 
polar  coordinates  as  in  Eq.  (2.11)  which  has  the  first  derivative  given  by  Eq.  (2.12). 

Eor  example,  consider  the  cost  to  travel  the  straight  line  segment  from  the  point 
(ri,0i)  to  the  point  (r2,6*2) 


Jl,2 


9, 


dO 


Substituting  Eqs.  (2.11)  and  (2.12)  into  the  cost  function,  Eq.  (5.1)  yields 


(5.1) 


Jl,2 


\J ri^  +  r2^  —  2rir2  cos  A6 
4  (rir2  sin  AO)^ 


[4rir2A0  (cos  AO 


sin  AO) 


—  (ri^  -|-  r2^)  {2A0  —  sin2A0)] 


where 


AO  =  01-  02 


Thus  we  have  eliminated  any  dependence  upon  0  and  the  cost  of  any  given  line  segment 
is  explicitly  determined  for  a  given  pair  of  points  (ri,0i),  {r2,02)-  The  total  cost  of  the 
optimal  path,  for  N  line  segments,  is  the  summation  of  the  costs  for  each  segment,  viz.. 


N 


r  = 


2  =  1 


2  +  1 


In  lieu  of  a  more  sophisticated  gridpoint  allocation  algorithm,  the  approximation  strategy 
employed  is  to  increase  the  number  of  segments,  N,  until  a  user  specified  termination 
criterion  is  met,  e.g.,  increase  N,  until 


N 

2=1 


7V-1 

^  ^  Jk,k+1 
k=l 


<  e 


(5.2) 
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Figure  5.1  Comparison  of  Unconstrained  Numerical  and  Analytic  Optimization  Results 
for  1  and  3  Waypoints,  Rf  /Ro  =  1,  0 f  =  45°. 

where  e  is  the  user  defined  tolerance  on  the  accuracy  of  J* .  A  similar  strategy  can  be  used 
to  approximate  the  path  length,  where 

N-l 

I*  =  ^  yj ri‘^  +  r2‘^  -  2rir2  cos  {Oi  -  ^j+i) 
i=l 

The  piecewise  linear  approximation  accurately  portrays  the  discrete  approximation 
to  the  cost  function.  This  permits  the  variational  problem  to  be  approximated  using  non¬ 
linear  programming  techniques^.  The  optimization  results  that  follow  are  performed  in 
Matlab  using  a  Sequential  Quadratic  Programming  (SQP)  routine  developed  at  AFIT 
[23],  as  well  as  the  fmincon.m  function  from  the  Matlab  Optimization  Toolbox  [13]. 
Several  optimization  scenarios  are  considered. 

5.2.1  Unconstrained  Path  Length.  Consider  the  optimization  problem  as  stated 
in  Eqs.  (3.1)-(3.3).  Since  Of  <  60°,  no  constraint  on  the  path  length  need  be  imposed.  The 
purpose  of  this  optimization  is  to  compare  the  analytic  and  numerical  optimization  results, 
viz.,  calibrate  the  numerical  results.  Fig.  5.1  shows  the  results  of  numerical  optimization 


^This  approach  also  facilitates  the  development  of  a  solution  using  the  finite  element  method. 
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e,  [deg] 

Figure  5.2  Comparison  of  Unconstrained  Numerical  and  Analytic  Optimization  Results 
for  17  Waypoints,  Rf/Ro  =  I,  Of  =  45°. 

using  the  piecewise  linear  approximation  developed  above.  A  radar  exposure  minimization 
problem  is  solved  for  the  cases  of  1  and  3  waypoints  and  the  analytic  solution  is  also 
presented  for  comparison  purposes.  Interestingly,  the  numerically  optimized  paths  fall 
mainly  inside  the  optimal  curve. 

Employing  this  strategy  and  Eq.  (5.2)  with  e  =  0.0001,  a  trajectory  of  17  waypoints 
was  determined  and  is  plotted  in  Eig  5.2.  The  same  approximation  strategy  is  used  to 
determine  the  path  length  of  the  function.  Table  5.1  compares  the  costs  and  path  lengths 
of  the  1,3  and  17  waypoint  solutions  with  the  optimal  result  for  the  case  of  Rf/Ro  =  1, 
Of  =  45°.  We  can  see  that  this  piecewise  linear  approximation  converges  from  above  to 
the  optimal  solution  as  the  number  of  waypoints  is  increased. 

Difficulties  can  present  themselves  with  the  unconstrained  problem  when  numeri¬ 
cal  optimization  is  employed.  Clearly,  obtaining  good  initial  guesses  for  the  value  of  the 
discrete  waypoints  will  factor  into  the  convergence  and  performance  of  the  numerical  op¬ 
timization  method.  As  Of  7r/3,  numerical  methods  have  difficulty  keeping  up  with  the 
explosive  nature  of  the  solution  -  see,  e.g.,  Novy  [30]. 
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Table  5.1  Comparison  of  Cost  and  Path  Length  for  Unconstrained  Radar  Problem 
Rf/Ro  =  1,  6f  =  45° 


Waypoints 

Cost 

Length 

1 

0.681329 

1.212458 

3 

0.631784 

1.292690 

17 

0.616710 

1.321193 

Optimal 

0.615920 

1.322899 

5.2.2  Constrained  Path  Length.  As  we  have  shown,  for  >  60°,  a  path  length 
constraint  must  be  included  to  make  the  optimization  problem  well  posed.  Clearly  the 
minimum  path  length  from  points  A  to  R  is  the  length  of  the  line  segment  AB.  The 
problem  can  be  stated  as 

/■i  \/fP  +  R2 

minimize  J[R{6)]  =  /  - -r - dO 

Jo  ^ 

subject  to:  R(0)  =  Rq  =  1, 

R{9f)  =  Rf  =  l, 

I  <  RoOf  =  ^ 

O<0<j 

4 

The  inequality  constraint  I  <  RoOf,  was  chosen  such  that  the  path  would  be  constrained 
about  a  circle  of  radius  Rq.  Moreover,  for  9f  <  60°,  if  the  constrained  path  length  I  >  I* 
from  Eq.  (3.6),  then  the  constraint  is  not  active. 

Fig.  5.3  shows  the  results  of  the  constrained  numerical  optimization  for  Rq  =  Rf, 
9f  =  45°.  The  discretization  was  performed  with  8  waypoints  evenly  spaced  over  the 
domain  of  6.  As  expected,  the  constrained  solution  is  a  segment  of  a  circle  with  radius 
Ro  =  1.  The  unconstrained  analytic  solution  is  also  shown  for  comparison  purposes. 
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Figure  5.3 


Comparison  of  Constrained  Numerical  Optimization  and  Unconstrained  An¬ 
alytic  Optimization  for  Rf  /Ro  =  1,  9f  =  45° 


Lastly,  we  consider  a  case  where  Of  >  60° 


minimize  J[R{9)] 


IStt 

18 


^4 


dO 


subject  to:  i?(0)  =  i?o  =  1, 


R{9f)  =  Rf  =  l, 


0<9< 


ISvr 


and  either:  max  R(9)  <  4 


or 

I  <  ^maxR{d)<4 


First,  we  use  our  numerical  optimization  method  to  solve  the  maximum  range  constrained 
problem.  That  is,  we  impose  only  max  R(9)  <  4.  Using  an  iterative  scheme,  the  number 
of  line  segments  is  increased  until  the  user  defined  tolerance  on  the  cost,  e  <  0.0001.  The 
length  of  the  maximum  range  constrained  trajectory,  lmaxR{e)<4  =  12.667841,  and  the  cost 
is  0.687642.  Next,  the  path  length  constrained  problem  is  solved  using  a  similar  iteration 
scheme,  imposing  the  constrained  path  length  /max_R(6»)<4  previously  calculated.  The  cost 
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for  the  latter  trajectory  is  calculated  as  0.684614.  Fig.  5.4  depicts  the  results  for  the  two 
different  types  of  constraints. 


Figure  5.4  Comparison  of  Constrained  Numerical  Optimizations  for  Rq  =  Rf,  Of  =  130° 


5.3  Shooting  Method 

The  Euler  equation  is  a  necessary  condition  for  a  functional  of  the  form 

Av]  =  [  F{x,y,y)  dx  (5.3) 

J  a 

to  be  extremized,  assuming  J[y]  is  defined  on  the  set  of  functions  with  continuous  first 
derivatives  on  [a,  b]  and  satisfy  the  boundary  conditions  y{a)  =  A,  y(b)  =  B.  Solving  the 
Euler  equation  for  the  function(s)  that  extremize(s)  the  functional  in  Eq.  (5.3)  results  in  a 
two  point  boundary  value  problem.  Often  these  problems  are  non-linear  and  closed  form 
solutions  are  rare. 

The  shooting  method  is  a  numerical  technique  used  to  solve  boundary  value  prob¬ 
lems.  Essentially  the  shooting  method  involves  searching  for  the  missing  initial  conditions 
and  once  they  are  found,  proceeds  as  if  the  problem  were  an  initial  value  problem.  The 
algorithm  is  as  follows: 

1.  Develop  guesses  for  the  missing  initial  conditions. 
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2.  Numerically  integrate  the  differential  equation  over  the  range  [a,b]  and  save  the 
values  of  the  function  at  the  end  points. 

3.  Determine  the  difference  between  the  computed  end  points  and  the  desired  end 
points. 

4.  Modify  the  guesses  for  the  missing  initial  conditions  based  upon  the  discrepancy. 

5.  Iterate  until  the  computed  and  desired  end  points  agree  to  some  user  defined  toler¬ 
ance. 

This  algorithm  describes  a  “forward”  shooting  method.  It  is  possible  to  integrate 
the  equations  backwards  from  guesses  of  the  missing  final  conditions.  Similarly,  “multi¬ 
shooting”  methods  exist  to  integrate  from  both  sides  simultaneously,  requiring  the  solutions 
to  meet  somewhere  in  the  middle. 

Obtaining  good  guesses  for  the  initial  condition  in  step  1  is  essential  for  obtaining 
valid  solutions  from  the  integrator  in  step  2.  Clearly,  a  poor  guess  for  the  initial  condi¬ 
tion  can  result  in  numerical  instability  and  failure  of  the  differential  equation  integration 
process.  A  robust  root  finding  method  is  required  to  cope  with  failures  in  the  integration 
process  as  well  as  the  sensitivity  of  the  problem  to  small  changes  in  initial  conditions  (i.e. 
stiffness).  Although  there  are  challenges  involved  with  the  shooting  method,  it  is  an  im¬ 
portant  and  practical  method  for  solving  non-linear  boundary  value  problems.  Additional 
information  on  shooting  methods  is  available  from  many  sources  including  [11,  40,  36,  7]. 


5.3.1  Unconstrained  Path  Length.  Once  again,  consider  the  problem  as  stated 
in  Eqs.  (3.1)-(3.3),  where  we  wish  to  minimize  the  exposure  to  a  radar.  The  first  integral 
form  of  the  Euler  equation,  used  in  the  proof  of  Theorem  3.3.1,  resulted  in  the  non¬ 
linear  ordinary  differential  equation  given  in  Eq.  (3.9).  Taking  into  account  the  results  of 
Section  3.7,  we  can  define  the  following 


m 


^l/C2-i?6(0) 

WW}  ’ 

0, 


30  +  (p<  ^ 
30  +  (/.  =  f 
36  +  (p  >  I 


(5.4) 
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where 


1 


>  i?®  >  0 


and  C  is  a  constant  of  integration  dependent  upon  the  parameters  Ro/Rf  and  Of.  The 
shooting  method  can  be  applied  to  Eq.  (5.4),  where  the  constant  of  integration  C  is  the  un¬ 
known  initial  parameter.  However,  there  are  several  numerical  difficulties  in  this  approach. 
For  example,  the  power  of  the  R^{6)  term  can  present  numerical  sensitivity  issues  in  re¬ 
lation  to  the  other  terms  in  the  expression.  Additionally,  in  the  neighborhood  of  i?  =  0, 
numerical  roundoff  and  truncation  errors  can  result  in  the  term  —  R^{9)  yielding 

a  complex  result.  This  issue,  as  well  as  the  sign  change  at  the  point  i?  =  0,  can  cause 
numerical  integration  schemes  to  fail  to  converge.  Lastly,  except  for  symmetric  problems, 
knowing  the  value  of  0  where  we  have  the  sign  change  is  not  possible  without  first  having 
the  analytic  result  for  0  in  hand. 


If  we  apply  the  Euler  equation 


A. 

dO' 


rather  than  the  first  integral  form,  we  obtain  the  following  second  order  non-linear  ordinary 
differential  equation. 

R{e)R{9) +2R{e) +  3R^{9)  =0  (5.5) 

which  can  be  expressed  as  the  following  system  of  first  order  equations 


Yl{9) 

nid) 

with  boundary  conditions 


(5.8) 

Y\{Sf)  =  Rf 

(5.9) 

=  y2io) 

=  3Y^{9)  + 


2Yi{9) 

Yi{9) 


(5.6) 

(5.7) 
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Figure  5.5  Top:  Comparison  of  the  Analytic  Result,  Shooting  Method  and  the  Direct 
Method  for  Rf/Ro  =  10  and  0f  =  10°;  Bottom:  Error  in  the  Two  Numerical 
Methods 


I 


Fig.  5.5  is  a  comparison  of  the  shooting  method  with  the  analytic  result  and  the 
direct  method  described  in  the  previous  section.  The  problem  parameters  Rf/Ro  =  10 
and  Of  =  10°.  The  difficulty  the  direct  method  has  in  matching  the  results  of  the  analytic 
solution  can  be  attributed,  in  part,  to  the  location  of  the  grid  points  in  the  mesh.  The 
direct  method  chosen  here  employed  an  initial  grid  of  20  evenly  spaced  waypoints,  where 
the  number  of  grid  points  was  adaptively  increased  until  the  cost  function  converged  to  an 
acceptable  tolerance.  A  non-evenly  spaced  grid  should  be  able  to  capture  the  shape  of  the 
optimal  curve  more  effectively.  In  comparison,  the  shooting  method  uses  an  adaptive  step 
size  Runga-Kutta  integrator. 

Another  difficulty  faced  by  the  direct  method  is  the  relative  flatness  of  the  cost 
function  in  the  vicinity  of  the  optimal  solution.  As  Of  ^  the  direct  method  experiences 
increased  difficulty  in  differentiating  between  the  cost  of  adjacent  solutions.  The  shooting 
method,  on  the  other  hand,  is  more  robust  in  this  regard.  The  optimal  cost,  given  by 
Eq.  (3.20),  the  cost  of  the  solution  by  the  shooting  method  and  the  cost  of  the  solution  by 
the  direct  method  are  as  follows: 


Analytic 

Direct 

Shooting 

Cost 

Error 

0.33304469 

0.33306025 

1.55589E-05 

0.33304497 

2.77249E-07 

It  is  often  useful  to  employ  a  technique  known  as  continuation  when  dealing  with 
problems  involving  stiffness  or  when  having  difficulty  obtaining  good  guesses  for  the  initial 
conditions.  Continuation  is  a  technique  where  solutions  are  first  obtained  in  an  area  of 
the  problem  space  where  solutions  are  easily  found.  Next,  the  parameters  that  causing 
the  stiffness  are  slowly  adjusted  in  the  direction  of  the  desired  solution.  The  results  of  the 
successful  attempts  at  the  shooting  method  are  then  used  as  the  initial  guesses  for  the  next 
attempt. 

For  complicated  problem  descriptions,  direct  methods  have  an  advantage  in  their 
flexibility  to  specify  a  variety  of  constraints.  Methods  based  upon  variational  principles 
are  limited  in  this  regard,  however,  subsidiary  conditions  can  be  readily  imposed,  i.e.,  path 
length  constraints. 
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5.3.2  Constrained  Path  Length. 


The  constrained  path  length  radar  exposure 


minimization  problem  is  given  by 


.  VR^  + 

mm  /  dO 

R(e)  Jo  R^ 

(5.10a) 

L[R{e)]  =  [  ^  Vr^  +  R‘^  d9  =  l 

Jo 

(5.10b) 

O 

II 

o' 

(5.10c) 

R{9f)  =  Rf,  0<9<9f 

(5.10d) 

where  L[R{9)]  is  a  functional  representing  the  isoperimetric  constraint. 

To  solve  the  constrained  path  length  problem,  we  will  employ  the  Lagrange  multiplier 
method  whereby  the  constrained  problem  is  formulated  as  an  equivalent  unconstrained 
problem.  We  can  express  the  augmented  cost  functional  as 

where  A  is  a  real  valued  constant,  i.e.,  the  Lagrange  multiplier.  The  resulting  Euler  equation 
yields  the  non-linear,  ordinary  differential  equation 

R(e)[R[e)  +  \R^{e)]  +  -  XR^{e)  +  2[R{e)f  -  2XR^{e)[R{e)f  =  o 


+  xVr^  + 


de 


(5.11) 


which  can  be  expressed  as  the  following  system  of  first  order  equations 


ym  =  y2{o) 

_  3Y,\e)  -  xrm  +  2Yi{e)  -  2Ay/(g)y|(g) 

^  Yi{e)  +  XY^^iO) 


with  boundary  conditions 


>1(0)  =  i?o 

yiiOf)  =  Rf 


(5.12) 

(5.13) 


(5.14) 

(5.15) 
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Since  we  don’t  know  what  value  of  A  achieves  a  given  path  length,  we  must  introduce 
an  additional  level  of  iteration  in  the  shooting  method  algorithm.  That  is,  we  must  guess 
a  value  of  A,  solve  the  boundary  value  problem  using  the  shooting  method,  compute  the 
resulting  path  length,  compare  the  result  to  our  desired  path  length,  adjust  our  guess  of  A 
and  repeat  until  the  desired  path  length  is  achieved. 

The  numerical  integrator  employed  by  the  shooting  method  provides  as  its  output 
a  finite  number  of  points  along  the  curve  R*{9).  In  order  to  evaluate  the  path  length, 
the  piecewise  linear  approximation  discussed  in  Sec.  5.2  is  made.  The  number  of  solu¬ 
tion  steps  the  numerical  integrator  takes  is  increased  until  the  path  length  approximation 
converges  within  some  user  defined  tolerance.  Fig.  5.6  depicts  a  comparison  of  an  opti¬ 
mal  constrained  trajectory  calculated  by  the  shooting  method  and  by  the  direct  method 
previously  described. 


4. 


3. 

0) 

O) 

c 

CO 

CC 

2. 


1. 


Figure  5.6  Comparison  of  the  Shooting  Method  and  the  Direct  Method  for  Ro  =  Rf, 
9f  =  130°  and  the  Path  Length  I  =  12.667841 

5.3.2. 1  Existence  Issues  for  the  Variational  Constrained  Path  Length  Problem. 
We  are  concerned  about  the  existence  of  solutions  for  the  constrained  radar  exposure 
minimization  problem.  Specifically,  for  a  given  set  of  initial  conditions  {Rf /Ro,9f),  can 
an  optimal  trajectory  be  determined  for  any  path  length  1?  Equivalently,  in  the  Lagrange 
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multiplier  formulation,  what  are  the  restrictions  on  the  choice  of  A  and  how  does  that  effect 
the  range  of  the  path  lengths  that  can  be  achieved? 

Recall  that  for  the  completely  unconstrained  case,  Eq.  (3.1),  there  exists  an  optimal 
solution  with  a  corresponding  optimal  path  length,  I*.  This  solution  exists  iff  9f  <  vr/S. 
For  tt/S  <  Of  <  TT,  there  does  not  exist  an  optimal  solution,  the  cost  is  bounded  from  below 
and  is  not  attained.  The  path  length  of  the  optimizing  trajectory  approaches  infinity. 

The  constrained  problem  is  solved  using  the  Lagrange  multiplier  method  previously 
described.  There  exist  restrictions  on  the  selection  of  A  which  will  impact  the  achievable 
path  length  1.  Similarly,  there  are  restrictions  of  the  range  of  achievable  path  length. 

Consider  the  path  length  constrained  problem  for  the  case  where  0  <  Of  <  60°. 
Without  loss  of  generality,  we  have  Rq  =  1  and  Rf  >  Rq-  Furthermore,  momentarily  ^ 
assume  Rf  =  1. 

Remark  5.3.1.  For  the  constrained  radar  exposure  minimization  problem,  an  optimal 
solution  does  not  exist  for  I  <  2sm(0j/2). 

Proof.  The  Euclidean  distance  between  the  points  A  and  B  is  given  by 

R"^  +  Rj,  —  2  Ro  Rf  cos  Of 

Thus,  for  Ro  =  Rf  =  I,  A  I  <  2sm(0j/2),  the  isoperimetric  constraint  cannot  be  satisfied 
and  an  optimal  solution  does  not  exist.  □ 

To  achieve  the  minimum  possible  path  length,  imin  =  2sin(0j/2),  allow  A  ^  oo.  The 
optimal  path  length  is  achieved  by  letting  A  =  0.  Thus,  for  A  G  [0,oo],  the  optimal  path 
length  I  G  [2  sin(0//2),  T].,  where  I*  is  given  by  Eq.  (3.6).  Fig.  5.7  shows  a  plot  of  optimal 
constrained  path  length  trajectories  with  A  =  0,0.1, 1,^  oo,  for  a  case  where  Rf/Ro  =  1 
and  Of  =  45°.  For  this  case,  the  optimal  path  length,  I*  =  1.3229,  and  the  minimum 
possible  path  length  is  0.7654. 

^This  restriction  is  imposed  for  mathematical  convenience  and  can  be  removed. 
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0.8 


Figure  5.7  Optimal  Trajectories  for  a  Class  of  Path  Length  Constrained  Problems,  A  = 
0,  0.1, 1,  — >  oo 

Fig.  5.8  depicts  the  curve  of  optimal  path  length  versus  the  Lagrange  multiplier  A, 
for  the  case  Rf/Ro  =  1  and  9f  =  45°,  for  A  >  0. 

Now  we  consider  achieving  a  path  length  greater  than  the  optimal  path  length,  viz., 
we  seek  I  >  I* (Rf /Ro,9f).  In  aviation  parlance,  this  is  the  “loitering”  case.  We  again  refer 
to  the  augmented  cost  functional,  Eq.  (5.11),  and  to  achieve  a  loitering  trajectory  we  must 
let  A  <  0. 

Remark  5.3.2.  Let  0  <  9f  <  60°  be  given.  For  A  <  0,  there  does  not  exist  a  global 
minimum  for  the  augmented  uneonstrained  radar  exposure  minimization  problem,  given  by 
the  eost  funetional,  Eq.  (5.11). 

Proof.  Allow  an  extremal  trajectory  to  run  out  to  oo.  We  know  the  exposure  is  bounded, 
yet  the  path  length  is  infinite.  This  implies  that  the  cost  goes  to  — oo,  and  thus  a  (global) 
minimum  cannot  exist  for  the  augmented  unconstrained  radar  exposure  minimization  prob¬ 
lem.  □ 

Remark  5.3.2  does  not  preclude  the  existence  of  a  local  minimum  in  the  augmented 
unconstrained  minimization  problem.  Moreover,  this  local  minimum  could  be  a  global 
minimum  in  the  constrained  minimization  problem  given  by  Eqs.  (5.10a)-(5.10d). 
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Figure  5.8  Sample  Path  Length  as  a  Function  of  A  for  2sin(0j/2)  <  I  <  I* 

Proposition  5.3.1.  For  the  functional  given  in  Eq.  (5.11),  for  0  <  Of  <  60°,  there  exists 
a  critical  value  of  the  Lagrange  multiplier,  Ac  <  0,  such  that  for  A  <  Ac,  no  extremizing 
trajectory  is  a  minimum. 


Proof.  According  to  the  Legendre  condition  for  a  weak  (and  hence  strong)  local  minimum, 
it  is  necessary  that  the  inequality  F f^[6 ,  R{9) ,  R{9)]  >  0  be  satisfied  at  every  point  along 
the  extremal  trajectory  R{9),  where  F[9,  R{9),  R{9)]  is  the  integrand  of  the  functional 
given  in  Eq.  (5.11) 


Ff,f,[9,R{9),m] 


l  +  XR^{9) 


R‘^{9)  R^{9)  +  R^{e) 


Fjirt  —  as  the  condition  1  +  \R'^{9)  >  0  holds,  or  equivalently 


^  -  R^{9)  ’ 


yo<9<9f 


Under  these  conditions,  Legendre’s  condition  holds.  Thus,  there  exists  some  Ac  <  0  such 
that  if  A  <  Ac,  a  minimizing  trajectory  cannot  exist.  For  Of  >  60°,  we  require  A  >  0  to 
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render  the  problem  well  posed-see,  e.g.,  Theorem  3.3.1.  Thus,  this  result  applies  only  for 
0  <  0/  <  60°.  □ 

Summarizing,  we  see  that  there  are  two  limitations  to  the  existence  of  minimizing 
trajectories;  see,  e.g..  Fig.  5.9:  First,  the  path  length  cannot  be  shorter  than  the  Euclidian 
distance  between  the  points  A  and  B.  Secondly,  for  Of  <  60°,  the  Lagrange  multiplier 
is  constrained  at  the  point  A  <  Ac  and  consequently  there  does  not  exist  a  minimizing 
trajectory  when  the  path  length  is  constrained  to  be  of  length  I  <  Ic- 

It  is  interesting  to  note  that  in  the  radar  exposure  minimization  problem  we  intro¬ 
duced  the  path  length  constraint  in  order  to  deal  with  the  lack  of  existence  of  solutions 
when  Of  >  60°.  This  solved  the  problem  for  Of  >  60°,  however  it  introduced  a  new  set 
of  problems  when  constrained  solutions  are  sought  for  Of  <  60°.  Similar  insights  into  the 
radar  exposure  maximization  problem  are  provided  in  Appendix  B. 


Figure  5.9  Upper  and  Lower  Limits  on  the  Path  Lengths  of  Some  Constrained  Trajec¬ 
tories 

Fig.  5.10  depicts  curves  of  optimal  cost  versus  path  length  for  the  situation  Rf  /Ro  =  1 
and  for  final  angle  Of  =  45°,  50°, . . . ,  70°.  Bounding  these  curves  on  the  left  is  a  plot  of 
the  cost  and  path  length  of  a  straight  line.  Also  plotted  is  the  curve  of  the  unconstrained 
optimal  cost  versus  path  length.  We  notice  that  for  small  final  angles  the  optimal  cost  and 
path  length  are  essentially  the  same  as  the  straight  line  cost  cost  and  path  length,  viz., 
there  is  very  little  curvature  in  R*{0)  for  shallow  final  angles. 
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Figure  5.10  Optimal  Trajectories  of  Path  Length  Constrained  Problems  for  Rq  =  Rf, 
Of  =  45°, 50°, 55°,  60°,  65°,  70° 

Fig.  5.11  is  a  zoomed  view  of  Fig.  5.10.  Generating  the  data  required  for  these  curves 
can  be  numerically  intensive.  The  shooting  method  using  continuation  is  employed.  The 
constrained  path  length  curves  for  9f  >  60°,  are  generated  by  allowing  A  to  be  very  large, 
such  that  the  optimal  trajectory  is  practically  a  straight  line  and  the  initial  conditions  are 
easy  to  guess.  The  value  of  A  is  decreased  using  continuation  until  A  — >  0,  and  the  path 
length  increases.  For  Of  <  60°,  we  employ  a  similar  strategy,  except  we  allow  A  to  reach 
0  and  continue  such  that  A  <  0.  We  see  the  convex  shape  of  the  constrained  path  length 
curves  in  Fig.  5.11,  and  the  lowest  cost  path  is  obtained  where  the  Lagrange  multiplier, 
A  =  0.  These  curves  are  bounded  on  the  left  by  the  Euclidean  distance  limit,  and  on  the 
right  by  the  existence  issues  for  A  <  0. 

5-4  Summary 

In  general,  it  is  rare  that  closed  form  solutions  to  optimal  control  problems  are 
obtained.  In  this  chapter,  two  numerical  methods  were  developed  to  address  two  point 
boundary  value  problems.  A  piecewise  linear  approximation  is  developed  which  converts 
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Figure  5.11  Zoomed  Image  of  Optimal  Trajectories  of  Path  Length  Constrained  Prob¬ 
lems  for  Ro  =  Rf,  9f  =  45°,  50°,  55°,  60°,  65°,  70° 

the  functional  optimization  problem  to  a  parameter  optimization  problem.  While  this 
method  seems  more  sensitive  to  numerical  issues  than  the  shooting  method,  it  readily  han¬ 
dles  constraints.  The  shooting  method  demonstrates  excellent  numerical  properties  for  the 
radar  exposure  minimization  problem  and  is  suitable  for  the  path  length  constrained  prob¬ 
lem.  The  shooting  method  will  be  the  primary  numerical  tool  employed  in  the  remainder 
of  this  dissertation. 

In  the  next  chapter,  several  applications  of  the  methods  developed  thus  far  will  be 
presented.  First,  the  cooperative  control  problem  of  isochronous  rendezvous  will  be  ad¬ 
dressed  using  the  path  length  constrained  formulation  of  the  radar  exposure  minimization 
problem.  Next,  the  problem  of  one  vehicle  minimizing  exposure  to  two  radars  is  addressed 
numerically  and  analytically.  Finally,  a  suboptimal  method  is  developed  which  extends 
the  two  radar  problem  to  n-radars. 
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VI.  Radar  Exposure  Minimization:  Application  and  Extension 

6. 1  Introduction 

In  this  chapter,  we  examine  the  application  and  extension  of  the  results  of  the  pre¬ 
vious  chapters  on  radar  exposure  minimization.  First,  we  examine  the  cooperative  con¬ 
trol  problem  of  multiple  vehicle  isochronous  rendezvous  in  the  context  of  radar  exposure 
minimization.  The  solution  of  a  representative  two-ship  cooperative  control  problem  is 
presented.  Next,  the  problem  of  radar  exposure  minimization  is  extended  to  the  case  of 
two  radars.  Numerical  and  analytical  results  are  sought.  Finally,  a  sub-optimal  algorithm 
is  developed  which  handles  the  case  of  radar  exposure  minimization  for  n-radars  and  it’s 
limitations  are  discussed. 

6.2  Multiple  Vehicle  Isochronous  Rendezvous 

In  order  to  extend  the  single  vehicle  radar  exposure  minimization  problem  to  multiple 
vehicles,  we  must  allow  for  optimal  paths  of  prespecified  length.  Although  a  closed  form 
solution  of  the  Euler  equation  for  the  path  length  constrained  radar  problem  is  out  of 
reach,  numerical  methods  developed  in  Chapter  V  are  available  to  study  this  problem. 
The  shooting  method  is  employed  for  solving  this  class  of  problems. 

Consider  the  problem  of  coordinating  the  isochronous  rendezvous  of  two  or  more  air 
vehicles  such  that  they  reach  the  target,  at  the  same  time,  while  minimizing  their  cumu¬ 
lative  exposure  to  a  radar.  A  hierarchical  control  structure  is  envisioned.  Two  decision 
levels  are  considered. 

At  the  lower  level,  for  each  air  vehicle  a  constrained  path  length  optimal  control 
problem  is  considered.  The  costs  for  various  path  lengths  are  separately  calculated  for 
each  vehicle,  from  a  minimum  length  path,  i.e.,  a  straight  line,  to  some  maximum  length 
path.  It  is  envisioned  that  tables  of  cost  vs.  time  (path  length)  are  then  provided  by  each 
vehicle  to  the  higher  level,  a  central  decision  maker,  who  determines  what  time  of  arrival 
minimizes  a  composite  cost  function;  the  latter  is  the  cumulative  exposure  to  radar  of 
the  (two,  or  more)  air  vehicles.  The  air  vehicles  are  then  each  provided  with  an  optimal 
joint  time  of  arrival  by  the  central  decision  maker,  and  each  air  vehicle  selects  the  optimal 
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(minimum  exposure)  path  to  perform  the  rendezvous,  given  the  joint  time  of  arrival  (path 
length) . 

For  example,  consider  a  ship  initially  at  (ri,0i)  =  (1,0°).  Consider  a  second  ship  at 
(?’2,02)  =  (1-5, 10°).  With  a  radar  located  at  the  origin,  the  ships  must  rendezvous  at  the 
point  {rf,6f)  =  (1,45°);  see,  e.g..  Fig  6.1.  For  this  problem,  we  assume  the  vehicles  move 
at  the  same  constant  speed,  thus  time  of  arrival  can  be  equated  to  path  length. 


Figure  6.1  Cooperative  Isochronous  Rendezvous 

Determining  a  series  of  optimal  paths  of  varying  length  can  be  accomplished  numer¬ 
ically,  in  this  case  by  employing  the  shooting  method,  with  Eqs.  (5.12)-(5.15). 

Remark  6.2.1.  The  problem  of  the  seeond  ship  ean  be  solved  by  transforming  to  the 
equivalent  problem,  Rq  =  1,R/  =  1.5cosl0°  and  Of  =  35°.  The  path  would  then  be 
traversed  baekwards  in  time  and  shifted  by  45°,  i.e., 

9^=45° -6j,  i  =  l,2,...,N,  j  =  N,N-l,...,l 

In  Fig.  6.2  we  see  a  curve  representing  optimal  path  length  versus  cost  for  both 
vehicles.  We  note  that  the  optimal  path  length  for  vehicle  1  is  1.3229.  The  optimal 
path  length  for  vehicle  2  is  1.1646.  Fig.  6.3  reflects  the  cumulative  exposure  cost  function 
Ji  *  +J2-  Each  vehicle  is  commanded  by  the  central  decision  maker  to  plan  an  optimal 
path  of  optimal  length  1.2556  and  the  resulting  trajectories  are  shown  in  Eig.  6.4.  The 
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optimal  path  length  I*  =  1.2556  is  determined  by  the  central  decision  maker  to  minimize 


the  cumulative  exposure  +  J^- 


Vehicle  1  (Ro:  1.000,  Rf:  1.000,  0^45°) 


Figure  6.2  Curves  of  Optimal  Cost  vs.  Path  Length  for  Vehicles  1  and  2 

Proposition  6.2.1.  The  existence  of  an  optimal  solution  to  the  isochronous  rendezvous 
problem  is  not  guaranteed.  However,  the  proposed  algorithm  constructively  settles  the  ex¬ 
istence  issue,  and  when  an  optimal  solution  exists,  its  synthesis  is  provided. 

Proof.  Clearly,  if  the  individual  cost  curves  for  vehicles  attempting  an  isochronous  ren¬ 
dezvous  do  not  share  a  common  path  length,  no  composite  cost  curve  can  be  constructed. 
This  occurs  when  the  maximum  (minimum)  achievable  path  length  for  one  vehicle  is  less 
(greater)  than  the  minimum  (maximum)  achievable  path  length  of  the  second  vehicle  -  see, 
e.g..  Proposition  5.3.1  and  Remark  5.3.1.  □ 


6.3  Exposure  Minimization:  Two  Radars 

In  Chapter  III,  optimal  trajectories  for  radar  exposure  minimization  were  derived. 
The  scenario  comprised  of  one  radar  and  one  air  vehicle  (target).  In  this  section,  optimal 
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Figure  6.3  Composite  Cost  Curve  Identifying  the  Optimal  Time  of  Arrival 


Figure  6.4  Trajectories  for  the  Coordinated  Isochronous  Rendezvous  Example 
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trajectories  for  minimizing  the  exposure  to  two  radars  are  considered.  This  is  an  additional 
stepping  stone  in  the  development  of  optimal  trajectories  for  multiple  threat  avoidance. 
The  Voronoi  diagram  [33],  constructed  from  the  positions  of  the  threat  radars,  has  been 
suggested  as  a  means  to  plan  trajectories  for  multiple  threat  radar  avoidance  [12].  In 
Fig.  6.5,  the  Voronoi  diagram  for  a  12  randomly  located  threats  is  shown.  Path  planning 
involves  devising  a  path  onto  one  of  the  Voronoi  segments,  performing  a  search  using  the 
Voronoi  edges  to  build  a  minimum  cost  polygonal  path,  and  lastly,  exiting  the  Voronoi 
segment  to  arrive  at  the  destination.  The  “minimum”  exposure  path  connecting  points 
A  and  B  is  shown  in  Fig.  6.5.  Understanding  the  consequences  of  non-globally  optimal 
solutions  for  complex  optimization  problems,  such  as  the  control  of  autonomous  air  vehicles, 
is  essential. 


Figure  6.5  Voronoi  Diagram  for  Representative  Multiple  Threat  Avoidance  Path  Plan¬ 
ning  Scenario 

In  our  problem,  a  Voronoi  edge  is  the  locus  of  all  points  of  equal  power  separating 
the  two  radars.  Specifically,  for  the  case  of  two  equal  power  radars,  the  Voronoi  edge  is  the 
perpendicular  bisector  of  the  line  connecting  the  two  radars.  When  the  two  radars  have 
unequal  power,  the  Voronoi  edge  is  an  Apollonius  circle  [32,  19].  In  this  section,  we  are 
concerned  with  the  optimality  of  the  Voronoi  path  for  the  two  radar  exposure  minimization 
problem.  The  question  is  addressed  whether  it  is  optimal  to  take  the  Voronoi  path  between, 
or  go  around,  the  two  radars. 
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Consider  an  aircraft  exposed  to  illumination  by  two  tracking  radars,  as  depicted  in 
Fig.  6.6.  Polar  coordinates  are  used.  Without  loss  of  generality,  the  first  radar  is  located 
at  the  origin,  O,  while  the  second  radar  is  located  at  the  point  (1,0°).  We  wish  to  travel 
from  point  A  at  {Ro,0o)  to  point  B  at  {Rf,6f)  while  minimizing  the  cumulative  exposure 
to  the  two  radars. 


Figure  6.6  Optimal  Avoidance  of  Two  Radars 
As  shown  in  Fig.  6.6,  the  distances  from  the  radars  to  the  aircraft  are  given  by 

n  =  Rie) 


and 

r2  =  ^R^{e)-2R{e)  cos  9  +  1 

Furthermore,  if  we  consider  that  the  two  radars  have  a  power  ratio  of  (a  :  1)  we  have  the 
weighted  sum 


G[e,R{e)]  =  —^  +  \ 

r+  r2 

a  1 

~  R^{d)  ^  (R2(0)  _2R(0)  cos  6  +  1)^ 


(6.3) 
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We  wish  to  find  R{0)  which  minimizes  a  cost  functional  of  the  form 


r^f 

j=  /  F[e,R{e),R{e)]d9 

JOo 

f 


G[e,R{9)]jR^{e)+R^{e)  d9 


(6.4) 


Jdo 

and  satisfies  the  boundary  conditions  Ri9o)  =  Ro-,  R{9f)  =  Rf. 

A  necessary  condition  for  the  functional  Eq.  (6.4)  to  be  extremized,  is  that  the  the 
Euler  equation 


A 


FR--^Ff,  =  ^ 


Fr  FF^jj  FFfj^  —  0 


RR 


RR 


(6.5) 


be  satisfied.  We  can  write 


R 


Fj.  =  Gf,VR^  +  R‘^  +  G-  _ 

VR^  +  R^ 


=  G 


R 


ViF  +  R? 
R  „ 


+  i?2 


Thus,  we  obtain  the  partial  derivatives 


^  R  ^  R^-R^  ^ 

~  i?2  +  i?2^R  +  “  ^ 


i?2  +  i?2 


R^ 

iF  +  R^ 


-2RR  ^  R  ^ 
Frr  —  ~  ^F  +  .  Fr 

i?2  +  R2\  +  R^ 


F  =  ^ 

IF  +  i?2 


Fe 


(6.6) 

(6.7) 

(6.8) 
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We  can  now  substitute  Eqs.  (6.6)-(6.8)  into  the  Euler  Equation,  Eq.  (6.5),  yielding 


— Fr 

i?2  +  i?2 


R 

+  i?2 


Fg- 


r(rr-2rA 

— - =  0 

(.R2  +  i?2^ 


(6.9) 


From  Eq.  (6.4)  we  have 


F  =  G\/R‘^  + 


(6.10) 


and  we  obtain 


Fg  =  GeVii^  +  R^ 


(6.11) 


Fr  =  GrVR^  +  +  G 

Vr^  +  R^ 


(6.12) 


Substituting  Eqs.  (6.10)-(6.12)  into  Eq.  (6.9),  and  since  R>0  and  R^  >  0,  we  obtain 


R^Gr  -  RGg 


{r  +  r)+2RG  =  {) 
i?2  +  i?2  V  ) 


From  Eq.  (6.3)  we  have 


Gr 


Gg 


4  (cos  6  —  R)  4a 

(R^-2Rcos0  +  lf 
4  R  sin  6 

{R^  -  2Rcose  +  lf 


(6.13) 


(6.14) 

(6.15) 


Eq.  (6.13)  can  be  expressed  as 


Now, 


R^ 


G  (  R 


Gg  Vi?2  +  i?2 


R  +  R]  -  2]  -R^^+R  =  0 
Gg 


G  _  R^  —  2R  cos 9  +  1 
Gg  4R  sin  9 


1  + 


a  {^R?  —  2R  cos  0  +  l) 

If 


21 


(6.16) 


(6.17) 
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and 


Gr 

Ge 


1 


R  sin  0 


cos  9  —  R  —  a 


{R‘^  -  2Rcos9  +  iy 


Substituting  Eqs.(6.17)  and  (6.18)  into  Eq.  (6.16),  we  obtain 


(6.18) 


1  R 

2  R2  +  R2 


(^r  +  r) 


2  R  cos  9  +  1  +  a 


(^R?  —  2  R  cos  6*  +  l) " 


+  R^  —  2  R  sin  9  —  1  +  a 


R^  J 

{R^  -  2Rcos9  +  iy 

^4 


0  (6.19) 


An  analytic  solution  to  Eq.  (6.19)  appears  out  of  reach.  Hence,  numerical  methods 
must  be  applied  to  examine  the  nature  of  the  extremal  trajectories.  Solving  Eq.  (6.19)  for 
R{9),  we  obtain 

2{R‘^  +  R^)  (  R^{l-R^  +  2Rsm9)-a{R^ -2Rcos9  +  lf\ 

R  =  —R  H - ;r -  - (6.20) 

R  \R*{R^ -2Rcos9  +  l)  +  a{R‘^ -2Rcos9  +  lf  J 

which  is  numerically  integrated  using  the  shooting  method. 

As  an  example,  let  each  radar  be  of  equal  power,  i.e.,  a  =  1.  Suppose  we  wish 
to  travel  from  point  A,  located  at  (—1/2,— 1/2)  to  point  B,  located  at  (1/2, 1/2),  while 
minimizing  the  cumulative  exposure  to  both  radars.  The  resulting  optimal  trajectory  is 
shown  in  Eig.  6.7;  we  see  that  while  a  swerve  maneuver  is  initially  performed,  the  solution 
of  the  Euler  equation  is  essentially  “attracted”  to  the  Voronoi  edge,  i.e.,  the  perpendicular 
bisector  of  the  line  joining  the  two  radars. 

Several  scenarios  were  useful  in  validating  the  algorithm  used  to  generate  trajectories 
for  the  two  radar  case.  Eirst,  both  radars  were  placed  at  the  origin  and  the  start  and 
endpoints  were  configured  such  that  9f  —  9o  <  60° .  The  results  of  the  numerical  integration 
were  compared  to  the  analytic  result  from  Eq.  (3.4).  Several  cases  were  generated  and  the 
trajectories  and  path  lengths  agreed  with  the  analytic  results  to  the  numerical  precision 
specified.  The  cost  calculated  by  the  numerical  scheme  was  twice  that  of  the  one  radar 
case.  Additionally,  several  cases  were  examined  where  the  line  segments  connecting  the  two 
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0.5 


Figure  6.7  Trajectory  for  Two  Radar  Exposure  Minimization,  Obtained  by  the  Shooting 
Method,  Radars  Located  at  (0,0)  and  (1,0) 

radars  and  the  points  A  and  B  lay  on  the  Voronoi  edge.  As  expected.  Fig.  6.8  shows  that 
the  extremizing  curve  obtained  by  the  shooting  method  is  a  straight  line.  The  optimality 
of  the  Voronoi  edge  is  discussed  in  the  next  section. 

From  the  analytical  results  obtained  in  the  one  radar  exposure  minimization  problem 
in  Chapter  III,  it  is  also  evident  that  for  the  case  of  multiple  radars,  a  path  length  constraint 
would  be  required  to  make  some  problems  well  posed.  For  example,  if  we  move  the  second 
radar  in  Fig.  6.7  from  (1,0)  to  (0,0),  we  recover  the  one  radar  exposure  minimization 
problem  with  the  angle  6f  included  between  the  radials  to  points  A  and  B,  180°  >  9f  >  60°. 
We  already  know,  for  this  case,  the  unconstrained  problem  is  not  well  posed.  Other 
geometries  could  also  be  ill  posed. 

For  the  constrained  formulation,  we  wish  to  find  the  optimal  trajectory  R*(9)  which 
minimizes  exposure  to  both  radars,  connecting  points  A  and  B  while  the  subsidiary  con¬ 
dition 

L[R{9)]  =  sj  R^  {9)  +  R^  {9)  =  I 
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Figure  6.8  Trajectory  for  Two  Radar  Exposure  Minimization,  Obtained  by  the  Shooting 
Method,  Radars  Located  at  (0,0)  and  (1,0) 


applies.  Thus,  we  wish  to  find  R*{6),  which  minimizes  the  augmented  functional 

= r  {im  ^  ^ ^ 

where  A  is  a  Lagrange  multiplier.  Thus, 


G[e,R{e)] 


a  1 

R^iO)  ^  {R2{d)-2R{e)  cos  6  +  1)"^ 


+  A 


(6.22) 


depends  only  on  9  and  R{9),  and  Eqs.  (6.16)  applies.  The  partial  derivatives  Gr  and  Gg 
remain  unchanged.  Hence  the  ratio  ^  given  in  Eq.  (6.18)  applies  and 


R^  —  2  R  cos  0-1-1 
4  R  sin  9 


1  +  (^  +  a)  (r2_2Rcos0  +  1)^ 


(6.23) 
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Substituting  Eqs.(6.23)  and  (6.18)  into  Eq.  (6.16),  we  obtain 


\^^^{r  +  r)  [r^-2Rcos9  +  1+[^  +  x)  {R^-2Rcose  +  lf) 

+  R^  -2Rsme  -1+  (R^  -2Rcos9  +  l)^  =  0  (6.24) 

Solving  Eq.  (6.19)  for  R(9),  we  obtain 

2(R^  +  R^)  (  R*{1- R^ +  2Rsme)  +  {R^\-a)  {R^  -  2  R  cosO  +  if  \ 

~  ^  R  \R‘^{R^ -2Rcose  +  l)  +  {a  +  R^\){R^ -2Rcose  +  lfj 

(6.25) 

Eq.  (6.25)  is  evaluated  using  constrained  shooting  technique  described  in  Sec.  5.3.2. 

6.3.1  Local  vs.  Global  Optimality  of  Solutions.  A  fundamental  question  in  the  two 
radar  exposure  minimization  problem  is  whether  the  best  policy  is  to  travel  on  the  Voronoi 
edge,  or  around,  the  two  radars.  Insights  into  this  question  are  obtained  by  examining  the 
relationship  between  locally  and  globally  optimal  solutions. 


6. 3. 1.1  Equal  Power  Radars.  Consider  the  symmetric  exposure  minimiza¬ 
tion  problem  for  two  equal-power  radars,  when  the  points,  A  and  B,  lie  on  the  perpendic¬ 
ular  bisector  of  the  line  joining  the  two  radars.  Cartesian  coordinates  are  now  used,  and 
the  two  radars  are  located  on  the  y-axis,  each  a  distance  a  =  1/2  from  the  origin  -  see,  e.g., 
Eig.  6.9.  In  Appendix  C,  it  is  shown  that  the  perpendicular  bisector  is  a  local  minimizer 
of  the  functional,  defined  by  Eqs.  (6.3)  and  (6.4),  given  that  the  points  A  and  B  lie  on  the 
bisector. 

The  cost  incurred  in  travelling  along  this  line  from  point  A  to  point  B  is 


J|  =  2 


=  2 


r-R 


■f  v^l  +  y'{xf 


'ro  {of  -h 
fRf  1 


dx 


Iro  {a'^  +  x'^y 


dx 
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since  y'{x)  =  0  along  the  perpendicular  bisector  in  our  chosen  coordinate  frame.  Let 

x  =  atan6  dx  =  a sec^  6  d6  (6.26) 


Thus,  the  integral 


2  /■Arctan(%/a) 

=  —  /  cos^  6  dO 


^  Arctan(_Ro/«') 


1 


I'd 


f 


=  ^  /  (1  +  cos  20)  dO 

J-Oo 

=  ^  ^0/  +  0O  +  ^  sin  20/  +  ^  sin  200 


(6.27) 


is  the  total  cost  of  traversing  the  path  from  point  A  to  point  B. 

Now  take  an  alternate  route  “around”  the  radars  by  following  the  path  from  point 
A  to  — oo  to  oo  to  B.  The  cost  is  obtained  by 


T  —  OO 


7^—2 


dx  + 


Rj  (a2+x2)2  (a2  +  x2)2 


dx 
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Using  the  transformation  (6.26),  we  evaluate 


T-oo  ,  jB  _  _ 

-'a  ^  -'oo  -  ^3 


cw/2 

^7r/2 

/  (1  -|-  COS  29)  d9  + 

/  (1-1-  COS  29)  d9 

hf 

ho 

=  ^  ( TT  -  9 f  -  6o  -  Ir  sm29 f  -  sm26o 

Qj  \  Zi 

—  'J  A 


(6.28) 


Hence, the  cost  of  traversing  the  entire  perpendicular  bisector,  i.e.,  from  — oo  to  oo,  is 

+  {Ja^  +  =  73 


Given  a  =  1/2,  if  the  cost  of  a  segment  between  any  points  A  and  B  along  the  perpendicular 
bisector  exceeds  dvr,  the  segment  cannot  be  globally  optimal. 

Lemma  6.3.1.  Consider  the  fully  symmetrie  ease:  9o  =  9f  =  9.  There  exists  a  eritieal 
angle  of  separation,  9c  =  29  ^  47.6535°,  beyond  whieh  the  path  AB  between  the  two  radars 
is  no  longer  globally  optimal. 

Proof.  Given  9o  =  9f  =  9,  the  cost  to  take  the  path  between  the  two  radars,  is 

=  — g  {29  +  sin  29) 


Similarly,  the  cost  to  “go  around”  the  radars  is 

«/a“  +  ^  -  20  -  sin20) 

Equating  and  JJ°°  +  yields 


26»  +  sin26' =  7r/2  (6.29) 

Solving  the  transcendental  Eq.  (6.29)  for  9  is  accomplished  numerically,  and  results  in  the 
unique  solution 

9  ^  23.8268° 
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to  four  significant  digits.  Thus,  the  critical  angle  of  separation  is 


9^  =  29^  47.6535° 


When  9  >  9c,  we  have  thus  the  path  between  the  radars  cannot  be 

globally  optimal.  □ 

Moreover,  the  following  holds. 


Theorem  6.3.1.  For  minimizing  the  exposure  to  two  equal  power  radars,  where  the  start 
and  end  points  of  the  trajeetory  lie  upon  the  perpendieular  biseetor  separating  the  two 
radars,  define  the  ratio 


V  = 


T-oo  I  JB 
^  A  '  ^  OQ 

9f  +  9o  +  ^sm29f  +  ^  sin  20o 
TT  —  9f  —  9o  —  ^  sin  29 f  —  |  sin  29o 


(6.30) 


where  V[9o,9f),  9o,9f  G  [0,7r/2],  is  a  deeision  variable.  If  V{9o,0  f)  >  1,  then  the  globally 
optimal  poliey  is  to  go  around  the  two  radars  and  ifV{9o,0f)  <  1,  then  the  globally  optimal 
poliey  is  to  go  between  the  two  radars. 


Proof.  For  the  case  where  the  start  and  end  points  lie  on  the  perpendicular  bisector,  viz., 
the  line  which  is  the  locus  of  all  equal  power  points,  the  trajectory  that  minimizes  the 
radar  exposure  cannot  deviate  from  that  line.  The  only  two  candidate  minimizers  involve 
following  the  path  from  A  to  B,  i.e.,  “between”  the  two  radars,  or  the  path  from  A  to  oo 
and  from  — oo  to  B,  i.e.,  “around”  the  two  radars.  Thus,  by  definition,  the  function  T> 
discriminates  between  the  globally  optimal  trajectory  and  a  locally  optimal  trajectory.  □ 

Corollary  6.3.1.  For  the  problem  of  Theorem  6.3.1, 


T>  <  1  <  dvr 


Proof.  This  follows  from  the  previous  discussion. 


□ 
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Figure  6.10  Go  Around  or  Go  Between  Decision  Boundary  for  the  Gase  of  Equal  Power 
Radars  and  Endpoints  on  the  Perpendicular  Bisector 

Evaluating  Eq.  (6.30)  over  the  entire  domain  of  9o  and  Of  yields  a  boundary  curve 
that  separates  the  ‘space’  defined  by  {9o,0f)  into  two  distinct  regions.  For  any  {9o,9f) 
pair,  the  decision  to  “go  around”  or  “go  between”  the  two  radars  can  be  made  by  referring 
to  the  curve  V[9o-,  9f),  as  depicted  in  Fig.  6.10.  As  can  be  seen,  for  the  majority  of  problem 
geometries,  it  is  better  to  go  around  the  radars  than  between  them.  Theoretically,  going 
around  the  two  radars  involves  following  a  path  of  infinite  length.  In  practical  terms  air 
vehicle  paths  are  limited  by  mission,  fuel  and  time  constraints.  The  go- around  path  is  still 
very  useful,  as  it  provides  a  lower  bound  of  the  best  achievable  cost. 

Endpoints  Off  the  Voronoi  Edge.  It  is  also  of  interest  to  examine  the 
relationship  between  the  “go  around”  and  “go  between”  paths  for  cases  where  one  or  both 
of  the  endpoints  do  not  lie  on  the  perpendicular  bisector.  Analytic  results  are  difficult  to 
obtain  without  a  closed  form  solution  for  the  extremizing  trajectory. 

Instead,  a  study  of  curves  obtained  from  the  Euler  equation  (6.25)  is  performed.  We 
consider  the  following  geometry  expressed  in  Gartesian  coordinates: 

•  Points  A:  (—0.5,0^),  B:  (0.5, 0.5) 
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•  Radars  :  (0,0)  and  (1,0) 

It  was  observed  that  this  geometry  provides  a  means  to  control  whether  the  solution  of 
the  Euler  equation  (obtained  numerically  by  a  shooting  technique),  went  “around”  or 
“between”  the  two  radars.  When  the  y-axis  coordinate  of  the  point  A  is  positive,  albeit 
infinitesimally  small,  viz.  O"*",  the  solution  to  the  Euler  equation  will  result  in  a  path 
“around”  the  two  radars.  Conversely,  if  the  y-axis  coordinate  of  the  point  A  is  negative, 
yet  albeit  infinitesimally  small,  viz.  0“,  the  solution  to  the  Euler  equation  will  result  in  a 
path  “between”  the  two  radars.  In  order  to  make  the  optimization  problem  well  posed,  a 
constraint  is  required  to  bound  the  path  length,  viz.,  A  >  0  in  Eq.  (6.25). 

Eig.  6.11  shows  some  of  the  paths  “around”  and  “between”  the  radars,  beginning 
with  the  shortest  possible  distance  A  ^  oo  and  then  relaxing  the  path  length  constraint. 
Eor  the  “go  around”  case,  when  A  ^  oo  we  obtain  a  straight  line  connecting  the  points  A 
and  B.  Eor  the  geometry  described  above,  this  results  in  a  cost  approaching  145.782  and 
a  path  length  of  \/2/2.  Alternatively,  for  the  “go  between”  case,  as  A  ^  oo,  we  obtain  two 
straight  line  segments  connecting  the  point  A  with  the  point  B  while  going  through  the 
radar  at  the  origin.  While  the  cost  in  this  case  becomes  infinite,  the  path  length  approaches 

{l  +  V2)/2. 

As  with  the  one  radar  exposure  minimization  problem,  although  the  path  length 
grows  without  bounds,  the  cost  appears  to  be  bounded.  Eig.  6.12  depicts  the  curves  of 
cost  vs.  path  length  for  the  “go  around”  and  “go  between”  trajectories.  While  it  may  seem 
obvious  from  the  geometry  alone,  the  analysis  clearly  shows  that  going  around  the  two 
radars  always  results  in  the  lower  cost,  for  the  given  geometry.  It  is  easy  to  imagine  that 
for  other  geometries,  the  two  curves  may  intersect,  that  is,  there  may  be  cases  when  the 
decision  to  “go  around”  or  “go  between”  is  not  as  obvious. 

6. 3. 1.2  Unequal  Power  Radars.  Consider  the  radar  exposure  minimization 
problem  for  two  unequal  power  radars.  The  radars  are  located  at  {xi,yi),  for  i  =  1,2,  in 
the  xy- plane.  The  power  ratio  is  a  >  0.  The  locus  of  all  points  of  equal  power,  viz., 

a  1 

Rx\d)  ^  R2\B) 
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Figure  6.11  Path  Length  Constrained  Trajectories  Around  and  Between  Two  Radars  for 
A  =  {0.1,1,10,100,^  00} 
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Figure  6.12  Comparison  of  Cost  vs.  Path  Length  for  Constrained  Trajectories  Around 
and  Between  Two  Radars 


from  the  two  radars  is  an  Apollonius  circle,  centered  at  {xa,ya),  with 


Xa 


X2  —  CXl 

1  —  c 


Va  = 


y2  -  cyi 
1  —  c 


and  radius  r, 


r  = 


(6.31) 

(6.32) 


(6.33) 


where 

d  =  \l {xi  -  X2f  +  (2/1  -  2/2)^ 
is  the  distance  between  the  two  radars  and 

1 


is  determined  by  the  two  radar’s  power  ratio. 

We  wish  to  traverse  a  segment  of  the  Apollonius  circle,  from  point  A  to  point  B  as 
shown  in  Fig.  6.13.  The  cost  can  be  expressed  as  the  summation  of  the  cost  from  the  first 
radar  with  the  cost  from  the  second  radar,  i.e.. 


r<^B  yj R,\e)  +  R^\e)  r^B  ^J R2\e)  +  R2\e) 

Ri\e)  Je^  R2H9) 


However,  our  chosen  path  from  A  to  H  lies  on  the  locus  of  equal  power  points.  That  is, 
the  cost  to  traverse  the  path  from  A  to  H  due  to  the  first  radar  is  equivalent  to  the  cost 
imposed  by  the  second  radar.  Thus,  we  let  R{6)  =  R2{9)  and  write 


J  =  2 


Sb  Ji?2(0) 


IdA 


R\e) 


dO 


(6.34) 
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Figure  6.14  Solution  Triangle 

To  calculate  R{9),  we  evaluate  the  solution  triangle  in  Fig.  6.14  where 


7  =  Itt  —  0| 


and  the  distance  between  the  center  of  the  Apollonius  circle  and  the  second  radar  is 


d'  =  \l {Xa  -  X2f  +  iva  “  ^2)' 


(6.35) 


We  have 
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R^{9)  =  —  2  r  d'  cos  7 

=  d'"^  +  —  2r  d!  cos  [vr  —  9\ 

=  d''^  +  r"^  —  2r  d'  cos  (tt  —  9) 

=  d'"^  +  +  2r  d' cos9  (6.36) 


Substituting  Eqs.  (6.31)  and  (6.32)  into  Eq.  (6.35)  yields 


Substituting  Eq.  (6.37)  into  Eq.  (6.36)  we  obtain 

R^{9)  =  ^  +  2  ^ ^j— — ^rdcosd 

Subsituting  Eq.  (6.33)  into  Eq.  (6.38)  yields 

R?{9)=(r^ — 1  - ^ — 2  d^  + -^  ( T~~')  d^cos9 

VI -cy  (i-c)2  Vd\i-cj 

=  ( — ^  d?  l  +  -  +  -^cosd 
V  1  -  c  y  [  C  y/c 

Thus,  we  obtain  the  following  expression  for  R{9), 


(6.37) 


(6.38) 


(6.39) 


(6.40) 
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Differentiating  yields 


yields  an  expression  for  the  cost  to  traverse  an  arc  AB  of  the  Apollonius  circle. 

Without  loss  of  generality,  let  d  =  1  and  consider  9a  G  [0,7r)  and  9b  G  (7r,27r).  We 
define  the  decision  variable 


= 

T27r  jB 

such  that  if  P  <  1 ,  then  the  optimal  strategy  is  to  travel  from  A  to  i?  in  a  counterclockwise 
fashion,  going  between  the  two  radars.  Conversely,  if  H  >  1,  then  the  optimal  strategy  is 
to  go  around  the  two  radars  clockwise. 
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Fig.  6.15  depicts  sample  curves  of  the  decision  variable  T)  vs.  9a^  for  the  case  where 
9b  =  137r/12  and  c  =  {2,  3, 4, 10}.  We  notice  the  curves  of  the  decision  variable  V  appear 
to  be  monotonic,  and  the  break  even  point,  viz.,  where  P  =  1,  is  unique. 


Oa  (deg) 


Figure  6.15  Decision  Variable  for  Two  Radars  of  Unequal  Power,  a  =  2,3,4,10  and 
9b  =  137r/12 

Consider  a  fully  symmetric  problem,  depicted  in  Fig.  6.16,  where  we  require  tt  —  9a  = 
9b  —  For  a  given  a,  we  wish  to  find  some  critical  angle  9c  such  that  the  cost  to  go 
between  the  two  radars  is  equivalent  to  the  cost  to  go  around  the  radars,  viz.,  T)  =  1. 
The  results  are  shown  in  Table.  6.1.  It  is  important  to  note  that  the  angle  9c  is  measured 
from  the  center  of  the  Apollonius  circle.  For  the  case  of  equal  power  radars  mentioned 
previously,  the  critical  angle  was  measured  from  the  origin. 

Table  6.1  Critical  Angles  for  the  Fully  Symmetric  Problem:  Unequal  Power  Radars 


a 

Critical  Angle,  9c 

2 

14.37° 

3 

21.91° 

4 

27.05° 

10 

43.40° 

By  employing  a  root  finding  algorithm,  the  break  even  points  are  calculated  for 
various  values  of  the  power  ratio  parameter  a,  as  9a  and  9b  are  allowed  to  vary  over  their 
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Figure  6.16  Depiction  of  a  Symmetric  Problem  where  D  =  1 

domains.  Fig.  6.17  depicts  such  a  locus  of  break  even  points  for  a  =  2,3,4,10.  From 
studying  Figs.  6.15  and  6.17,  we  see  that  the  strategy  of  going  between  the  two  radars  is 
generally  only  favored  when  Oa  is  close  to  180°.  Values  of  9a  significantly  smaller  than 
180°  requires  0s  to  approach  180°  as  well,  viz.,  the  points  A  and  B  are  close  to  each  other. 
As  expected,  as  9b  — >  360°,  the  break  even  point  occurs  when  9a  — >  180°.  Similarly,  as 
9a  0°  the  break  even  point  is  found  as  0s  — >  180°. 


Ob  (deg) 

Figure  6.17  Locus  of  Break  Even  Points  for  Two  Radars  of  Unequal  Power,  a  =  2, 3, 4, 10 
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6-4  Suboptimal  Technique  for  Multiple  Radar  Avoidance 

In  the  previous  section,  the  optimal  radar  exposure  minimization  trajectory  for  the 
two  radar  case  was  developed.  Extending  this  approach  to  more  than  two  radars  would 
be  impractical.  In  this  section,  an  alternative,  albeit  suboptimal,  method  is  explored  for 
multiple  radar  exposure  minimization. 

One  approach  to  the  multiple  radar  exposure  minimization  problem  is  to  minimize 
exposure  to  the  strongest  (maximum)  single  source  of  radar  energy.  We  consider  the 
multiple  radar  problem  to  be  a  series  of  single  radar  exposure  problems,  viz., 

min  max  \Ji\R(6)]} 

R{e)  Ji[Rm 

such  that 

L[R{e)\  =  [  ^  Vr^  +  de  <  i 

Jo 

where  Ji[R{0)]  is  the  cost  of  the  trajectory  to  avoid  the  radar  and  the  appropriate 
boundary  conditions  are  satisfied. 

Consider  the  following  algorithm: 

1.  Identify  the  radar  with  the  highest  individual  cost 

2.  Plan  a  path  to  avoid  the  radar  which  meets  the  path  length  constraint 

3.  Progress  along  the  planned  trajectory  some  distance  Al 

4.  Repeat  until  the  terminal  point  is  reached 

This  algorithm,  while  suboptimal,  is  easily  implemented  and  can  scale  to  any  number 
of  radars.  However,  for  most  problems,  a  path  length  constraint  will  need  to  be  imposed 
to  keep  the  single  radar  problems  feasible. 

Fig.  6.18  depicts  the  first  three  iterations  of  the  minimax  algorithm.  We  have  two 
equivalent  radars,  one  at  the  origin  and  the  second  at  (1,0).  The  start  point  is  (1/2,- 
1/2)  and  the  end  point  is  (1/2, 1/2).  The  algorithm  begins  by  planning  a  path  to  avoid 
the  radar  at  the  origin.  The  aircraft  traverses  a  distance,  Al  =  0.2,  denoted  by  a  solid 
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line.  The  dotted  lines  indicate  the  remainder  of  the  planned  trajectory.  The  aircraft 
then  determines  that  the  second  radar  has  greater  intensity  and  plans  a  new  route  around 
it.  Fig.  6.18  depicts  the  first  three  segments  of  the  trajectory  generated  by  the  minimax 
technique.  The  Al  parameter  is  chosen  to  be  extremely  coarse  to  illustrate  the  algorithm. 


Figure  6.18  Portion  of  a  Minimax  Generated  Trajectory  for  Two  Radar  Exposure  Min¬ 
imization,  Radars  Located  at  (0,0)  and  (1,0),  Al  =  0.2 


Fig.  6.19  shows  the  same  problem  solved  by  the  minimax  algorithm  using  a  A/  =  0.01. 
We  see  that  the  minimax  algorithm  converges  to  the  Voronoi  edge  separating  the  two 
radars,  as  expected. 

While  the  minimax  approach  provides  a  simple  means  to  deal  with  a  complex  prob¬ 
lem,  there  are  drawbacks  to  this  suboptimal  approach.  In  the  minimax  approach,  solutions 
that  lie  on  the  Voronoi  edges  are  often  preferred.  Consider  the  minimax  algorithm  under 
the  scenario  depicted  in  Fig.  6.20.  Under  this  scenario,  travelling  from  point  A  to  point  B 
would  place  the  aircraft  in  between  a  sequence  of  two  radars  where  the  distance  separating 
the  two  continues  to  decrease.  By  making  sequential  locally  optimal  decisions,  the  globally 
optimal  decision  is  not  considered.  The  preferred  strategy  would  be  to  avoid  such  a  trap. 
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Minimax  Generated  Trajectory  for  Two  Radar  Exposure  Minimization, 
Radars  Located  at  (0,0)  and  (1,0),  Al  =  0.01 
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Figure  6.20  Trapping  the  Minimax  Algorithm 
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6.5  Summary 

In  this  chapter,  several  applications  and  extensions  of  the  fundamental  results  ob¬ 
tained  in  earlier  chapters  were  presented.  A  hierarchical  cooperative  control  law  was  de¬ 
veloped  for  the  problem  of  multiple  vehicle  isochronous  rendezvous  for  radar  exposure 
minimization.  The  algorithm  was  applied  to  a  representative  two-ship  problem.  The  con¬ 
structive  nature  of  the  algorithm  renders  the  existence  of  these  solutions  readily  verifiable. 

The  optimality  of  the  Voronoi  edge,  for  the  problem  of  an  aircraft  exposed  to  illu¬ 
mination  by  two  threat  radars,  has  been  investigated.  An  analytic  solution  to  the  two 
radar  problem  is  presented  for  the  case  where  the  end  points  of  the  trajectory  lie  upon 
the  Voronoi  edge.  The  issue  of  local  versus  global  optimality  is  investigated  and  a  deci¬ 
sion  criterion  is  identified  for  when  going  around  the  two  radars  is  preferable  than  going 
between  the  two  radars.  This  decision  criterion  is  developed  for  the  cases  of  equal  and 
unequal  power  radars.  We  have  shown  that  while  the  Voronoi  edge  is  locally  optimal,  the 
globally  optimal  decision,  in  most  cases,  is  not  to  go  between  the  two  threats.  Thus,  an 
awareness  of  the  consequences  of  embracing  non-globally  optimal  solutions  for  large  scale 
optimization  problems  such  as  the  cooperative  control  of  UAVs  is  essential. 

Finally,  a  suboptimal  technique  is  presented  which  is  scalable  to  address  the  n-radar 
exposure  minimization  problem.  The  next  chapter  develops  optimal  trajectories  for  target 
identification  and  classification  missions  performed  by  autonomous  air  vehicles. 
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VII.  Optimal  Trajectories  for  Autonomous  Target  Classification 

7. 1  Introduction 

A  projected  role  for  autonomous  uninhabited  air  vehicles  is  to  classify  and  subse¬ 
quently  attack  time  critical  targets,  as  well  as  perform  battle  damage  assessment  after  an 
attack  [21,  25].  Thus,  the  problem  of  determining  optimal  look  angles  for  automatic  target 
recognition/classification  is  addressed  first.  Next,  minimum  time  trajectories  for  this  mis¬ 
sion,  for  a  vehicle  with  a  minimum  turning  radius  are  constructed.  Lastly,  an  algorithm 
for  performing  cooperative  classification  and/or  battle  damage  assessment  involving  more 
than  one  air  vehicle,  is  presented. 

1.2  Optimal  Look  Angles  for  ATR 

The  targets’  universe  of  discourse  is  stipulated  to  consist  of  rectangles  located  in  a 
plane  and  with  an  arbitrary  orientation.  We  assume  that  the  probability  of  successfully 
classifying  a  rectangular  target  with  sides  a  and  6,  using  an  Automatic  Target  Recognition 
(ATR)  algorithm,  is  directly  proportional  to  the  projection  of  the  visible  sides  of  the 
rectangular  object  onto  a  line  perpendicular  to  the  aspect  angle  6  of  the  viewing  sensor, 
as  shown  in  Fig.  7.1. 

Without  loss  of  generality  assume  a  >b.  The  probability  of  classification  p(6)  is  then 
calculated  as  follows: 

/ 

a  COS  0+6  sin  0  C\  ^  ^  n 

a+6  ^  ^  ^  2 

—  a  cos  0+6  sin  0  't:  ^  a  ^ 

p{e)  =  \  ^  (7.1) 

—  a  cos  0—6  sin  0  ^  ^  a  37r 

- ^+5 - 

a  COS  0— 6  sin  0  Sn  ^  n  ^ 

- ^  <  27r 

\ 

A  plot  of  p{6)  for  — TT  <  6*  <  TT  is  provided  in  Fig.  7.2. 
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Figure  7.1  Geometry  for  the  Optimal  Look  Angle  ATR  Problem 

Consider  the  function  p{6)  where  9  G  (0, 7r/2).  The  first  and  second  derivatives  of 
p{9)  are 


P{0) 

p\e) 


—a  sin  6  +  6  cos  9 
a  +  b 

—a  cos  9  —  6  sin  9 
a  +  b 


(7.2) 

(7.3) 


Evidently,  we  have  an  extreme  value  of  p  at  the  aspect  angle 


=  Arctan  (  — 
a 


Substituting  Eq.  (7.4)  into  Eq.  (7.3)  yields  the  result 


p"(0‘)  =  <  0 


(7.4) 


i.e.,  we  have  a  relative  maximum  at  the  extreme  point  9  = 
Eq.  (7.1)  finally  yields 

+  6^ 
a  +  b 

This  maximum  probability  is  repeated  at  the  aspect  angles 


9*.  Substituting  Eq.  (7.4)  into 


»•  =  ||  +  Arctan  (1)  ,  -  Arctan  (^)  .  -f  -  Arctan  Q)  | 
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p[e] 


e 


and 


Figure  7.2  Parametric  Plot  of  Probability  of  Classification  for  — tt  <  9  <  tt 
At  9  =  7rj2,  we  have  p  =  However,  p'{7rf2)  does  not  exist.  Since 

p{9)  <  0  for  Arctan  ^  ^ 

p{9)  >0  for^<6<<^  +  Arctan  |  -  j 

^  ^  \(XJ 


we  conclude  pijx j 2^  is  a  relative  minimum.  Similarly,  we  find  that  /)(— 7r/2)  is  a  relative 
minimum.  Lastly,  we  note  that 


7.2.1  Multiple  Look  Classification.  For  multiple-look  classification,  where  a  pre¬ 
determined  probability  of  correct  classification  threshold  needs  to  be  achieved  for  a  target 
to  be  classified,  possibly  employing  more  than  one  air  vehicle,  we  consider  the  probability 
of  correct  classification,  p{9i)  for  a  look  i,  where  the  aspect  angle  is  0*,  to  be  an  indepen¬ 
dent  event.  Thus,  the  probability  p  of  correctly  identifying  a  target  after  n  independent 
snapshots  have  been  taken  is 

n 

^  =  1  _  n  [1  -  p{9i)]  (7.5) 

i=\ 
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and  for  the  special  case  of  two  looks,  n  =  2,  we  have  the  probability  of  correctly  having 
classified  the  target 

P  =  piOi)+p{02)-p{Oi)p{92)  (7.6) 

The  function  p(9)  is  the  probability  of  being  able  to  classify  the  target  when  a 
snapshot  of  the  target  is  taken  from  an  aspect  angle  of  9.  Strictly  speaking,  a  snapshot  of 
the  target  is  taken  and  is  being  reviewed,  after  which  a  binary  (yes/no)  decision  is  made: 

•  The  target  can  be  identified/classified 

•  The  target  cannot  be  identified/classified 

Then,  if  n  snapshots  of  the  target  are  taken  from  aspect  angles  9i,92,  ■  ■  ■  ,9n,  the  probability 
that  the  target  has  been  correctly  identified/classified,  is  given  by  Eq.  (7.5).  Here,  the  n 
snapshots  are  taken,  and  only  thereafter,  are  they  examined. 


7.3  Optimal  Angular  Separation  for  Second  Look 


Consider  the  scenario  where  an  initial  snapshot  of  the  target  has  been  taken  at  an 
unknown  aspect  angle,  9.  We  wish  to  find  the  optimal  change  in  aspect  angle  A  that 
maximizes  the  average  probability  of  identifying  the  target  in  two  passes.  A  directly 
translates  into  a  change  in  the  vector  of  approach  to  the  target.  Without  loss  of  generality, 
assume  9  £  [0,7r/2],  A  G  [0,7r/2]  and  0  <  |A  —  0|  <  |.  Thus,  the  optimization  problem  is 
posed 


^  I  /  +  A)  -  p{9)p{9  +  A))  d9 1 

A  7r/2 

Substituting  Eq.  (7.6)  into  Eq.  (7.7),  we  define  the  following  cost  function 


(7.7) 


J  = 


a  cos  9  +  bsm  9  ^  a  cos(A  +  9)  +  b  sin(  A  +  9) 


a  +  b 


a  +  b 


{a cos  9  +  6sin6*)  (acos(A  +  9)  +  6sin(A  +  9)) 
{a  +  6)^ 


d9 
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Figure  7.3  Example  Cost  Function  for  Second  Target  Identification  Attempt 


Integration  yields 


,  ,  —  +  6^)  TT  cos  A  2  cos  4  ,  A  ,  ,  A 

J(  A)  = - - ({a+b)  .in  -  +  {a-b)  co.  - 


+ 


a?  sin  A  —  2  a  6  cos  A  —  6^  sin  A 


2  (a  +  hy 


Equivalently, 


J(A)  = 


[4  a  6  —  {a?  +  6^)  (vr  —  4)]  cos  A  +  4  (a  +  6)  (a  +  6)  —  ^  (a  —  6)  sin  A] 


4  (a +  6)2 

This  function  is  plotted  in  Fig.  7.3  for  the  aspect  ratio  a/6  =  1/2,  and  0  <  A  < 
We  can  find  the  extreme  value  for  this  function  by  setting 


dJ 

~dK 


=  0 


which  yields 


—2  {a?  —  6^)  cos  A  +  (a^  (vr  —  4)  +  6^  (vr  —  4)  —  4  a  6)  sin  A 

4  (a +  6)2 
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Table  7.1  Optimal  Change  in  Aspect  Angle  (degrees)  for  the  Second  Target  Identification 
Pass  for  Various  Aspect  Ratios 


b 

a 

A*  (deg) 

1 

0 

2 

26.0180 

3 

37.8579 

4 

44.4394 

5 

48.5994 

6 

57.3850 

We  have  an  extreme  value  of  J  at 

A*  =  Arctan  (  — ^ ^ — - — -  )  (7.8) 

l(a2  +  62)(7r-4) -4a6  y  ^  ^ 

Substituting  Eq.  (7.8)  into  the  second  derivative  of  J  yields 


d\j  _  (vr  -  4)  -  4a6]  ^ 

9A2  4  (a +  5)2 

which  holds  for  all  a,  6  >  0.  Thus,  the  extreme  value  obtained  at  A*  is  a  relative  maximum. 
As  an  example,  for  the  case  shown  in  Fig.  7.3  where  the  aspect  ratio  a/b  =  1/2,  the 
maximum  value  of  J  is  achieved  at 


A* 


Arctan 


6 


28  -  57r 


=  26.018° 


It  is  important  to  note  that  A*  exists  in  the  first  quadrant,  i.e.,  0  <  A*  <  7r/2,  only 
if  6  >  a.  For  all  a  >  b,  J(A)  is  maximized  on  the  interval  [0, 7r/2]  at  the  boundary  A*  =  0. 
Thus,  for  any  rectangular  object  with  a  >  b,  the  optimal  change  in  aspect  angle  for  the 
second  look,  given  an  unknown  aspect  angle  0,  is  zero  -  see,  e.g.,  the  illustration  in  Fig.  7.4. 
Furthermore,  for  a  square  object,  i.e.,  a  =  b,  the  optimal  change  in  aspect  angle  is  zero. 
Table  7.3  presents  the  optimal  change  in  aspect  angle  for  varying  aspect  ratios  for  the  case 
b  >  a. 
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b 


A*  =  0 


A*  =  Arctan 


2(l-r^) 


Figure  7.4  Geometry  for  the  Optimal  Look  Angle  ATR  Problem 


7.3.1  Feedback  Control.  If  the  ATR  algorithm  can  provide  an  estimate  of  the 
aspect  angle,  0,  then  the  second  pass  should  be  flown  such  that  the  probability  of  classifi¬ 
cation  of  the  second  pass  will  be  maximized,  i.e., 

6*  =  Arctan 


Thus,  invariably,  the  target  identification  probability  after  the  second  look  is  given  by 
Eq.  (7.4).  Let  the  aspect  ratio  r  =  b/ a.  We  have 


p{e*) 


Vr"^  +  1 
r  +  1 


The  probability  of  classification  after  two  looks  is 


p{9)  =  p{9)  + 


\/r2~+T 
r  -|-  1 


piO) 


\/r2~+T 
r  -|-  1 


1 

r  -|-  1 


\/r2  -|-  l\ 

r  +  l  J 
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P  [r] 


Figure  7.5  Average  Classification  Probability  after  Two  Looks  when  First  Pass  Aspect 
Angle  is  Known 


Thus,  the  average  classification  probability  according  to  the  feedback  strategy  is 


P{r) 


\/r2  +  1 
r  +  1 


\/r2  +  l\ 

r  +  1  j 


p{e) 


dO 


\/r2  +  1  1  (  \/r2  +  1  \  f  2  a  cos  9  +  bsm  9 

TTT  +  f  -  TTT  j  X  - 


Hence,  the  average  probability  of  classification  after  two  looks,  as  a  function  of  the  aspect 
ratio,  is 


P(r)  =  -+[l-- 


TT 


TT 


\/r2  +  1 
r  +  1 


A  plot  of  P(r)  is  shown  in  Fig.  7.5.  Note  that  P(oo)  =  1. 


7-4  Minimum  Time  Trajectories  with  a  Minimum  Turning  Radius  Constraint 

In  this  section  we  consider  minimum  time  optimal  trajectories  where  an  air  vehicle 
travelling  in  the  plane  with  a  constant  velocity  is  constrained  to  have  a  minimum  turning 
radius,  R. 

7.4.1  Specified  Terminal  Point.  The  first  case  examined  entails  a  specified  ter¬ 
minal  point.  The  vehicle  is  initially  at  the  origin,  O  with  some  initial  heading  angle  and 
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Figure  7.6  Minimum  Time  Trajectory  Problem  with  Specified  Terminal  Point  Outside 
the  Minimum  Turning  Radius 

the  terminal  point,  Pf,  is  outside  the  minimum  turning  radius.  Without  loss  of  generality, 
assume  Pf  £  RH  plane,  see,  e.g..  Fig.  7.6. 

Proposition  7.4.1.  The  minimum  time  trajeetory,  for  the  ease  where  the  speeified  termi¬ 
nal  point  Pf  lies  outside  the  minimum  turning  radius  eirele,  eonsists  of  a  hard  turn  into 
Pf  until  the  hearing  to  Pf  is  0°,  followed  by  a  straight  line  dash  to  Pf. 

Proof.  This  problem  can  be  viewed  as  a  special  case  of  a  two  player  differential  game, 
where  the  second  player  is  considered  stationary  and  the  first  player  is  restricted  to  move 
at  constant  velocity  with  a  limited  maximum  turning  radius.  The  problem  described  is 
a  special  case  of  the  Homicidal  Chauffeur  Game,  and  the  minimum  time  trajectory  of  a 
hard  turn  into  Pf  until  the  bearing  to  Pf  is  0°,  followed  by  a  straight  line  dash  to  Pf  is 
the  optimal  strategy  for  the  first  player  -  see,  e.g.,  [22].  □ 

The  second  case  considered  is  the  situation  where  the  final  point,  Pf,  is  inside  the 
minimum  turning  radius  circle.  The  position  of  Pf  is  specified  by  the  distance,  d,  from  the 
center  of  the  left  minimum  turning  radius  circle,  and  by  the  angle  ip  -  see,  e.g..  Fig.  7.7. 

The  parameters  d  and  (p  are  constrained  as  follows.  The  equation  of  the  right  mini¬ 
mum  turning  radius  circle  is 

(x  -  2Rf  +y^  =R^ 
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Figure  7.7  Minimum  Time  Trajectory  Problem  with  Specified  Terminal  Point  Inside  the 
Minimum  Turning  Radius 

The  final  point  Pf  has  coordinates  {x,y)  =  (dcos  99,  dsin  (/j).  Since  Pf  is  inside  the  right 
minimum  turning  radius  circle,  we  have 


{d  cos  (p  —  2  R)"^  +  dF  sin^  <  R? 

d^  —  A  Rd  cos  (^  +  3  <  0 


Solving  for  the  parameter  yields  the  constraint 


max(0,  2  cos  tp  —  yj  4  cos^  f  —  ^)  <  79  <  2  cos  (p  +  \/A~co^Pp^^  (7-9) 

R 


Given  Pf  is  contained  within  the  right  minimum  turning  radius  circle,  the  angle  ip  is 
constrained  by  the  two  lines  tangent  to  the  right  minimum  radius  turning  circle  which 
pass  through  the  origin.  By  constructing  right  triangles  with  the  following  three  points: 
the  origin,  the  center  of  the  right  minimum  turning  radius  circle,  and  each  of  the  points  of 
tangency,  it  is  evident  that 


TT 

6 


TT 


-  -  6 


(7.10) 


In  Fig.  7.8  we  plot  the  constraints  (7.9)  as  a  function  of  ip. 


The  solution  of  the  optimal  control  problem  entails  the  construction  of  a  circle  of 
radius  R  which  is  tangent  to  the  left  minimum  turning  radius  circle  and  which  passes 
through  the  point  Pf.  This  requires  construction  of  a  triangle,  given  three  sides:  d,  R  and 
2  R.  The  construction  is  illustrated  in  Fig.  7.9(a). 
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Figure  7.8  Plot  of  vs  (/? 


Figure  7.9  Construction  of  Circle  Required  to  Solve  the  Minimum  Time  Trajectory  Prob¬ 
lem  with  Specified  Terminal  Point  Inside  the  Minimum  Turning  Radius 


7-11 


The  shown  trajectory  entails  a  hard  turn  to  the  left  followed  by  a  hard  turn  to  the 
right,  viz.,  a  swerve  maneuver.  The  length  of  the  trajectory  is  determined  by  the  angles  of 
the  constructed  triangle,  which  is  shown  in  Fig.  7.9(b).  Note:  Since  d  <3  R,  the  solution 
triangle  can  always  be  constructed.  From  the  law  of  cosines  we  have 


Similarly 


=  A  R?  +  R?  —  A  cos  a 


R^  =  AR^  +  d^  -  ARd  cos  (5 


(3  =  Arccos 


3  ( d 


The  path  length,  Z,  for  the  swerve  maneuver  is 


Hence, 


—  =  /?  +  (/?  + 27r -  a 
R 


I  A  /  1 

—  =  ZTT  +  09  +  Arccos  - 
R  \  4 


Arccos  - 

\  4 


-'I 


(7.11) 


The  construction  shown  in  Fig.  7.9(b)  will  “fail”,  viz.,  the  point  at  which  the  right  hand 
constructed  circle  is  tangent  to  the  left  hand  minimum  radius  circle  will  fall  below  the 
initial  position,  if 

If  <  0 

Similarly,  if 


ip  >  30° 


the  constructed  circle  cannot  simultaneously  contain  points  inside  the  right  minimum  turn¬ 
ing  radius  circle  and  have  a  point  tangent  to  the  left  minimum  turning  radius  circle.  How¬ 
ever,  these  conditions  need  never  occur. 
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Figure  7.10  Plot  of 

Proposition  7.4.2.  All  final  points,  Pf,  in  the  interior  of  the  right  hand  minimum  radius 
eirele  ean  be  reaehed  by  the  eonstrueted  eirele  if  the  inequalities,  0  <  (/?  <  tt/G  and  /?  >  0 
are  satisfied. 

Proof.  Let  0  <  (^  <  tt/G  such  that  A  is  any  point  in  the  interior  of  the  right  minimum 
turning  radius  circle,  some  distance,  d,  from  the  center  of  the  left  minimum  turning  radius 
circle,  O  -  see,  e.g..  Fig.  7.11.  3  a  pair  of  triangles,  each  with  a  common  side  OA  a  side 
of  length  2R  and  a  side  of  length  R.  Let  the  points  of  the  triangles  opposite  OA  be  C 
and  C' .  Construct  two  circles  of  radius  R  at  C  and  C'.  Both  circles  will  be  tangent  to 
the  left  minimum  radius  circle  at  one  point.  The  angle  between  the  point  of  tangency  and 
center  of  the  left  turning  radius  circle  is  (5.  Associate  the  circle  at  C  with  /3  >  0  and  the 
circle  at  C  with  (3  <t).  The  circle  for  which  /?  <  0  must  be  disregarded  since  this  requires 
movement  in  the  opposite  direction  of  the  initial  heading  angle.  Thus,  the  constructed 
circle  with  /3  >  0,  is  the  only  solution  which  intersects  point  A  and  permits  movement  in 
the  direction  of  the  initial  heading  vector. 

□ 

We  also  must  investigate  the  possibility  of  a  hard  turn  to  the  left  followed  by  a 
straight  line  dash,  as  shown  in  Fig.  7.12.  Now, 
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R 


1 

Arccos 

.(l)J 

-  Arccos 

1 

m 

1  R  +  Vd^  + 


+ 


I'- 


(7.12) 


Equating  Eqs.  (7.11)  and  (7.12)  yields  the  following  transcendental  equation,  in  terms 
of  a  single  parameter,  ^ 


Arccos  - 

\  4 


- 


Arccos  - 

\  4 


5-  - 


Arccos 


(I) 


+ 


-1  (7.13) 


A  solution,  of  Eq.  (7.13)  must  also  satisfy  the  inequality  (7.9).  Obviously,  1  < 
^  <  3.  Further,  if  Pf  is  restricted  to  the  interior  of  the  right  minimum  turning  radius 
circle,  then  1  <  ^  <  3  holds. 

Proposition  7.4.3.  The  transcendental  equation  (7.13)  does  not  have  a  solution  which 
satisfies 

Proof.  Suppose  is  a  solution  of  the  transcendental  equation  (7.13).  Evidently,  ip  must 
satisfy  </?  +  /?  <  0,  which  implies 


(p<(p 


Arccos  - 

V  4 


- 


Should  a  solution  of  Eq.  (7.13)  exist,  there  would  be  a  boundary  in  the  right  minimum 
turning  radius  circle  separating  the  optimal  policies  of  two  circular  turns  and  one  circular 
turn  with  a  straight  line  dash  -  see,  e.g..  Fig.  7.13.  Consider  that  the  final  point,  Pf,  lies 
in  the  unshaded  areas  of  the  right  minimum  turning  radius  circle  shown  in  Fig.  7.13.  The 
question  arises,  what  would  then  be  the  optimal  policy  in  those  unshaded  areas?  If  no 
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Figure  7.13  Hypothetical  Boundary  Separating  Two  Optimal  Policies 


Figure  7.14  Hypothetical  Boundary  Separating  Two  Optimal  Policies 


unshaded  areas  exist,  e.g.,  see  Fig.  7.14,  then  we  have  the  case  where 


d 

R 


=  V3 


Inserting  ^  =  a/S  into  Eq.  (7.13),  we  can  check  to  see  if  this  can  be  a  solution.  For  the 
left  hand  side  we  have 


Arccos 


+  Arccos  (  ^ 
.\/3 


For  the  right  hand  side  we  have 


Arccos 


+  \/2 
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Since 


and 


Arccos 


\/2  > 


vr 


6 


>  Arccos 


Arccos 


we  see  that 


^  =  \/3  is  not  a  solution  of  the  transcendental  Eq.  (7.13). 


□ 


Theorem  7.4.1.  The  minimum  time  trajeetory,  for  the  ease  where  the  final  point,  Pf,  is 
inside  the  minimum  turning  radius  eirele,  entails  a  swerve  maneuver  whieh  eonsists  of  a 
hard  turn  to  the  left  of 

(p  +  Arccos 


followed  by  a  hard  turn  to  the  right  of 


2ti  —  Arccos 


Proof.  This  result  proceeds  from  the  previous  discussion. 


□ 


7.5  Minimum  Time  Trajeetories  for  Target  Classifieation 

Consider  a  target  classification  problem  where  the  sensor  footprint  is  circular,  of 
radius  r,  and  it  is  offset  from  the  sensor  platform  by  a  distance  d.  The  sensor  platform 
moves  at  a  constant  velocity  and  has  a  minimum  turning  radius  of  R.  We  define  the 
distance 

u  =  d  +  r 

and  stipulate  that  the  target  circle  defined  by  a  target  at  the  center,  Pf,  with  radius  n,  be 
approached  orthogonally  from  the  outside.  See,  e.g..  Fig.  7.15. 

7.5.1  Type  1  Problems.  Consider  the  case  where  Pf  is  outside  the  circle  of  radius 
z/',  which  is  concentric  with  the  right  minimum  turning  radius  circle. 

P  =  Vi/2  +  i?2 
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Figure  7.15  Minimum  Time  Trajectories  for  Target  Classification 
The  critical  circle  is  shown  in  Fig.  7.16. 

Proposition  7.5.1.  When  Pf  is  outside  the  eritieal  eirele  the  minimum  time  trajeetory 
entails  a  hard  turn  into  Pf  until  a  hearing  of  0°  to  Pf  is  established.  Thereafter,  a  straight 
line  path  is  followed  until  the  target  eirele  is  met,  as  shown  in  Fig.  7.17. 

Proof.  See  proof  of  Proposition  7.4.1.  □ 


7.5.2  Type  2  Problems.  Next  we  consider  the  case  where  Pf  is  inside  the  circle 
of  radius  P  which  is  concentric  with  the  right  minimum  turning  radius  circle,  see,  e.g.. 
Fig.  7.18.  This  yields  the  following  condition 

If  u'  >  2R,  i.e.,  +  R'^  >  AR^,  i.e.,  n  >  \/3R,  then 
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Critical  Circle 


Figure  7.18  Pf  Inside  Critical  Circle 


If  u'  <  2R,  i.e.,  v  <  y/SR,  then 


(  y  \  (  y 

—  Arccos  —  <  (/?  <  Arccos  — 

V  2R )  \ 2R 


,  (  y/y^  +  R^\  ,  ( y/y^  +  R? 

—  Arccos  - — -  <  </?  <  Arccos  - — - 

\  2R  I  \  2R 


(  I  /  y\^  \  (  1  f  y\^ 

^  —  Arccos  +!  <</?<  Arccos  7: \/  I  75  I  +  1 

\  2  V  Vii/  /  V  2  V  yRy 


For  Pf  to  be  inside  the  critical  circle,  we  have 


{x  -  2Rf  +  y‘^  <  y^  =  +  R^ 


with  X  =  dcos(p,  y  =  dsinc/?,  thus 


(dcos  (f  —  2R)^  +  (f  sin^  (p  <  y"^  +  R^ 


cf  —  A  Rd cos  (/?  +  3  <  0 
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d 


Figure  7.19  Plot  of  ^  vs  (p 


We  have  the  following  inequality 


max  0,2cos(/?  — 


(7.14) 


Sinoe 


2  008  (/9 


and 


2  00s  p 


Eq.  (7.14)  beoome8 


4  008^  p  + 


-3  <  0, 


4  OOs2  p  + 


-3  >  0, 


y  u  <  VsR 


y  u  >  VsR 


0  <  ^  <2  00s  p  +  y^4oos2  p  +  —  3, 

2  00s  p  —  y^4oos2  p  +  (;|)^  —  3  <  <2  00s  (/5  +  \J Acos^  p  +  (;|)^  —  3, 


z/  <  \/3i7 
z/  >  ^/^R 


In  Fig.  7.19  we  plot  the  eon8traint  (7.14)  a8  a  funotion  of  p. 

Con8ider  the  8ituation  where  Pf  i8  in8ide  the  oirole  of  radiu8  z/',  whioh  is  oonoentrio 
with  the  right  minimum  turning  radius  oirole,  the  minimum  time  trajeotory  involves  a 
swerve  maneuver.  The  swerve  is  just  enough  for  Pf  to  be  reaohed  by  the  boundary  of  a 
new  oritioal  oirole  of  radius  oonoentrio  with  a  new  right  minimum  turning  radius  oirole. 
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New  Critical  Circle 


Figure  7.20  Construction  of  New  Critical  Circle 


Figure  7.21  Solution  Triangle 

This  is  depicted  in  Fig.  7.20.  We  need  to  construct  a  solution  triangle,  e.g.,  see  Fig.  7.21, 
given  its  three  sides  d,  2R  and  We  calculate 

=  (f +4  R^  —  4Rd  cos  /3 


Thus, 


,/2 


cos  (3  = 


(P  +  AR'^-u^ 
iRd 
_  3R^  + 

~  iRd 

6  =  Arccos  I  - 

\  4 


'd\ 

kR) 


+  ^  - 


'  \  2 
kR) 

ii) 
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Similarly, 


(f  =  1/'^  +  4  —  4  Ru'  cos  /3 


Thus, 


cos  a  = 


u'^  +  4R^-£ 
4Ru' 

5i?2  +  z/2  _  (^2 


c^^2 


5-(i) +(^r 


/3  =  Arccos  - 

'  4 


a/i  +  (1)^ 


For  the  turn  to  the  left,  the  change  in  heading  is 


(p  +  P 


For  the  turn  to  the  right,  the  change  in  heading  is 


27r  —  a  —  Arctan  (  — 
^  R 


The  path  length  is  then  expressed  as 


—  =  P  +  ip  +  2TT 


a  +  Arctan  (  — 
V  R 


Equivalently, 


—  =  2tt  +  p  —  Arctan  ( — 'j  +  Arccos  (  - 
R  \RJ  V  4 


3 


ii)  ii) 


—  Arccos  I  - 


a/i + i^y 


Remark  7.5.1.  The  target  eirele  is  penetrated  orthogonally  by  eonstruetion. 
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Figure  7.22  Two  Candidate  Minimum  Time  Trajectories  After  An  Initial  Swerve 

Theorem  7.5.1.  When  Pj  is  inside  the  eirele  of  radius  v' ,  whieh  is  eoneentrie  with  the 
right  minimum  turning  radius  eirele,  the  minimum  time  trajeetory  involves  a  swerve  ma¬ 
neuver.  The  swerve  is  just  enough  for  Pj  to  be  reaehed  by  the  boundary  of  a  new  eritieal 
eirele  of  radius  v' ,  eoneentrie  with  a  new  right  minimum  turning  radius  eirele. 

Proof.  Consider  the  situation  after  an  initial  swerve  which  puts  Pf  on  the  boundary  of  a 
circle  of  radius  n' ,  centered  at  the  new  position  of  the  right  minimum  turning  radius  circle. 

From  this  time  onward,  the  “optimal”  trajectories,  (a)  and  (b),  are  indicated  in 
Fig.  7.22.  Similarly  to  the  proof  in  Proposition  7.4.1,  we  realize  that  the  optimal  trajectory 
is  (a).  □ 

7.6  Cooperative  Target  Classifieation 

Lastly,  consider  the  problem  of  cooperative  target  classification.  As  in  Section  7.5, 
the  sensor  footprint  is  circular,  of  radius  r,  and  offset  from  the  sensor  platform  by  a  distance 
d.  Additionally,  there  is  a  commanded  approach  angle  which  corresponds,  for  instance, 
to  an  optimal  second  look  at  the  target  directed  by  the  results  of  Sec.  7.3.  Since  our  target 
is  assumed  rectangular,  symmetry  provides  four  candidate  angles  of  approach,  as  shown 
in  Fig.  7.23. 
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Figure  7.23  Cooperative  Target  Classification  Problem 


The  first  step  in  developing  the  optimal  trajectory  to  perform  the  cooperative  target 
classification  is  to  construct  minimum  turning  radius  circles  such  that  the  entry  into  the 
target  circle  is  orthogonal.  Eight  circles  are  constructed,  see,  e.g..  Fig.  7.24. 

All  candidate  trajectories  consist  of  a  hard  turn,  a  straight  line  dash  and  a  hard 
turn.  The  straight  line  dash  is  to  be  constructed  tangent  to  the  two  minimum  turning 
radius  circles.  Depending  on  the  location  of  Pf  and  the  sizes  of  the  target  circle  and  the 
minimum  turning  radii  circles,  a  swerve  maneuver  may  be  considered  as  a  portion  of  a 
candidate  trajectory.  Fig.  7.25  depicts  several  candidate  minimum  time  trajectories. 

Lastly,  the  candidate  trajectories  are  ranked  according  to  path  length,  and  the  min¬ 
imum  length  trajectory  is  selected  as  optimal. 

7. 7  Summary 

In  this  chapter,  a  cost  function  was  developed  for  enhancing  the  probability  of  au¬ 
tomatic  target  recognition  and  the  corresponding  optimization  problem  was  solved.  This 
requires  approaching  the  target  from  specified  directions  and  making  sure  that  the  air 
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Figure  7.25 


Construction  of  Some  Candidate  Minimum  Time  Trajectories  for  Coopera¬ 
tive  Target  Classification 
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Figure  7.26  Characterization  of  the  Minimum  Time  Cooperative  Classification  Trajec¬ 
tory 

vehicle’s  sensor  covers  the  target.  Minimum  time  trajectories  for  air  vehicles  with  a  mini¬ 
mum  turning  radius  were  constructed  for  these  tasks,  and  also  for  target  attack  and  battle 
damage  assessment.  Lastly,  the  single  vehicle  results  are  encapsulated  in  a  methodology 
addressing  cooperative  target  classification  and  attack. 
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VIII.  Conclusion  and  Recommendations 


8. 1  Introduction 

The  objective  of  this  research  was  to  develop  a  mathematically  rigorous  approach  to 
a  problem  in  air  vehicle  path  planning,  seeking  fundamental  truths  concerning  autonomous 
air  vehicles  and  their  cooperative  control.  The  objective  was  motivated  by  the  intention 
to  develop  a  significant  military  force  consisting  of  UAVs,  UCAVs  and  autonomous  mu¬ 
nitions  and  the  desire  for  this  technology  to  make  an  effective  contribution  to  air  power. 
Specifically,  this  research  considered  the  air  vehicle  path  planning  problem  in  the  context 
of  radar  exposure  minimization. 

In  Chapter  1  of  this  dissertation,  an  overview  of  air  vehicle  path  planning  was  pre¬ 
sented,  along  with  historical  background,  problem  statement  and  key  results.  Chapter  2 
provided  the  essential  elements  of  radar  technology  and  optimization  theory  utilized  in  this 
research.  In  Chapter  3,  the  single  vehicle  radar  exposure  minimization  problem  was  ad¬ 
dressed,  including  derivations  of  the  optimal  trajectory,  it’s  cost  and  path  length.  A  similar 
analysis  of  the  passive  sensor  exposure  minimization  problem  was  delineated  in  Chapter 
4.  Chapter  5  presented  a  description  of  the  numerical  methods  developed  and  employed 
in  this  research.  Applications  and  extensions  of  the  work  in  the  previous  chapters  were 
presented  in  Chapter  6.  A  hierarchical  cooperative  control  algorithm  was  formulated  for 
the  problem  of  isochronous  rendezvous  of  multiple  air  vehicles.  Furthermore,  the  problem 
of  a  single  vehicle  minimizing  exposure  to  two  radars  was  studied.  In  Chapter  7,  the  pro¬ 
cess  of  automatic  target  recognition  was  modelled.  An  optimal  minimum  time  control  law 
was  developed  for  the  problem  of  autonomous  target  classification  and  was  extended  to 
address  cooperative  target  classification.  In  the  present  chapter,  the  significant  advances 
of  this  research  are  summarized  and  and  suggestions  for  further  research  are  provided. 

8.2  Summary  of  Results 

The  analysis  conducted  in  this  research  has  led  to  the  following  key  results: 

•  The  problem  of  determining  the  flight  path  connecting  the  point  of  departure  and  the 

point  of  arrival,  such  that  the  exposure  of  an  aircraft  to  illumination  by  a  tracking 
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radar  located  at  the  origin  is  minimized,  has  been  solved,  in  closed  form,  using  the 
Calculus  of  Variations. 

—  The  solution  is  shown  to  satisfy  necessary  and  sufficient  conditions  for  a  weak 
local  minimum  as  well  as  the  Weierstrass  necessary  conditions  for  a  strong  local 
minimum. 

—  The  solution  is  shown  to  exist  if,  and  only  if,  the  angle  included  between  the 
departure  and  destination  points  is  less  than  60°. 

—  A  closed  form  expression  for  the  optimal  heading  angle  control  law  is  obtained. 

—  A  closed  form  expression  for  the  optimal  cost  is  obtained. 

—  An  expression  for  the  optimal  path  length,  not  available  in  closed  form,  is  derived 
using  elliptic  integrals. 

—  A  method  for  determining  optimal  path  length  constrained  solutions  to  the 
single  vehicle  radar  exposure  minimization  problem  is  presented  and  conditions 
for  the  existence  of  these  solutions  are  identified. 

•  A  similar  analysis  is  conducted  for  the  single  vehicle  passive  sensor  exposure  mini¬ 
mization  problem  to  include  a  closed  form  solution  and  closed  form  expressions  for 
the  optimal  path  length  and  cost.  A  comparison  between  the  solutions  to  the  passive 
sensor  and  radar  exposure  minimization  problems  is  made. 

•  It  is  verified  that  the  Voronoi  edge  is  the  locally  optimal  solution  for  exposure  mini¬ 
mization  against  two  radars.  The  issues  of  global  versus  local  optimality  are  explored 
yielding  analytic  results  identifying  the  conditions  where  going  around  (between)  two 
radars  is  preferable  to  going  between  (around)  two  radars. 

—  A  suboptimal  algorithm  for  n-radar  exposure  minimization  is  developed  and 
examined. 

•  A  hierarchical  cooperative  control  algorithm  is  formulated  to  determine  optimal 
trajectories  minimizing  radar  exposure  for  two  (or  more)  air  vehicles  performing 
isochronous  rendezvous.  The  constructive  nature  of  the  algorithm  renders  the  exis¬ 
tence  of  these  solutions  readily  verifiable. 
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•  A  methodology  for  designing  optimal  trajectories  for  an  autonomous  air  vehicle  is 
presented.  The  automatic  target  recognition  process  is  formulated  as  a  novel  op¬ 
timization  problem  and  solved.  Minimum  time  trajectories  for  air  vehicles  with  a 
minimum  turning  radius  were  covered.  Finally,  the  problem  of  target  classification 
was  formulated  using  these  concepts  and  a  solution  to  the  autonomous  air  vehicle 
search  problem  is  presented. 

8.3  Recommendations  for  Future  Research 

The  research  in  this  dissertation  has  revealed  several  mathematical  facts  concerning 
the  radar  exposure  minimization  problem.  While  this  research  has  answered  many  impor¬ 
tant  questions,  many  other  challenges  have  been  uncovered  or  remain  unanswered.  These 
issues  are  outlined  in  the  following  paragraphs. 

•  In  this  work,  a  particular  deterministic  metric  was  utilized  to  measure  radar  exposure. 
Future  work  could  consider  other  radar  metrics. 

—  Probability  based  radar  metrics  such  as  probability  of  detection  and  probability 
of  tracking  could  be  developed. 

—  A  radar  cross  section  model  could  be  added  to  the  cost  functional  in  Eq.  (3.1) 
such  that  radar  cross  section  becomes  a  function  of  the  air  vehicle  attitude.  A 
rigorous  mathematical  analysis  would  complement  previous  work  in  this  area 
by  Pendelton  [34]. 

•  Consider  the  problem  of  minimizing  exposure  to  a  radar  where  the  final  point  is  not 
fixed  but  lies  on  some  curve  7.  For  example,  one  may  wish  to  flee  from  “hostile” 
territory  into  “friendly”  territory  where  the  border  is  defined  by  a  line  segment  -  see, 
e.g..  Fig.  8.1.  Alternatively,  one  may  consider  the  lethality  zone  of  an  air  defense 
radar  to  be  a  circle  and  it  is  desired  to  escape  from  inside  the  circle  while  minimizing 
exposure  to  the  radar. 

Similarly,  one  might  consider  the  variable  end  point  problem  as  a  means  to  address 
the  multiple  radar  exposure  minimization  problem.  In  this  case,  the  target  curve,  7, 
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Figure  8.1  Application  of  Single  Radar  Exposure  Minimization  with  Variable  End  Point 

is  a  Voronoi  edge,  and  the  problem  is  to  minimize  exposure  to  a  single  radar.  Both 
McLain  [28]  and  Novy  [30]  presented  ad-hoc  methods  of  getting  onto  the  Voronoi  edge 
and  transitioning  between  Voronoi  edges.  An  open  question  remains  as  to  what  is  a 
computationally  efficient,  possibly  suboptimal,  approach  to  multiple  radar  exposure 
minimization.  Using  a  variational  approach  and  applying  a  transversality  condition 
may  result  in  a  closed  form  solution,  thus  being  extremely  computationally  efficient. 
However,  satisfying  the  transversality  condition  may  result  in  the  extremal  curve 
and  the  Voronoi  edge  meeting  orthogonally.  In  this  case,  an  inequality  constraint 
imposed  with  slack  variables  could  result  in  a  solution  that  approaches  the  Voronoi 
edge  tangentially  as  opposed  to  orthogonally  -  see,  e.g.,  Eig.  8.2.  In  the  end,  a 
comparison  could  be  made  between  these  ad-hoc  methods  and  the  optimal  method 
presented  in  this  dissertation. 


Figure  8.2  Variable  End  Point  Problems  as  a  Means  to  Address  Multiple  Radar  Exposure 
Minimization 

•  It  is  possible  that  a  solution  for  the  path  length  constrained  radar  exposure  mini¬ 
mization  problem  can  be  obtained,  i.e.,  in  terms  of  elliptic  integrals,  similar  to  what 
was  accomplished  for  the  constrained  passive  sensor  problem  in  Sec.  4.6.  Such  a 
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formulation  would  be  a  great  improvement  over  the  numerically  challenging  shooting 
method  employed  in  this  research. 

•  A  numerical  study  could  be  undertaken  to  determine  the  “go-around”  versus  “go- 
between”  decision  boundary  for  the  two  radar  exposure  minimization  problem  for 
the  case  where  the  start  (end)  point  is  not  on  the  Voronoi  edge. 
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Appendix  A.  Characterizing  the  Radar  Exposure  Minimization  Extremal 

In  this  Appendix  we  show  that  the  extremal  identified  in  (3.4)  satisfies  the  necessary  and 
sufficient  conditions  for  a  weak  local  minimum.  Recall  that  our  cost  function  is  of  the  form 


Cf 

J  =  /  F{e,  R,  R)  de 

Jo 


where 


F{9,R,R) 


R4 


(A.l) 


We  define  the  partial  derivatives  we  will  need  to  evaluate  the  conditions  for  weak  and 
strong  extrema: 
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(A.2b) 


(A.2c) 


(A.2d) 


(A.2e) 


The  following  conditions  are  sufficient  for  the  cost  functional  (A.l)  to  have  a  weak 
local  minimum  for  R  =  R{0)  [16]. 


1.  The  curve  R  =  R{0)  is  an  extremal,  i.e.,  it  satisfies  the  Euler  equation. 

2.  The  strengthened  Legendre  condition 

PiO)^^-F^^{e,R,R)>0  (A.3) 

holds  along  the  curve  R  =  R{6). 

3.  The  strengthened  Jacobi  condition  holds,  i.e.,  the  interval  [a,b]  contains  no  points 
conjugate  to  a. 
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We  will  show  these  three  sufficient  conditions  are  met  by  examining  the  second 
variation  of  J,  which  can  be  expressed  as 

f-b 

5^J{h)  =  /  {Ph^  +  Qh'^)  d9  (A.4) 

J  a 

where  h{9)  is  an  increment  of  R{9)  satisfying  h{a)  =  h(b)  =  0  and 

m  = 
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2ii2(^2 
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2R^{R^  +  R^)I 
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In  section  3.3  we  showed  that  R{9)  extremizes  our  cost  functional.  Thus  we  have 
met  the  first  sufficient  condition  for  a  weak  local  minimum. 

The  Legendre  necessary  condition  for  a  weak  minimum  is  that  P{9)  >  0.  Clearly 
this  is  satisfied  in  (A. 5)  for  any  extremal  satisfying  the  Euler  equation  since  R^  >  0  and 
R^  >  0.  Furthermore,  the  strengthened  Legendre  condition  P{9)  >  0  is  met. 

Now  we  examine  the  Jacobi  condition  to  show  there  are  no  points  conjugate  to  0  in 
the  interval  [0,0/].  If  we  can  show  that  the  quadratic  functional  (A.4)  is  nonnegative  for 
all  h{9)  where  P{9)  >  0,  then  the  interval  [a,  b]  contains  no  points  conjugate  to  a  [16].  It  is 
then  sufficient  (albeit  conservative)  for  the  second  variation  to  be  nonnegative,  if  Q{9)  >  0. 

We  wish  to  show  that  Q{9)  >  0.  Clearly,  the  denominator  of  Q{9)  is  strictly  positive, 
since  >  0  and  R^  +  R^  >  0.  Similarly  the  first  three  terms  of  the  numerator  of  Q{9) 
are  positive,  that  is  12R'^  +  15R^R^  +  GR"^  >  0.  Thus,  for  Q{9)  >  0  we  need  only  show 
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Proof. 
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Since  <5(0)  >  0  we  have  demonstrated  all  sufficient  conditions  for  our  extremal  (3.4)  to  be 
a  weak  local  minimum.  □ 

To  be  considered  a  strong  local  minimum,  it  is  necessary  that  the  Weierstrass  E- 
function  condition  be  satisfied,  that  is 

F{e,  R,  R*)  -  F{e,  R,  R)  -  F^{e,  R,  R){R*  -R)>0  (A.7) 

For  our  cost  functional,  the  Weierstrass  E-function  (A.7)  becomes 

\/r*^  +  R‘^  _  \/r^  +  i?2  ^  ^(i7*  -  R) 

+  R^ 
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Obtaining  a  common  denominator 


(A.8) 

Notice  that  the  quantities  {RR*)"^ ,{RR)^ ,{R* and  are  all  guaranteed  to  be 
non-negative.  In  fact,  the  only  quantity  that  could  be  negative  in  (A.8)  is  RR* .  Thus  we 
can  neglect  the  denominator  of  (A.8)  and  to  show  that  the  condition  (A. 7)  holds,  we  need 
only  show  that  the  numerator  of  (A.8)  is  greater  than  or  equal  to  zero.  That  is,  we  wish 
to  show 
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Proof. 
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Thus  we  have  satisfied  the  Weierstrass  E-function  condition,  a  necessary  condition  for  the 
extremal  (3.4)  to  have  a  strong  local  minimum.  □ 
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Appendix  B.  Radar  Exposure  Maximization 

B.  1  Introduction 

In  this  Appendix,  the  problem  of  maximizing  exposure  to  radar  is  examined.  From 
an  operational  standpoint,  the  idea  is  nonsensical.  However,  analysis  of  the  maximization 
problem  provides  insight  into  the  dual  problem  of  radar  exposure  minimization. 

Given  a  radar  at  the  origin,  we  wish  to  find  the  planar  trajectory,  R*{9),  that  max¬ 
imizes  the  reflected  RF  energy  reflected  from  the  target.  The  trajectory  connects  two 
prespecified  points  A  and  R,  respectively  located  a  distance  Rq  and  Rf  from  the  origin, 
and  separated  by  an  included  angle  Of  -  see,  e.g..  Fig.  B.l. 


Figure  B.l  Radar  Exposure  Maximization 


B.2  Unconstrained  Radar  Exposure  Maximization 

As  in  Chapter  III,  the  problem  is  posed  in  the  Calculus  of  Variations.  We  wish  to 
maximize  the  cost  functional 


^  VR^  +  R2 


(B.l) 
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Without  loss  of  generality,  we  can  write  the  boundary  conditions  as 


R(0)  =Ro  =  l 

(B.2) 

R{ef)  =  Rf>i,  o<e<ef. 

(B.3) 

Remark  B.2.1.  For  the  unconstrained  radar  exposure  maximization  problem,  an  optimal 
solution  does  not  exist. 

Proof.  Starting  at  point  A,  any  trajectory  that  maximizes  J  will  run  to  the  origin  O,  and 
then  to  the  point  B.  Hence,  J*  ^  oo  and  there  does  not  exist  an  optimal  solution  to  the 
unconstrained  radar  exposure  maximization  problem.  □ 

B.3  Constrained  Radar  Exposure  Maximization 

Assume^  Rf  =  1.  In  order  to  render  the  maximization  problem  well  posed,  we 
introduce  the  isoperimetric  constraint 

+  R^  de  =  l  (B.4) 

Remark  B.3.1.  For  the  constrained  radar  exposure  maximization  problem,  an  optimal 
solution  does  not  exist  for  I  >  2  or  I  <  2sin(0j/2). 

Proof.  If  /  >  2,  maximizing  J  requires  running  to  the  origin  O.  Hence,  as  in  Remark  B.2.1, 
an  optimal  solution  does  not  exist. 

The  Euclidean  distance  between  the  points  A  and  B  is  given  by 

R"^  +  Rj,  —  2  Ro  Rf  cos  9f 

Thus,  for  Ro  =  Rf  =  I,  A  I  <  2sin(0j/2),  the  isoperimetric  constraint  cannot  be  satisfied 
and  an  optimal  solution  does  not  exist.  □ 

'^This  restriction  is  for  mathematical  convenience  and  can  be  relaxed. 
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Now  consider  the  case  where  2sin(0j/2)  <  I  <  2.  The  constrained  maximization 
problem  described  by  Eqs.  (B.l)  and  (B.4),  can  be  reformulated  as  an  unconstrained 
problem  using  the  Lagrange  multiplier  method.  That  is,  we  form  the  augmented  cost 
functional 

Ja  =  de  (B.5) 

and  choose  a  value  for  A  that  results  in  the  desired  path  length  1.  Clearly,  if  we  choose  a 
value  of  A  such  that  the  trajectory  runs  to  the  origin  and  Ja  ^  oo,  there  does  not  exist 
an  optimal  solution.  Thus,  questions  arise  as  to  whether  the  Lagrange  multiplier  method 
can  be  used  at  all  in  this  maximization  problem  and  how  should  one  solve  the  constrained 
maximization  problem  given  by 
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(B.6b) 

(B.6c) 

(B.6d) 


We  note  that  there  does  not  exist  a  global  maximum  in  the  unconstrained  optimization 
problem  with  the  Lagrange  multiplier,  e.g.,  Eq.  (B.5).  This  does  not  forego  the  possibility 
of  a  local  maximum.  Interestingly,  it  is  precisely  such  a  local  maximum  in  the  unconstrained 
maximization  problem  with  a  Lagrange  multiplier  that  provides  a  globally  optimal  solution 
to  the  constrained  maximization  problem,  Eqs.  (B.6a)-(B.6d). 

Furthermore,  a  range  of  Ac  >  A  >  0  can  exist  where  a  local  maximizing  solution  to 
the  unconstrained  maximization  problem  with  Lagrange  multipliers  does  not  exist.  Here 
lambduc  will  be  determined  by  the  Legendre  necessary  condition  for  a  local  maximum 


^R,R  —  ^ 


Thus  we  require  the  condition 


1 

R^{e) 
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be  satisfied  in  addition  to  our  restrictions  on  I  for  a  maximizing  solution  to  exist.  As  these 
are  necessary  conditions,  there  may  exist  additional  constraints  on  A  and  thus  I  as  well. 

Fig.  B.2  is  a  radar  exposure  maximizing  trajectory  obtained  using  the  Lagrange 
multiplier  method  for  Rf/Ro  =  1,  9f  =  45°  and  A  =  —8.9.  The  resultant  path  length  is 
approximately  0.8418. 


Figure  B.2  Example  Trajectory  for  Radar  Exposure  Maximization 


B.4  Summary 

We  see  that  the  radar  exposure  maximization  problem  is  somewhat  simpler  than 
the  radar  exposure  minimization  problem  of  Chapter  III.  In  the  radar  exposure  mini¬ 
mization  problem,  an  unconstrained  solution  exists  for  O<0/<7r/3  and  a  path  length 
constraint  must  be  imposed  for  Of  >  n/ 3.  In  the  maximization  problem,  an  unconstrained 
solution  does  not  exist  and  only  constrained  solutions  are  sought.  In  addition  to  the  Leg¬ 
endre  necessary  condition  condition  on  the  Lagrange  multiplier,  the  path  length  constraint 
2sin(0j/2)  <l  <2  may  further  constrain  the  range  of  A  we  can  choose. 
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Appendix  C.  Local  Optimality  of  the  Voronoi  Edge  in  Two  Radar  Exposure 

Minimization  Problems 

C.  1  Introduction 

We  are  concerned  with  demonstrating  the  local  optimality  of  a  segment  of  the  Voronoi 
edge  in  the  two  radar  exposure  minimization  problem.  The  Voronoi  edge,  for  this  problem, 
is  the  locus  of  all  points  of  equal  power  separating  the  two  radars.  Specifically,  for  the  case 
of  two  equal  power  radars,  the  Voronoi  edge  is  a  straight  line.  In  this  Appendix,  we  show 
that  the  Voronoi  edge  is  an  extremizing,  and  locally  minimizing^,  path  for  the  two  radar 
exposure  minimization  cost  functional. 

C.2  Equal  Power  Radars 

Consider  the  two  radar  exposure  minimization  problem  where  the  radars  are  located, 
without  loss  of  generality,  on  the  x-axis  at  (0,0)  and  (1,0).  Polar  coordinates  are  used. 
For  the  case  of  equal  power  radars,  the  Voronoi  edge  is  the  perpendicular  bisector  of  the 
line  connecting  these  two  radars  -  see,  e.g..  Fig.  C.l. 

Let  the  initial  point  A  on  the  Voronoi  edge  be  at  a  distance  Rq  from  the  origin,  such 

that 

i?o  =  i?(0o)  =  ^sec0o  (C.l) 

Similarly  for  the  final  point  B  we  have 

Rf  =  Ri.0f)  =  \^ecef  (C.2) 

Thus,  for  6  £  (— 7r/2, 7r/2),  the  equation  of  the  candidate  extremal  trajectory  connecting 
A  and  B,  viz.,  a  segment  of  the  Voronoi  edge,  is 

R{e)  =  ^ sec 0,  eo<e<9f  (c.s) 

'^Conditions  concerning  the  global  optimality  of  segments  of  the  Voronoi  edge  are  addressed  in  Sec.  6.3. 
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Figure  C.l  Two  Radar  Exposure  Minimization  and  the  Voronoi  Edge 

Lemma  C.2.1.  Any  segment  of  the  Voronoi  edge  extremizes  the  eost  funetional  assoeiated 
with  the  (equal  power)  two  radar  exposure  minimization  problem. 


Proof.  Let  AB  be  any  segment  of  the  Voronoi  edge,  Eq.  (C.3)  with  endpoints  given  by 
Eqs.  (C.l)  and  (C.2).  Erom  Sec.  6.3,  the  Euler  equation  for  the  unconstrained  two  radar 
exposure  minimization  problem  is 


1  R 

2  i?2  + 


(r  +  r) 


—  2a  R  cos  9  +  a 


(i?^  +  —  2  a  i?  cos  9) " 

R4 


+  R^  —  —  2aR  sin  0  +  a 


(R^  +  a'^  —  2aR  cos 0) " 

R4 


0  (C.4) 


Let  a  =  1  and  a  =  1.  Differentiating  Eq.  (C.3),  we  have 

R{9)  =  -  sec  9  tan  9  (C.5) 

R{9)  =  I  (3  —  cos  29)  sec^  9  (C.6) 


Substituting  the  candidate  extremal,  Eq.  (C.3),  and  its  derivatives,  Eqs.  (C.5)  and  (C.6), 
into  Eq.  (C.4)  satisfies  the  Euler  equation.  Thus,  any  segment  of  the  Voronoi  edge  is  an 
extremizing  path.  □ 
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Lemma  C.2.2.  Any  segment  of  the  Voronoi  edge  satisfies  the  Legendre  neeessary  eondition 
for  a  weak  heal  minimum,  viz., 


m  =  ^Frr  >  0 


holds  at  every  point  along  the  trajeetory  R{6)  given  by  Eq.  (C.3). 
Proof.  From  Sec.  6.3,  we  have  the  following: 


F-  ■  — 
^RR  ~ 


R^  +  R^ 


where  the  integrand  of  the  cost  functional  is  given  by 

F  =  gVr"^  + 


(C.7) 


(C.9) 


and 

Q  —  ^  _| _ \ _ 

R^{0)  {R^{0) +  a^ -2aR{e)cosef 

Substituting  Eqs.  (C.3)  and  (C.5)  into  Eq.  (C.IO)  and  Eq.  (C.9),  letting  a 
we  have 


(C.IO) 


1  and  a  =  1, 


G  =  32  cos^  6 
F  =  16cos^  6 

Similarly,  we  find 

P  =  32cos®0  (C.ll) 

Thus,  the  Legendre  necessary  condition  \F^p^  >  0  is  satisfied.  □ 

Indeed,  for  any  closed  interval  6  G  [9o,0f],  where  6o  >  — 7r/2  and  Of  <  7r/2,  the 
strengthened  Legendre  condition  \Ff^f^  >  0  is  satisfied. 
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Lastly,  we  show  that  the  strengthened  Jacobi  condition  holds,  namely,  that  the  in¬ 
terval  [6o,9f],  where  6o  >  — 7r/2  and  Of  <  7r/2,  contains  no  points  conjugate  to  the  point 


Oo. 

The  Jacobi  equation,  Eq.  (2.7),  is 

-^(Ph')+Qh  =  0  (C.12) 

du 

where  the  function  P{9)  >  0  is  given  by  Eq.  (C.7)  and 

0  =  (3(9)  =  (C.13) 


Evaluating  the  partial  derivatives  and  substituting  Eqs.  (C.3),  (C.5)  and  (C.6)  into  Eq.  (C.13), 
yields 

Q  =  16  cos®  0  (11 -|- 31  cos  20)  (C.14) 

Substituting  Eqs.  (C.ll)  and  (C.14)  into  the  Jacobi  equation,  we  have 


-16  cos®  0  [2  cos^  0  h''{9)  —  8  sin  29  h' (9)  —  (11  -|-  31  cos  29)  h{9)]  =  0 
2  cos^  9  h"{9)  —  8  sin  20  h'{9)  —  (11  -|-  31  cos  20)  h{9)  =  0 


2  cos^  0 


h'V)  -  H'(«)  -  ^ 

cos^  0  2  cos^  0 

h" {9)  —  8  tan  0  h' {9)  —  (31  —  10  sec^  0)  h{9)  =  0 
h"{9)  —  8  tan  0  h' {9)  +  (lO  tan^  0  —  2l)  h{9)  =  0 


(C.15) 


This  second  order  linear  differential  equation  is  of  the  form 

h”  -|-  a  tan  9  h'  +  [a  tan^  0  +  7)  h  =  0 


which  is  a  special  case  of  the  equation 


h”  +  {a tanO  +  b cot  0)  /i'  +  (a tan^  9  +  P cot^  ^  +  7)  h  =  0 


(C.16) 
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where 


a  =  —8 
6  =  0 
a  =  10 
/3  =  0 
7  =  -21 

Also,  h{9o)  =  h{6f)  =  0.  According  to  [35],  Eq.  (C.16)  has  an  exact  solution.  Following 
the  methodology  in  [35],  pp.  188,  Sec.  2.1.6,  Eq.  (55),  we  employ  the  transformation 


^  =  sin^  0 
h  =  w  tan  0  sec^  0 


Applying  the  transformation  above  leads  to  the  equation: 

c  (C  - 1)  w"{0  +  (?  -  ^)  ^'(0  +  ^  HO  =  0 

with  boundary  conditions  w{^i)  =  0,  rc(^2)  =  0,  G  (0,1)-  Note:  we  assumed  that 
xi  =  sin^  9o  and  X2  =  sin^  0f.  Hence,  this  transformation  lead  to  the  Gauss  hypergeometric 
equation  given  in  [35],  pp.  152,  Sec.  2.1.2,  Eq.  (158), 

X  {x  -  l)y”  +  [{a  +  /3  +  1)  X  - 'y]yH  aPy  =  0  (C.17) 


with  the  boundary  conditions  y{xi)  =  y{x2)  =  0  and  the  parameters 


a  =  ±i 


VIE 


.  \/l5  _ 

P  =  =  a 
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Since  7  is  not  an  integer,  the  general  solution  of  the  Gauss  hypergeometric  equation, 
Eq.  (C.17),  is  [35] 

y  =  CiF(a,  (3, 7;  x)  +  C2  x^~'^F{a  -  7  +  1,  /3  -  7  +  1, 2  -  7;  x)  (C.18) 


where 


F{a,l3,r,x) 


r(7) 

-|  I  i^)k 

k 


(1  —  t)"''  ^  ^{l  —  tx)  ^  dt 

{a)^  =  a  {a  +  1)  . . .  {a  +  k  —  1) 


is  a  hypergeometric  function  and  r(/3)  is  the  gamma  function.  Without  loss  of  generality^, 
let  a  =  and  /?  =  —i'^.  We  have 


(C.19) 


Applying  the  first  boundary  condition  y(xi)  =  0  and  solving  Eq.  (C.19)  for  the  constant 
Cl,  yields 


-  _  C2F{-\+i^,-\-i^,\;xi) 

Substituting  Eq.  (C.20)  into  Eq.  (C.19),  and  imposing  the  second  boundary  condition 
y{x2)  =  0,  yields  the  equation  in  C2 


C2 


X2)  E(— ^  —  \ 

^iF{i^,-i^,l-xi) 


=  0  (C.21) 


^We  have  two  solutions  for  the  parameters  a  and  /?•  However,  P  =  a,  the  complex  conjugate  of  a. 
Hypergeometric  functions  have  the  property  F{a,  P,^-,x)  —  F{P,a,'y,x)  -  see,  e.g.,  [31].  Thus,  we  need 
only  consider  one  of  the  two  solutions  for  a  and  /3- 
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Proposition  C.2.1.  For  any  X2  /  xi,  xi,X2  G  (0,1),  the  expression 


^(-2 +^^>-5  -^^>25^2)  F(z^,-z#,i;x2)P(-5+z#,-5  -z^,5;xi) 


is  not  zero. 


Proof.  Momentarily  let  X2  >  xi.  Define 


fi{x) 

fsix)  =  y/xF{l^^, 
=  \/x  fl{x) 


It  can  be  shown  that  the  series  representations  of  the  hypergeometric  functions  /i(x),  f2{x) 
and  fsix)  are  absolutely  convergent  -  see,  e.g.,  [31].  Furthermore,  the  functions  /i(x),  f2{x) 
and  fs^x)  are  continuous  and  monotonically  increasing  on  the  interval  (0, 1)  -  see,  e.g. 
Fig.  C.2.  Evidently,  the  functions  have  many  properties,  including 


/i(x)  >  0 
f2{x)  >  0 
h{x)  >  0 
f2{xi)  >  fsixi) 

f2{x2)  >  h{x2) 

f2{x2)  >  f2{xi) 
f2{x2)  >  h{xi) 
f3{x2)  >  fsixi) 


Define 


9{x) 


fsjx) 

f2{x) 
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60 


50 


Figure  C.2  The  Hypergeometric  Functions  fi{x),f2{x)  and  fsix) 


Clearly, 


0  <  g{x)  <  1 


Furthermore,  g{x)  is  monotonically  increasing,  see,  e.g..  Fig.  C.3,  so  that 


9{x2)  >  g{xi) 


for  all  X2  >  xi-  Thus, 


/a  (3^2)  h{xi) 

f2{x2)  f2{xi) 

h{x2)  f2{xi)  >  h{xi)  f2{x2) 


fl{x2)  f2{xi)  >  ^/^  fl{xi)  f2{x2) 


f2{x2) 


fl{x2)  f2{xi) 

fl{x2)  f2{xi) 
V^fiixi) 


f2{x2) 

<  0,  y  Xi  €  (0, 1) 


(C.22) 
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Figure  C.3  The  Hypergeometric  Function  g{x) 

Or,  equivalently,  the  expression  Eq.  (C.22)  is  not  zero.  To  show  this  holds  for  xi  >  X2,  the 
same  procedure  is  followed,  resulting  in  the  inequality  Eq.  (C.22)  being  reversed.  There  is 
no  zero  crossing  for  any  segment  (xi,  X2)  where  xi  7^  X2.  □ 

Lemma  C.2.3.  The  strengthened  Jaeobi  eondition  holds,  viz.,  that  the  interval  [60, Of], 
where  Oq  >  — 7r/2  and  Of  <  7r/2,  eontains  no  points  eonjugate  to  the  point  Oq. 

Proof.  By  Proposition  C.2.1,  Eq.  (C.22)  is  not  equal  to  zero,  which  implies  the  constant 
C2  in  Eq.  (C.21)  must  be  zero.  Hence,  from  Eq.  (C.20),  Ci  =  0  and  the  solution  to  the 
differential  equation  resulting  from  the  Jacobi  equation  is  y  =  0.  Hence,  there  are  no 
conjugate  points.  □ 

Theorem  C.2.1.  For  the  two  radar  exposure  minimization  problem,  any  segment  of  the 
Voronoi  edge  is  loeally  minimizing.  □ 

In  conclusion,  in  view  of  the  above  analysis,  the  problem  with  the  Voronoi  edge  is 
the  latter’s  ceasing  to  be  a  global  minimum  -  this,  according  to  the  condition  imposed  by 
Theorem  6.3.1. 
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